similar to: nested aov: plot available?

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, I posted this question several days ago, but did not get any answer until now. Since I still have no clue about the source of this error message, I repost a description of the problem including some code: A student at our institute fitted an aov model, and got the following error message: > plot(p.aov) Hit <Return> to see next plot: Hit <Return> to see next plot:

Dear WizaRds, having met some of you in person in Vienna, I think even more fondly of this community and hope to continue on this route. It was great talking with you and learning from you. Thank you. I am trying to work through an artificial example in post stratification. This is my dataset: library(survey) age <- data.frame(id=1:8, stratum=rep(

Dear WizaRds, Working through sampling theory, I tried to comprehend the concept of stratification and apply it with Survey to a small example. My question is more of theoretic nature, so I apologize if this does not fully fit this board's intention, but I have come to a complete stop in my efforts and need an expert to help me along. Please help: age<-matrix(c(rep(1,5), rep(2,3),

Dear R-help, I'm facing a problem with defining a repeated measures anova with weighted data. Here's the code to reproduce the problem: # generate some data seed=11 rtrep <- data.frame(rt=rnorm(100),ti=rep(1:5,20),subj=gl (20,5,100),we=runif(100)) # model with within factor for subjects/repeated measurements, no problem aov(rt~ti + Error(subj/ti),data=rtrep) #model with weights

I believe you are right, but can you please explain why anyone would want to fit this model? It differs only in the coding from aov(y ~ a + b + Error(c), data=test.df) and merely lumps together the top two strata. There is a much simpler fix: in the line if(intercept) nmstrata <- c("(Intercept)", nmstrata) remove the condition (and drop the empty stratum later if you

Dear all, We are analysing some survey data and we are not sure if we are using the correct syntax for our design. The population of interest is a set of 4416 polygons with different sizes ranging from 0.003 to 45.6 ha, 7460 ha in total. Each polygon has a binary attribute (presence/absence) and we want to estimate the probability of presence in the population. We used sampling with replacement

Dear list, I got some surprising results when using the svytotal routine from the survey package with data containing missing values. Some example code demonstrating the behaviour is included below. I have a stratified sampling design where I want to estimate the total income. In some strata some of the incomes are missing. I want to ignore these missing incomes. I would have expected that

At any rate: Error(SUBJECT/IV) specifies two random effects: SUBJECT and SUBJECT:IV. This is most easily understood if you conceptually arrange your data in a SUBJECT x IV table: One effect is a set of random errors added to each row, the other is a set of effects added to each cell. If you have more than one observation within each cell, then you need a third set of errors to account for

Hi, I am playing with clogit(case~spontaneous+induced+strata(stratum),data=infert) from clogit help file. This line works. 1. But, why strata(stratum) doesn't have a coefficient like spontaneous and induced? 2. When I remove strata(stratum) from the command, this function seems to keep running forever. Why? 3. I think the equation for clogit looks like P=1/(1+

I think there are two bugs in aov() that shows up when the right hand side of `formula' contains both `-1' and an Error() term, e.g., aov(y ~ a + b - 1 + Error(c), ...). Without `-1' or `Error()' there is no problem. I've included and example, and the source of aov() with suggested fixes below. The first bug (labeled BUG 1 below) creates an extra, empty stratum inside

It is unclear to me how aov() handles non-categorical variables. I mean it works and produces results that I would expect, but I was under impression that ANOVA is only defined for categorical variables. In addition, help(aov) says that it "call to 'lm' for each stratum", which I presume means that it calls to lm() for every group of the categorical variable, however I

Dear R users: I am not entirely convinced that clogit gives me the correct result when I use pspline() and maybe you could help correct me here. When I add a constant to my covariate I expect only the intercept to change, but not the coefficients. This is true (in clogit) when I assume a linear in the logit model, but the same does not happen when I use pspline(). If I did something similar

Dear R-helpers, I would like to generate a binary random variable within a stratum's stratum. Here is a simple example. ## x is the first level strata index, here I have 3 strata. x=c(rep(1,5), rep(2,5), rep(3,5)) ## within x, there is a second strata indexed by t=0 and t=1 t=rep(c(0,0,1,1,1),3) ## and within strata i and t=0 and t=1, I generate the random binomial variable respectively,

Hello, I am analyzing the results of a survey of the student body on the use of certain technologies. The student body was divided in strata according to two criteria, college and user class. The resulting 24 strata were sampled. Subjects filled out a survey with a number of questions, most of them of the yes-no variety. I created a data.frame with the results for the survey. One row for

Greetings, The aov() function may mislabel the random effects as in the example below: Has anybody else noticed this? Cordially, Giles Crane, MPH, ASA, NJPHA gilescrane at verizon.net > m2 Call: aov(formula = y ~ ap + pe + Error(ju), data = d) Grand Mean: 77.50667 Stratum 1: ju Terms: ap Sum of Squares 4322.538 Deg. of Freedom 12 13 out of 25 effects not

Hello! I have something like this: test1 <- data.frame(intx=c(4,3,1,1,2,2,3), status=c(1,1,1,0,1,1,0), x1=c(0,2,1,1,1,0,0), x2=c(1,1,0,0,2,2,0), sex=c(0,0,0,0,1,1,1)) and I can easily fit a cox model: library(survival) coxph(Surv(intx,status) ~ x1 + x2 + strata(sex),test1) However, I want to

Hello, Are there any access functions to the various residual variables that should result from a repeated measures ANOVA ? MyAOVObject$residuals does not exist, and simply printing MyAOVObject gives a very long print of all fields in the result list, many of which I can't see what they are exactly : $error.qr$qraux, for instance. What I would like basically is to inspect those residuals

Hi, I’m new to R (just installed today) and I’m trying to figure out how to do stratified randomisation using it. My google search expedition has lead me to believe that blockrand package will most probably be the answer to it. I’ve played around with blockrand for awhile and tried the sample code: library(blockrand) ##stratified by sex male <- blockrand(n=100,

Hi all, Given a number of points (x,y) in a plane, I'd like to plot a map of polygons, so that 1) each polygon contains exactly one point 2) the polygon defines the area for which this specific point is closer than any other point. It's a bit like a map of areas "influenced" by that point, and it's obviously a matter of intersecting the perpendicular bisectors