rs3 {gnlm} | R Documentation |
rs3
fits a three-covariate power-transformed response surface by
iterating the function, glm
.
rs3(y, x1, x2, x3, power=c(1,1,1), weight=rep(1,length(y)), family=normal, iterlim=20)
y |
Response variable |
x1 |
First covariate |
x2 |
Second covariate |
x3 |
Third covariate |
power |
Initial estimates of the three power transformations |
weight |
Weight vector |
family |
glm family |
iterlim |
Iteration limit |
A list of class, rs
, is returned containing the model and the power
estimates.
J.K. Lindsey
x1 <- rep(1:4,5) x2 <- rep(1:5,rep(4,5)) x3 <- c(rep(1:3,6),1,2) #y <- rpois(20,1+2*sqrt(x1)+3*log(x2)+1/x3+4*x1+log(x2)^2+1/x3^2+ # 2*sqrt(x1)*log(x2)+sqrt(x1)/x3+log(x2)/x3) y <- c(9,11,14,33,11,19,20,27,22,32,24,24,20,28,25,41,26,31,37,34) rs3(y, x1, x2, x3, family=poisson)