Displaying 2 results from an estimated 2 matches for "reduceddata".
2017 Oct 10
1
Unbalanced data in split-plot analysis with aov()
...iduals 10?? 6013?? 601.3
Error: Within
????????????? Df Sum Sq Mean Sq F value?? Pr(>F)
nitro????????? 3? 20020??? 6674? 37.686 2.46e-12 ***
Variety:nitro? 6??? 322????? 54?? 0.303??? 0.932
Residuals???? 45?? 7969???? 177
---
Signif. codes:? 0 ?***? 0.001 ?**? 0.01 ?*? 0.05 ?.? 0.1 ? ? 1
> reduceddata <- data[-5,] # delete one observation
> summary(aov(yield ~ Block+Variety*nitro +
Error(Block/Variety),reduceddata))
Error: Block
????? Df Sum Sq Mean Sq
Block? 5? 16070??? 3214
Error: Block:Variety
????????? Df Sum Sq Mean Sq F value Pr(>F)
Variety??? 2?? 1820?? 910.0?? 1.373? 0.302...
2007 Dec 08
0
help for segmented package
...imated breakpoints are out of range.
Since my problem is time related problem. So I hink the
breakpoints should be all at the same time for all variables. Can I
achieve this requiement with this package? How?
Thanks
Best Regards
Peter
Here is the code I used and the data is also attached.
reduceddata <- read.table("reduceddata.txt")
detach()
attach(reduceddata)
library(segmented)
dat=data.frame(y=V2,x1=V3,x2=V4,x3=V5,
x4=V6,x5=V7,x6=V8,x7=V9,x8=V10,x9=V11,
x10=V12,x11=V13,x12=V14,x13=V15,x14=V16,x15=V16)
m=2
psi <- array(0, c(15,m))
for (i in 1:m)
{
psi[1,i]=range(...