similar to: anova error

Full_Name: jean coursol Version: 1.7.1, 1.8.0 OS: linux & Windows-XP Submission from: (NULL) (129.175.52.7) Binary MS-Windows akima module from CRAN (1.8.0 version) produces wrong results with some data. Installing akima source in linux, with same data: -with gcc-2.95.3 -O2 : give correct results (under R 1.7.1); -with gcc-3.2.3 -O2 : give wrong results (under R-1.7.1 and R-1.8.0); -with

Hi, I need a help. I am new in R and I need to run a nested anova with fixed and random factors (Mixed Model). I have a design with three factors: Day, Area and Plot and the dependent variable is density. The factors Day and Area are fixed while Plot is random, factor Area is nested in factor Day, and factor Plot is nested in Area. I can do it using aov by: mod1<-aov(density~ day +

I want to do a zoom with recordPlot(). I have problems with lists. (R-2.0.1 patched 2004-11-30 , various linux). I have problems with RecordPlot class structure. > plot(1:10) > saveP <- recordPlot() > dev.off() > sx <- saveP[[1]][[2]][[2]] > saveP[[1]][[2]][[2]] <- sx Error in "[[<-"(`*tmp*`, 1, value = list(list( .Primitive("plot.new")),

Hi, I'm try to use the SJava package. The install is OK. In R I have this error message: -------------------------- > library(SJava) Error in dyn.load(x, as.logical(local), as.logical(now)) : unable to load shared library "/opt/lib/R/site-library/SJava/libs/SJava.so": libRSNativeJava.so: cannot open shared object file: No such file or directory Error in

I currently work on a draft of an aquatic bioassessment. The conditions tested are the following: ER river water T dechlorinated water control 0.5 + 0.5mg / L of malate T + 1 dechlorinated water control + 1g / L of malate T ED dechlorinated water control SED + ER + river water sediment SED ED + sediment + water dechlorinated. It is the result of AChE in muscle (fillet of fish). The production of

Hi, I searched for this in the mailing list, but found no results. I have a large dataframe ( dim(mydata)= 1297059 16, object.size(mydata= 145280576) ) , and I want to perform some calculations which can be done by a factor's levels, say, mydata$myfactor. So what I want is to split this dataframe into nlevels(mydata$myfactor) = 80 levels. But I must do this efficiently, that is, I

Hello colleagues, I am a novice with R and am stuck with an analysis I am trying to conduct. Any suggestions or feedback would be very much appreciated. I am analyzing a data set of psi (ESP) ganzfeld trials. The response variable is binary (correct/incorrect), with a 25% base rate. I've looked around the documentation and other online resources and cannot find how I can correct for

I made an ANOVA with repeated mesures (aov(Mesure~Distance*Genre*Correct+Error(Sujet/(Distance*Genre*Correct)), data)) and I would like to extract the pvalues. The output is: ----------------------------------------------------------- Error: Sujet Df Sum Sq Mean Sq F value Pr(>F) Residuals 21 97.082 4.623 Error: Sujet:Distance Df Sum Sq Mean Sq F value Pr(>F) Distance

Hello! I want to create a messagebox whenever the user wants to destroy the window (e.g. <Alt-F4> or the 'x' in the right top corner) and ask if a modified file should be saved or not. If 'cancel' is chosen then nothing should happen and the windows still should be existing. This doesn't work. When I press cancel the window will be destroyed although. I also

Hi, there I'm trying to get the value of the Mean Square from the ANOVA model summary that comes from specifying the error term, and am wondering if one can actually do this ( I know it's possible when using anova(lm) objects and the like, but I'm having a tough time with it under this framework). There does appear to be some indexing in the output of this type, but perhaps not

Hi all, If we wanted to study the effect on the mean of the hourly data based on the hours within a day... and we wanted to do Anova analysis... We have two choices: Please see below: Why are these two approaches giving very different p-values? And which one shall I use? Thanks a lot! 1. treating the hours as double/floating numbers: anova(lm(hourlydata~as.double(hours_factors))) Df Sum

Hi all, I'd like to extract a value from an ANOVA table, but experience the following problem: ### This works: > s.pseudo <- summary(aov(m ~ block + mix*graz,data=split1)) > s.pseudo Df Sum Sq Mean Sq F value Pr(>F) block 2 1114.66 557.33 4.4296 0.04192 * mix 1 6.14 6.14 0.0488 0.82956 graz 2 1.45 0.72 0.0057 0.99427 mix:graz

Hi all, Anyone experienced in the LaTeX format? I'm trying to use the xtable package to create nice anova tables, but how do I do to produce a pdf from the resulting LaTeX table? I've tried WinShell and MiKTeX, but I couldn't get any of them working... Here's an example of the output in R: % latex table generated in R 2.9.2 by xtable 1.5-6 package % Tue Nov 24

Hello all, A really noddy question for you all: I''m trying without success to do some subhypothesis testing. Using simple anova model, with a toy dataset from a book. I have four factors A,B,C,D, and wish to test mu_C = mu_D. This is what I have tried: > contrasts(infants$group,how.many=1) <- c(0,0,1,-1) > contrasts(infants$group) [,1] A 0 B 0 C 1

I wonder if someone could explain the behavior of the anova() and lm() functions in the following situation: I have a standard 3x2 factorial design, factorA has 3 levels, factorB has 2 levels, they are fully crossed. I have a dependent variable DV. Of course I can do the following to get the usual anova table: > anova(lm(DV~factorA+factorB+factorA:factorB)) Analysis of Variance Table

Hi everyone, I'm witnessing an odd modelling phenomenon that I can't explain. If anyone has seen this before, or can explain what's going on would let me know, I'd be very grateful! Especially if I'm just being dim. I'm fitting a pair of continuous variates and their interaction to some residuals from another model. The sequential anova statement changes with the term

Hello dear R members. I have been learning the Anova syntax in order to perform an SS type III Anova with repeated measures designs (thank you Prof. John Fox!) And another question came up: where/what are the (between/within) residuals for my model? ############ Play code: phase <- factor(rep(c("pretest", "posttest", "followup"), c(5, 5, 5)),

hello all, Quick question, how do I get the p value out of the anova? Thanks, Paul > pb<-aov(as.numeric(diff[5,16:33]) ~ grF) > summary(pb) Df Sum Sq Mean Sq F value Pr(>F) grF 3 2.7860e+10 9.2867e+09 4.2236 0.02534 * Residuals 14 3.0783e+10 2.1988e+09 --- Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 '

Hi all! I have 4 treatments on 5 animals Treat1 Treat2 Treat3 Treat4 Animal1 36 37 35 39 Animal2 33 34 36 37 Animal3 37 35 33 38 Animal4 34 36 34 35 Animal5 35 36 33 36 I use an Anova and i try to verify calcul So i retrieve: DF SS

Hey! Could you please take a quick look at what I have done? Somehow I get wrong results using the anova(lm()) combination compared to doing a two way ANOVA by hand. Running: Data<-read.table("Data.txt"); g<-lm(ExM~S1*S2,Data); anova(g); Gives: Analysis of Variance Table Response: ExM Df Sum Sq Mean Sq F value Pr(>F) S1 1 4.3679