[R] concordance correlation coefficient using R

2008-12-07 Thread paul murima
Hi. I have data which i would want to assess the degree of agreement between two assays, e.g., to evaluate reproducibility or for inter-rater reliability. I have used the Pearson product-moment correlation coefficient. It looks good ranginging between 0.90 to 0.998. Though this looks good. I am t

[R] Post Hoc methods for anova in R

2008-12-03 Thread paul murima
Dear R community Does anyone know of some methods already programmed up open source in R for the post hoc methods in anova. I have some microarray data generated from different drug treatments 8 in particular, and a subset of 90 genes which i would like to analyseusing the above method. Thank y

[R] Null values In R.

2008-12-02 Thread paul murima
"); My problem is largely when i attempt to use correlation for my data... xcc <- cor(a); The error i get is as follows Error in cor(a) : missing observations in cov/cor In addition: Warning message: In cor(a) : NAs introduced by coercion Regards Paul Murima _

Re: [R] row labels in heatmap.2(package gplots)

2008-11-25 Thread paul murima
Try Exporting to a Large pdfpdf("filename.pdf", width=20, height=20) You can adjust the width and height accordingly On Wed, Nov 26, 2008 at 4:30 AM, David Winsemius <[EMAIL PROTECTED]>wrote: > I could not tell from the help file whether rotation of the x labels is > supported in heatmap.2. I di

[R] Error in image.default(1:nc, 1:nr, x, xlim = 0.5 + c(0, nc), ylim = 0.5 + : must have one more break than colour

2008-11-17 Thread paul murima
heatmap.2(m, margins= c(9,9), col = bluered(64), trace=c("none"), breaks=c(seq(-60,0,60/20)), symkey=TRUE, density.info="histogram", cexRow=1,) Error in image.default(1:nc, 1:nr, x, xlim = 0.5 + c(0, nc), ylim = 0.5 + : must have one more break than colour Hie all. I am working on scaling my d

[R] HELP ON SCALING GENE EXPRESSION DATA TO -1,0,1

2008-11-16 Thread paul murima
"histogram", cexRow=1) I thank you all. Paul -- BEST Paul Murima __ R-help@r-project.org mailing list https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read the posting guide http://www.R-project.org/posting-guide.html and provide commented, minimal, self-contained, reproducible code.

[R] plot.hclust with lots of objects

2008-11-07 Thread paul murima
t I'm unable to figure it out. Thanks in advance -- BEST Paul Murima __ R-help@r-project.org mailing list https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read the posting guide http://www.R-project.org/posting-guide.html and provide commente

[R] Clustering In R. (rookie)

2008-11-04 Thread paul murima
Hi all. I have alrge microarray dat set that i would like to analyze using hierarchical clustering. The problem is when i use the command below, > hc<- hclust(dist(array), "ave") i get get this feedback... Error in as.vector(x, mode) : cannot coerce type 'closure' to vector of type 'any' Can som

[R] R rockie

2008-10-29 Thread paul murima
Help. - Is there a beginner's manual for R? - How do i analyse gene expression data using R, to generate a dendrogram. I would greatlyy appreciate every bit of input. -- BEST Paul [[alternative HTML version deleted]] __ R-help@r-pr