search for: fitmaanova

...dTGR.Untreated dTGR were divided in two subgroups, rats with cardiac hypertrophy (n=18) and rats with THF (n=8).Rats were killed at age 7 weeks.Total RNAs (five animals per group) were isolated from the left ventricle.For all rgu34_a affymetrix chips were used. Iam using maanova package. In using fitmaanova function, we have to set variables to formula and random.Plz, suggest me what to take in formula and random. *fit.fix <- fitmaanova(rat, formula = ~Array+sample, random = ~Array)* Error in fitmaanova(rat, formula = ~Array + sample, random = ~Array) : Not enough degree of freedom to do regress...

...0' then read the data again 2: In read.madata("rmaexpr.dat", designfile = "design.dat") : Assume that intensity value is saved from the second column. Otherwise provide 'intensity' (first column storing intensity) information, and read the data again > fit.fix=fitmaanova(TGR,formula=sample,random=~1) Error in x$terms : object of type 'closure' is not subsettable > fit.fix=fitmaanova("TGR",formula=sample,random=~1) Error in fitmaanova("TGR", formula = sample, random = ~1) : The first input variable is not an object of class madata....

Dear R helpers, I am working with R 2.4.1 GUI 1.18 (4038) for MacOSX. I have a matrix of 10 000 genes and try to run the following commands: > model.mix<-makeModel (data=data, formula=~Dye+Array+Sample+Time, random=~Array+Sample) > anova.mix<-fitmaanova (data, model.mix) > test.mix<-matest (data, model=model.mix, term="Time", n.perm=100, test.method=c(1,0,1,1)) I get the following error message: Doing F-test on observed data ... Doing permutation. This may take a long time ... Error in La.svd(X) : error code 1 from Lapack routin...

...tion to calculate the adjusted P values for FDR control. I did a search for LAPACK and not found a package. Could you help me on how I could solve this problem? I am try to do this: library(maanova) read the data file and design ##### Make the full model based on the design anova.full.mix <- fitmaanova(data, formula=~ Var+ind+Trat+Time+ Var:ind+Var:Trat+Var:Time+ ind:Trat+ind:Time+ Trat:Time+ Var:ind:Trat+ Var:ind:Time+ Sample,random=~Sample) If I remove the Var:ind:Trat and Var:ind:Time from formula the script runs very well. Did is not possible to do that interactions (Var:ind:Trat and Var:in...

...Dimensions Observations(per gene): 36 Columns in X: 24 Columns in Z: 3 Warning message: number of rows of result is not a multiple of vector length (arg 2) in: cbind(1, Class, Levels, Effect) (I do no what this warning means. my be the error was here.) > anova.full.mix <- fitmaanova(fabiana, model.full.mix) Calculating variance components for fixed model... Fitting mixed effect model... Finish gene number 100 ... (...) Finish gene number 3500 ... Error in next.fix:(next.fix + ncols - 1) : NA/NaN argument I was inspecting the files of data and there's nothing wrong...

Hi, Thanks for reading. I am running a process in R for microarray data analysis. RedHat Enterprise Linux 4, dual AMD CPU, 6G memory. However, the R process use only a total of <200M memory. And the CPU usage is total to ~110% for two. The program takes at least 2 weeks to run at the current speed. Is there some way we can increase the usage of CPUs and memories and speed up? Any

Hi, I have some problem to install RMySQL package. Could you help me? I set up the path. [root at ep2 library]# export PKG_CPPFLAGS="-I/var/lib/mysql" [root at ep2 library]# export PKG_LIBS="-L/var/lib/mysql - lmysqlclient" [root at ep2 library]# R CMD INSTALL /home/younko/RMySQL_0.5- 5.tar.gz I got this bunch of error message. Please help me. * Installing *source*

...<AANLkTi=mQ5bhntKQbfUapHpWNWpVLn+[158]Hzat_FD6jpZF6 at mail.gmail.com> Content-Type: text/plain Hi, I am using maanova package for analysis. In my dataset, two fixed factors time and treatment and sample as random factor is there. Am able to get madata object and fitmaanova object. But, am unable to do f-test with two factors, but i have done f-test seperately for two factors. fit.full.mix <- fitmaanova(madata, formula = ~Sample+Time+Treatment, random = ~Sample) ftest.all = *matest*(madata, fit.full.mix, test.method=c(1, 1), shuffle.method=...