rs3: Three-factor Box-Tidwell Nonlinear Response Surface Models

View source: R/rs.r

rs3R Documentation

Three-factor Box-Tidwell Nonlinear Response Surface Models

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(x1)),
  family = gaussian,
  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, 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)
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)


gnlm documentation built on June 8, 2025, 11:32 a.m.

Related to rs3 in gnlm...