similar to: aov for a split plot design

I have used lme() on data from a between-within subjects experiment. The correct ANOVA table is known because this is a textbook example (Experimental Design by Roger Kirk Chapter 12: Split-Plot Factorial Design). The lme() F-values differ from the known results. Please help me understand why. d<-read.table("kirkspf2.dat",header=TRUE) for(j in 1:4) d[,j] <- factor(d[,j]) ### Make

##I want to remove the rows where the row sums are zero and this is as far as I have gotten ffg <- (structure(list(CD = c(0, 0, 0, 0, 3.125, 0, 0, 0, 0, 1.6, 3.125, 0, 0, 6.25, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 3.125, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1.6, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1.6, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 3.125, 0, 0, 0, 0, 0, 0, 0, 0,

Hi all Consider the classic data below from Darwin on the heights of 15 pairs of zea mays (corn) plants either cross-fertilized or self-fertilized, where the goal is to see if it makes a difference. > head(ZeaMays) pair pot cross self diff 1 1 1 23.500 17.375 6.125 2 2 1 12.000 20.375 -8.375 3 3 1 21.000 20.000 1.000 4 4 2 22.000 20.000 2.000 5 5 2 19.125

station month bas 190 5 0.000 190 7 1.563 190 10 0.000 190 11 0.000 202 4 18.750 202 5 18.750 202 7 6.250 202 10 4.800 202 11 3.125 198 4 18.750 198 5 31.250 198 7 3.125 198 10 3.200 198 11 12.500 205 4 0.000 205 5 0.000 205 7 0.000

Hi, Using package gclus in R, I have created some graphs that show the trends within subgroups of data and correlations among 9 variables (v1-v9). Being interested for more details on these data I have produced also the var-covar matrices. Question: From a pair of two subsets of data (with 9 variables each, I have two var-covar matrices for each subgroup, that differ for a treatment on one

#this is my stab at - I am sure that I am missing something. If this doesn't make sense then please ask for more details. #This may show my low level of programing knowledge hester. <- c(1,1,1,1,2,2,2,2,3,3,3,3,4,4,4,4) value <- rnorm(16) x <- data.frame(value, hester.) z <- (if(x[,"hester."]==1){ x[,"value"]*6.250} else if(x[,"hester."]==2){

Dear Rexperts, I am trying to add a '+' identifying the mean in a boxplot using the following sizelist <- split(size, grp) centers <- boxplot(sizelist, style.bxp = "att", medpch = "o", ylab = "Prostate Volume (cm3)") points(centers, unlist(lapply(sizelist, mean)), pch = "+") But, I get error Error in xy.coords(x, y) :

someone sent in a question earlier about doing something in 3D so i took a stab at it purely for educational purposes ( i'm not even sure that I understood the question actually ). Unfortunately, persp gives me an error that I don't understand because it says "object y not found". I'm sending y in as a parameter to persp similar to what ?persp shows in one of oits examples

Dear all, I am stuck with a problem that might be trivial for most of you (and therefore is a bit embarrassing for me...): I want to calculate a generalized least squares regression using two time series (Y depending on X) with an autoregressive correlation structure of order two (the data along time are given below). I use 'gls' from package 'nlme': Calib.gls <- gls(Y~X,

Hi, I have a situation here. I try this update: mmaa <- update(mma,biomass~qvartemp) but I have this message: Error in eval(expr, envir, enclos) : object "qvartemp" not found but this object exist: [1] "cont" "i" "levelsord" "mma" "qvar" "qvarmma" [7] "qvartemp" "test"

On Wed, Feb 6, 2019 at 3:28 PM Dan Espen <dan1espen at gmail.com> wrote: > > Ilia Mirkin <imirkin at alum.mit.edu> writes: > > > It would be useful to know how the screen is connected. Also please > > grab the monitor's EDID from /sys/class/drm/cardN-connector/edid and > > attach it here. It would also be interesting to get a boot with > >

Ilia Mirkin <imirkin at alum.mit.edu> writes: > On Wed, Feb 6, 2019 at 3:28 PM Dan Espen <dan1espen at gmail.com> wrote: >> >> Ilia Mirkin <imirkin at alum.mit.edu> writes: >> >> > It would be useful to know how the screen is connected. Also please >> > grab the monitor's EDID from /sys/class/drm/cardN-connector/edid and >> >

Hello R-users, I apologize for my question but I'm a newbie. I want to read a file which columns separator is "\t". At the end of each row there is a "\n" to go to the following line. The three first lines are remarks lines and the fourth contains columns titles ( variables names, ids, dates, calculated values, observed values...) . I do: read.table("myFile",

Dan Espen <dan1espen at gmail.com> writes: > Ilia Mirkin <imirkin at alum.mit.edu> writes: > >> On Wed, Feb 6, 2019 at 3:28 PM Dan Espen <dan1espen at gmail.com> wrote: >>> >>> Ilia Mirkin <imirkin at alum.mit.edu> writes: >>> >>> > It would be useful to know how the screen is connected. Also please >>> > grab

On Sat, Feb 16, 2019 at 12:49 PM Dan Espen <dan1espen at gmail.com> wrote: > > Dan Espen <dan1espen at gmail.com> writes: > > > Ilia Mirkin <imirkin at alum.mit.edu> writes: > > > >> On Wed, Feb 6, 2019 at 3:28 PM Dan Espen <dan1espen at gmail.com> wrote: > >>> > >>> Ilia Mirkin <imirkin at alum.mit.edu> writes:

Dear R fans, I would like to create a scatterplot showing the relationship between 4 continuous variables. I thought of using the package "scatterplot 3d" to have a 3-dimensional plot and then using the size of the symbols to represent the 4th variable. Does anybody know how to do this? I already tried to create this graph using the colour of the symbols, but I was unable to generate

Based on what follows, the most favourable construction is that the documentation, by failing to say that the function should be used only for completely balanced designs, is deficient. Even for such designs, there is a bug that needs correction. The discussion is lengthy because I think it important to document, at least wrt effects and means, exactly what model.tables() does. My preference is

-----BEGIN PGP SIGNED MESSAGE----- Hash: SHA1 I'm at a loss here. Details of my setup: OS: Solaris 10 gcc version: 3.4.3 postfix version 2.5.2 dovecot version 1.1.2 dovecot sieve version 1.1.5 How I installed: Postfix: make makefiles CCARGS='-DUSE_SASL_AUTH \-DDEF_SERVER_SASL_TYPE=\"dovecot\" -DUSE_TLS -I/opt/csw/include' \AUXLIBS="-R/opt/csw/lib -L/opt/csw/lib

read.table("bug.txt", header=TRUE) plotMeans(bug$tot, bug$station, error.bars="se", main="Total") The plot is ordered in this way HC RM190 RM198 RM202. I would like it to plot in this way HC RM202 RM198 RM190. thanks Stephen Is this okay or would you like the data to be copied and pasted into the message body. I am trying to optimize my postings. -- Let's

Hi, I am trying to find m breakpoints of a linear regression model. I used the segmented package. It works fine for small number of predicators and breakpoints.(3 r.v. 3 points). However, my model has 14 variables it even would not work even for just one breakpoints!. The error message is always estimated breakpoints are out of range. Since my problem is time related problem. So I