[R] Quick question: Omitting rows and cols with certain percents of missing values

2011-05-13 Thread Vickie S
Hi naive question. It is possible to get R command for omitting rows or cols with missing values present. But if i want to omit rows or cols with i.e . >20% missing values, I could´t find any package-based command, probably because it is too simple for anyone to do that manually, though not fo

[R] Error message while mapping probes with hgu133plus2.db

2011-06-02 Thread Vickie S
Hi, eset is an expressionset based on hgu133plus2 platform. Using hgu133plus2.db package, I want to map probes to UniGene cluster IDs. It results in an error message. >xx=hgu133plus2UNIGENE >uniaf1=xx[[as.character(featureNames(eset))]] Error in .checkKeys(value, Lkeys(x), x@ifnotfound) : v

[R] Question abt power calculation in multistrata data using Hmisc or other package

2011-04-16 Thread Vickie S
Hi, I have two unequal patient groups belonging to different geographic locations, treated by different drug regimens and their survival data. N1=150 and N2=60 In each group of data there are three strata for one independent variable- treatment status (2=two drugs administered, 1=only one conv

[R] profile plot in R

2011-06-17 Thread Vickie S
Hi friends, I have a matrix with following format. group var1 var2 ...varN c1 group1 1.2399 1.4990-1.4829 c2 group4 0.8989 0.7849.1.8933 ... ... c100 group10 . I want to draw a profile plot of each condition c1 to c100, which rows in above matrix and each line represen

Re: [R] profile plot in R

2011-06-20 Thread Vickie S
ta, grp) long.dat=melt(sampledata[,1:5], idvars=as.character(rownames(sampledata)), measure.vars=as.character(colnames(sampledata)[1:5]), variable_name="cols") ggplot(long.dat, aes(x = cols, y = value)) + +  geom_line(aes(colour = sampledata$grp)) + + facet_grid(~ rownames(sampledata

Re: [R] profile plot in R

2011-06-20 Thread Vickie S
8/en_US.UTF-8 attached base packages: [1] grid  stats graphics  grDevices utils datasets  methods   base other attached packages: [1] plyr_1.1  ggplot2_0.8.9 proto_0.3-8   reshape_0.8.3 loaded via a namespace (and not attached): [1] tools_2.10.1 Any suggestions ? Vickie S ---

Re: [R] profile plot in R

2011-06-20 Thread Vickie S
Hi Dennis and all friends, After upgrading to R 2.11.1, I managed to work this out perfectly. I also added group information here, since I wanted to color code the profiles according to the groups they were assigned to. library("reshape") library("ggplot2") sampledata <- matrix(c(1.002, 1.76,

[R] Question about converting list items in matrix

2011-07-20 Thread Vickie S
Hi friends, I have got a list where each element might have variable number of members. $`4213` [1] "214077_x_at" $`164832` [1] "225996_at" "235977_at" $`339010` [1] NA $`23410` [1] "221562_s_at" "221913_at" "49327_at" $`285386` [1] "229764_at" $`2099` [1] "205225_at" "211233_x_at" "2

Re: [R] RES: Question about converting list items in matrix

2011-07-21 Thread Vickie S
dta <- c() for (i in 1:length(lst)) dta <- rbind(dta, as.matrix(lst[[i]])) dta Hope it helps. Cheers, Filipe -Mensagem original- De: r-help-boun...@r-project.org [mailto:r-help-boun...@r-project.org] Em nome de Vickie S Enviada em: quarta-feira, 20 de julho de 2011 18:04 Para: r-

[R] filterMicroRna function: Sample replicates in preprocessing Agilent miRNA dataset

2011-07-28 Thread Vickie S
option for this type of filtering. Are there any suggestions ? Thanks Vickie S [[alternative HTML version deleted]] __ R-help@r-project.org mailing list https://stat.ethz.ch/mailman/listinfo/r-help PLEASE d

[R] samr package: missing values in RNAseq data

2012-08-10 Thread Vickie S
t(prop, c(0.25, 0.75)) :   missing values and NaN's not allowed if 'na.rm' is FALSE It seems like there is an inbuilt quantile normalization function. Though i don't know how to debug this error. I also like to receive any suggestion about interpretation of missing count values in th

[R] dropterm in MANOVA for MLM objects

2012-02-08 Thread Vickie S
Dear R fans, I have got a difficult sounding problem. For fitting a linear model using continuous response and then for re-fitting the model after excluding every single variable, the following functions can be used. library(MASS) model = lm(perf ~ syct + mmin + mmax + cach + chmin + chmax, dat

Re: [R] dropterm in MANOVA for MLM objects

2012-02-08 Thread Vickie S
ology > McMaster University > Hamilton, Ontario, Canada > http://socserv.mcmaster.ca/jfox > > > > > > -Original Message- > > From: r-help-boun...@r-project.org [mailto:r-help-bounces@r- > > project.org] On Behalf Of Vickie S > > Sent: February-08-1

Re: [R] dropterm in MANOVA for MLM objects

2012-02-09 Thread Vickie S
) > > mod <- lm(cbind(A, B, C, D, E) ~ syct + mmin + mmax + cach + chmin + > chmax, data=Data) > > Anova(mod) > > Type II MANOVA Tests: Pillai test statistic > Df test stat approx F num Df den Df Pr(>F) > syct 1 0.41622 18.395 5 129 9.31e-14 *** > m

[R] Hotelling T2 test extension for multigroup data

2012-02-09 Thread Vickie S
Hi all, I've got the following matrix :   mat <- matrix(rnorm(700), ncol=5, dimnames=list( paste("f", c(1:140), sep="_"), c("A", "B", "C", "D", "E"))) I can see that currently most of the multivariate Hotelling T2 tests are limited for application on two groups/samples. I wud appreciate if so

Re: [R] Hotelling T2 test extension for multigroup data

2012-02-09 Thread Vickie S
ickie > Subject: Re: [R] Hotelling T2 test extension for multigroup data > From: pda...@gmail.com > Date: Thu, 9 Feb 2012 14:08:22 +0100 > CC: r-help@r-project.org > To: is...@live.com > > > On Feb 9, 2012, at 12:11 , Vickie

Re: [R] dropterm in MANOVA for MLM objects

2012-02-09 Thread Vickie S
't > perform statistical tests. As Peter Dalgaard pointed out, the same confusion > was reflected in your subsequent question about Hotelling's T^2. Hotelling > T^2 is equivalent to MANOVA when there are two groups. > > Best, > John > > On Thu, 9 Feb 2012 09: