similar to: Difficult to set a quiet formula in maanova

Hi, I run a maanova analysis and found this message error: Error in ma.svd(X, 0, 0) : 0 extent dimensions I did a google search and found this: \item ma.svd: function to compute the sigular-value decomposition of a rectangular matrix by using LAPACK routines DEGSVD AND ZGESVD. \item fdr: function to calculate the adjusted P values for FDR control. I did a search for LAPACK and

Dear R-users: I am using the package MAANOVA to analyze microarray data and have encountered problems when trying to get interactions. I am a newbie in both, R and maanova, and I do not have good knowledge in statistical methods. I have four effects: Effects Levels Var 2 Ind 2 Trat 2 Time 4 Sample 3 <-- biological replicate Spot 4 <-- technical replicate I had

Dear R-users: I am using the package MAANOVA to analyze microarray data and have encountered problems when trying to plot data. I have tried emailing a MAANOVA discussion group, as well as the author of the package, and have not yet received a response so I am hoping that someone on this listserv can be of assistance. There are several functions in MAANOVA (riplot, resiplot) which call the

"c" an extremely commonly-used function. Functions are first-class objects that occupy the same namespaces that variables do, so they can obscure each other. In short, don't use variables called "c" (R is case sensitive, so "C" has no such problem). Wherever possible, avoid incremental concatenation like the plague. If you feel you must use it, at least

Maybe I wasn't clear about my query. I'm very familiar with pre-allocation and vectorization and I had already wrote an R code for this problem in this way. My question wasn't about the most efficient way to solve the problem. It was about whether in R it was possible to use the same index used in the loop to create a new variable and store the results in as in the example showed

... "(R is case sensitive, so "C" has no such problem)." Well, not quite. Try ?C To add to the previous comments, Dr. Gordon appears to need to do her/his homework and spend some time with an R tutorial or two before posting further here. There are many good ones on the web. Some recommendations can be found here: https://www.rstudio.com/online-learning/#R Cheers, Bert

Dear professors and collegues I need to perform a analysis of dates from a nested experimental design. From "Bioestatical Analysis" of Zar "Bimetry of Sokal" & Rohlf "Design and Analysis of Experiments" of Montgomery I have: Sum (mean(x)_i - mean(x)_T)2 / (a-1) -> var(epsilon) + n sigma2_B + n b (sum alfa_i)2 / (a-1) Sum (mean(x)_ij - mean(x)_i)2 /

Hello, Suppose the you need a loop to create a new variable , i.e., you are not reading data from outside the loop. This is a simple example in Matlab code, for i=1:5 r1=randn r2=randn r=[r1 r2] c(i,:)=r; % creation of each row of c , % the ":" symbol indicates all columns. In R this would be [i,] end The output of interest is c which I'm creating inside the "for" loop

I'm trying to install packages on windows XP and I have trouble with command Rcmd build (R version 1.8.1) : In the Windows console for package maanova for example, answer is : C:\Documents and Settings\dillies\Mes documents\ghis\packages>Rcmd build maanova * checking for file 'maanova/DESCRIPTION' ... OK * preparing 'maanova': * cleaning src *

Hi, I'm building Tukey's table using qtukey function. It happens that I can't get the values of Tukey's one degree of freedom and also wanted to eliminate the first column. The program is: Trat <- c(1:30) # number of treatments gl <- c(1:30, 40, 60, 120) # degree freedom tukval <- matrix(0, nr=length(gl), nc=length(Trat)) for(i in

Hi all, I'm working with a dataset from 10 treatments, each treatment with 30 subjects, each subject measured 5 times. The plot of the dataset suggests that a 3-parameter logistic could be a reasonable function to describe the data. When I try to fit the model using gnls I got the message 'Step halving factor reduced below minimum in NLS step'. I´m using as the initial values of the

Hi, I used the commands below to make Hartley's table, but some values are NA. require(SuppDists) trat = seq(2, 15, 1) gl = seq(2, 40, 1) har = matrix(0, nr=length(gl), nc=length(trat)) for(i in 1:length(gl)) for(j in 1:length(trat)) har[i,j] <- qmaxFratio(.95, df=gl[i], k=trat[j]) rownames(har) <- gl colnames(har) <- trat head(har) The output (head): 2

Oct 16 22:52:42 mta0 dovecot: lda(luciano.moreira at metal.eeimvr.uff.br): sieve: msgid=<a04f508ac9d299a0733a09ce714ae56f at localhost.localdomain>: stored mail into mailbox 'luciano/moreira at metal/eeimvr/uff/br' Oct 16 22:52:48 mta0 dovecot: lda(fabiana at vm.uff.br): sieve: msgid=< 4CB3CFE3.4060206 at wspc.com.sg>: stored mail into mailbox 'fabiana at vm/uff/br' Oct

I'm new using R im trying to do a tukey test, but when i see the results the p value results in NA im guessing its because i have missing values but im not sure how to fix it AnovaModel.2 <- aov(area ~ trat, data=apilados) > summary(AnovaModel.2) Df Sum Sq Mean Sq F value Pr(>F) trat 4 11847 2961.76 9.9905 1.500e-06 *** Residuals 76 22531 296.46

Hi, I try to plot bars with different colours in a barchart graphic. My idea is make that all X-Levels from trat var with different colour (grey scale). I search for a solution but dont find any. Any help? Thanks dados <- structure(list(Medias = c(0.994169096209855, 0.99416342412449, 0.974683544303797, 0.954430379746835, 0.669047619047619, 0.999999998475569, 0.994163424124514,

Hi!! Can anyone help me, i have problems to converge the following data with 5 nonlinears models that i evaluated. Firtly, i send my data (totalsinatipicos) that i just try to fit with the nonlinear models. Next, i have the following script where i called the data as totalsinatipicos. I made selfstarting each nonlinear model. ###Library library(NRAIA) ###Data d<-totalsinatipicos

Hi, I did the analysis of variance of a split-plot and the effect of treatment was significant. I would like compare treatment means using Tukey. I can't extract the mean square to apply HSD.test to use in agricolae package. anava = aov(ganhos ~ Blocos + Trat*Supl + Error(Blocos/Trat)) names(anava) summary(anava) require(agricolae) HSD.test(ganhos, Trat, df, MSerror, alpha = 0.05)

Hi, I am using "*Maanova* package" to do anova. I have created *datafile* with probeID as the first column, which is a tab limited text file and also created *designfile*. I have created *readma object* which is named as abf1. >From that readma object, i have to create data object by using *createData*function and also i hav to create model object by using *makemodel* function,

Hi, I am using "*Maanova* package" to do anova. I have created *datafile* with probeID as the first column, which is a tab limited text file and also created *designfile*. I have created *readma object* which is named as abf1. >From that readma object, i have to create data object by using *createData*function and also i hav to create model object by using *makemodel* function,

Hi R-users I'm trying to do multivariate analysis of variance of a experiment with 3 treatments, 2 variables and 5 replicates. The procedure adopted in SAS is as follow, but I'm having difficulty in to implement the contrasts for comparison of all treatments in R. I have already read manuals and other materials about manova in R, but nothing about specific contrasts were found in them,