Fit a Three-factor Box-Tidwell Nonlinear Response Surface

Usage

rs3(y, x1, x2, x3, power=c(1,1,1), weight=rep(1,length(y)),
	family=normal, iterlim=20)

Arguments

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

Description

rs3 fits a three-covariate power-transformed response surface by iterating the function, glm.

Value

A list of class, rs, is returned containing the model and the power estimates.

Author(s)

J.K. Lindsey

See Also

lm, glm, gnlr, gnlr3, fmr, rs2.

Examples

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)
rs3(y, x1, x2, x3, family=poisson)


[Package Contents]