#! /usr/bin/perl

use strict;
use warnings;

my $pcaRScript = "pca.R";

if(@ARGV != 2) {
	print "USAGE::: perl $0 <gene expression file with annotation> <PCA plot file name>\n";
	exit;
}

my ($f, $pcaOut) = @ARGV;

# Perform PCA analysis

`Rscript --vanilla $pcaRScript $f $pcaOut`;

