
Hello all,

 

I have the follow type of data:

 

 

farm <- factor(rep(1:2, each = 30))  ## two different locations
age <- factor(rep(1:3, each = 10, times = 2)) ## age of the secondary forests (1=young,2=intermediate,3=old)
block <- factor(rep(1:5, each = 2, times = 6)) ## experiemental blocks
trt <- factor(as.vector(replicate(30, sample(c(0, 1)))))     # treatments with randomizes order within pair
survival <- rbinom(60, 0:2, 0.5) ## survival of the two seedlings planted in the plots
death <- 2 - survival


And I want to fit a split-plot model, having age as the main-plot effect, block and trt as the split-plot effects, with all effects nested with farms.

 

Anyone could  help me on that?

 

 

Guilherme
 		 	   		  
	[[alternative HTML version deleted]]


