Displaying 2 results from an estimated 2 matches for "ssprod".
Did you mean:
siprod
2010 Jun 10
1
do faster ANOVAS
...1 A1 S1 R1 1 0 0 … 1
I want to realize 1 ANOVA per timepoint and per attribute, there are 101 timepoints and 8 attributes so I want to realize 808 ANOVAS. This will be an ANOVA with two factors :
Here is one example:
Aov(t1~Subject*Product,data[data$attribute==”A1”,])
I want to store for each ANOVA SSprod,SSsujet,SSerreur,SSinter and SStotal.
In fact I want the result in several matrices:
Ssprod matrice:
T1 t2 t3 t4 … t101
A1 ssprod(A1,T1)
A2
A3
…
A8
So I would like a matrice like that for ssprod, ssujet,sserreur,ssinter and sstotal.
And this is for one permutation, and I want to do 1000 permutation...
2003 Jul 30
2
Comparing two regression slopes
...e <-function(lm1, lm2) {
## lm1 and lm2 are two linear models on independent data sets
coef1 <-summary(lm1)$coef
coef2 <-summary(lm2)$coef
sigma <-(sum(lm1$residuals^2)+sum(lm2$residuals^2))/(lm1$df.residual +
lm2$df.residual-4)
SSall <-sum(lm1$model[,2]^2) + sum(lm2$model[,2]^2)
SSprod <-sum(lm1$model[,2]^2) * sum(lm2$model[,2]^2)
F.val <-(as.numeric(coefficients(lm1)[2]) - as.numeric(coefficients(lm2)
[2]))^2/((SSall/SSprod)*sigma)
p.val <-1-pf(F.val, 1, (lm1$df.residual + lm2$df.residual-4))
cat("\n\nTest for equality between two regression slopes\n\n")
c...