[R] NLME and limits on parameter space
    Karsten D Bjerre 
    kdb at kvl.dk
       
    Wed Dec 10 12:16:17 CET 2003
    
    
  
Dear R-users,
How can I impose limits on parameter estimates in NLME as indicated in
the second #'ed line below? 
Can I use the naPattern, some thing like this:  
       nlme(response = yftp(x,ii,cc,b,dd),... , naPattern=
is.na(response) ) ?
yftp<-function(x,ii,cc,b,d){
  tpx<-tpx(ii,cc,b)    # x value for top point of 2nd order polonomial
  tpy<-tpy(ii,cc,b)
  #  if (ii<0 || ii>1 || cc<0|| cc>1|| b<0) {print(NA); return(NA)}
  return( (ii+b*x)*(cc-x)*(x<=tpx) + tpy*(x>tpx)+ d*x)
}
Best whishes,
Karsten
    
    
More information about the R-help
mailing list