[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

[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

[R] Null values In R.

2008-12-02 Thread paul murima
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 __ R-help@r

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 did

[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 data

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

2008-11-16 Thread paul murima
(gplots) zcz- read.table(chitekete.csv, header = T, row.names=1, sep=,); xcc - cor(zcz); pdf(SIGA.pdf, height=10,width=9) heatmap.2(xcc, margins= c(9,9), col = bluered(64), trace=c(none), breaks=c(seq(-1,1,1/32)), symkey=TRUE, density.info=histogram, cexRow=1) I thank you all. Paul -- BEST Paul

[R] plot.hclust with lots of objects

2008-11-07 Thread paul murima
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 commented, minimal

[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 some