rs3 {gnlm}R Documentation

Fit a Three-factor Box-Tidwell Nonlinear Response Surface

Description

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

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

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]