search for: nonaddit

Is there code available to decompose interactions involving at least one nominal factor with more than 2 levels as described, e.g., by Tukey or by Mandel (1971, Technometrics, 13: 1-18)? Tukey's model: E(y[i,j]) = mu0 + a[i] + b[j] + c*a[i]*b[j], estimating a, b, and c so sum(a) = sum(b)= 0. Mandel essentially describes a singular value decomposition of the interaction. Thanks,

Dear R users: I have some difficulties analizing data with mixed effects NLME and the last version of R. More concretely, I have a repeated measures design with a single group and 2 experimental factors (say A and B) and my interest is to compare additive and nonadditive models. suj rv A B 1 s1 4 a1 b1 2 s1 5 a1 b2 3 s1 7 a1 b3 4 s1 1 a2 b1 5 s1 4 a2 b2 6 s1 2 a2 b3 7 s2 6 a1 b1 8 s2 8 a1 b2 9 s2 10 a1 b3 10 s2 3 a2 b1 11 s2 6 a2 b...

I would like to write a function that computes Tukey's 1 df for nonadditivity. Here is a simplified version of the function I'd like to write: (m is an object created by lm): tukey.test <- function(m) { m1 <- update(m, ~.+I(predict(m)^2)) summary(m1)$coef } The t-test for the added variable is Tukey's test. This won't work: data(BOD) m1 &l...

I'm having problems with environments and update() that I expect have a simple explanation. To illustrate, suppose I wanted to make a very primitive Tukey one-degree-of- freedom for nonadditivity test and naively wrote: nonadd <- function(formula){ f <- lm(formula) v <- f$fitted.values^2 g <- update(f, . ~ . + v) anova(f,g) } x <- rnorm(20) y <- rnorm(20) nonadd(y ~ x) Evidently, update is looking in the environment prod...

Dear listers I am a very newbie with graphs in R. I have a pulmonary function prediction equation in the form of PVC = 1.1 - 0.45*age in years + 0.011*height in cm. How can I draw the corresponding nomogram?. I read the help for the design.nomogram function but it is too difficult for me. Excuse my ignorance. Any direct help will be appreciated. Thanks in advance.

Hi eveybody. I've seen posted before the same question, but no answer. Is there anybody who has implemented the Hosmer-Lemeshow test in R? TIA, Stefano -.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.- r-help mailing list -- Read http://www.ci.tuwien.ac.at/~hornik/R/R-FAQ.html Send "info", "help", or "[un]subscribe" (in the