File:RbdBlock.png
Jump to navigation
Jump to search
RbdBlock.png (500 × 500 pixels, file size: 34 KB, MIME type: image/png)
Summary
set.seed(1234) block <- as.factor(rep(1:5, each=6)) treatment <- rep(c("A","B","C"),5) block_shift <- rnorm(5, mean = 0, sd = 2) treatment_shift <- c(A=0, B=4, C=2) random_effect <- rnorm(30, mean = 0, sd = 1) yield <- rnorm(30, mean = 10, sd = 2) + treatment_shift[as.integer(factor(treatment))] + block_shift[as.numeric(block)] + random_effect data <- data.frame(block, treatment, yield) summary(fm1 <- aov(yield ~ treatment + block, data = data)) # Df Sum Sq Mean Sq F value Pr(>F) # treatment 2 92.09 46.05 9.555 0.000954 *** # block 4 140.45 35.11 7.286 0.000607 *** # Residuals 23 110.84 4.82 interaction.plot(x.factor = data$block, trace.factor = data$treatment, response = data$yield, fun = median, col=c("green","red", "blue"), xlab="Block", ylab="Yield", trace.label="Treatment", lty=1:3, lwd=2.5 )
File history
Click on a date/time to view the file as it appeared at that time.
Date/Time | Thumbnail | Dimensions | User | Comment | |
---|---|---|---|---|---|
current | 21:51, 14 January 2023 | 500 × 500 (34 KB) | Brb (talk | contribs) | <pre> set.seed(1234) block <- as.factor(rep(1:5, each=6)) treatment <- rep(c("A","B","C"),5) block_shift <- rnorm(5, mean = 0, sd = 2) treatment_shift <- c(A=0, B=4, C=2) random_effect <- rnorm(30, mean = 0, sd = 1) yield <- rnorm(30, mean = 10, sd = 2) + treatment_shift[as.integer(factor(treatment))] + block_shift[as.numeric(block)] + random_effect data <- data.frame(block, treatment, yield) summary(fm1 <- aov(yield ~ treatment + block, data = data)) # Df Sum Sq Mean Sq... |
You cannot overwrite this file.
File usage
The following page uses this file: