We like to always keep the discussion on the list, to avoid having to
answer duplicate questions.
Collapsing all the patients into a single point defeats the purpose of
PCA: to see the distances between individual samples and groups of
samples. Showing just the mean for each group might mislead someone
looking at the plot into thinking the clusters are distinct, when the
samples might have high variance around that average point. I would
recommend instead just coloring the types of carcinoma.
On Wed, Jul 16, 2014 at 5:43 AM, deepak karthik <deepaksrna at="" gmail.com=""> wrote:
> Thanks for your reply.
> I have data of hundreds of patients from each carcinoma , consisting
> rnaseq expression with certain gene of interest. If i perform pca
> for numerous carcinoma , my pca plot would be clumsy difficult to
> the type of carcinoma are clustered together . so i would like to
> single point for a particular type of carcinoma with consideration
> rnaseq expression for my gene of my interest . Thanks in advance .
> with regards,
> On Tue, Jul 15, 2014 at 9:31 PM, Michael Love <michaelisaiahlove at="" gmail.com="">
>> You might get more feedback if you describe what kind of experiment
>> you have performed (microarray or RNA-Seq?).
>> The other reason you might not be getting response is that the
>> principal component functions are not implemented in Bioconductor,
>> in base R. So it's not necessarily a Bioconductor question, but a
>> statistics/R question.
>> The very basic code for making a PCA plot from an expression set
>> pc = prcomp( t ( exprs( e ) ) )
>> plot( pc$x[ , 1:2 ] )
>> On Tue, Jul 15, 2014 at 8:36 AM, karthik [guest] <guest at="" bioconductor.org="">
>> > Hi all,
>> > I am writing this mail for second time. I wanted perform
>> > analysis ,for each cancer type and genes of interest expression.
>> > wanted to plot only a single point which is able represent each
>> > their genes expression .Can you please explain me on it.( And
>> > gene basis should i take median or mean values to represent their
>> > expression). Thanks in advance.
>> > -- output of sessionInfo():
>> > pca()
>> > --
>> > Sent via the guest posting facility at bioconductor.org.
>> > _______________________________________________
>> > Bioconductor mailing list
>> > Bioconductor at r-project.org
>> > https://stat.ethz.ch/mailman/listinfo/bioconductor
>> > Search the archives:
>> > http://news.gmane.org/gmane.science.biology.informatics.conductor
> Deepak karthik
> PhD student
> Dr. Mali Salmon-Divon, Genomic Bioinformatics laboratory
> The Department of Molecular Biology
> Ariel University, Israel