library(NISTnls)
Chwirut <- Chwirut1[order(Chwirut1$x,Chwirut1$y),]
attach(Chwirut)
loc.f <- function(vP){
beta <- vP[1:3]
-beta[1]*x - log(beta[2] + beta[3]*x)
}
start <- c(b1=0.15, b2=0.005, b3=0.012)
fit <- ssym.nl(log(y),loc.f,start=start,ncs=x,family='Hyperbolic',xi=1,
local.influence=TRUE)
summary(fit)Run the code above in your browser using DataLab