R/pred.R

Defines functions predict.ddpcomp predict.ddp

Documented in predict.ddp predict.ddpcomp

predict.ddp<-function(object,newdata,alpha=0.05,tpred=NULL,...){
if(is.null(tpred)){
tpred<-object$predtime
}
tpred<-tpred/object$high.pct*10

xpred<-model.matrix(object$Terms,newdata,xlev = object$xlevels)
xpred<-xpred[,-1]
xpred<-(xpred-matrix(rep(object$xmean,times=nrow(xpred)),nrow=nrow(xpred), byrow=TRUE))/matrix(rep(2*object$xsd,times=nrow(xpred)),nrow=nrow(xpred), byrow=TRUE)	
predobject<-.Call('DPWeibull_predreg', PACKAGE = 'DPWeibull',
object$alpharec,object$lambdascaled,object$betarec,xpred,tpred,alpha)
class(predobject)<-"predddp"
predobject$tpred<-object$predtime
predobject$alpha<-alpha
predobject$d<-predobject$d/object$high.pct*10
predobject$h<-predobject$h/object$high.pct*10
predobject$dpred<-predobject$dpred/object$high.pct*10
predobject$dpredl<-predobject$dpredl/object$high.pct*10
predobject$dpredu<-predobject$dpredu/object$high.pct*10
predobject$hpred<-predobject$hpred/object$high.pct*10
predobject$hpredl<-predobject$hpredl/object$high.pct*10
predobject$hpredu<-predobject$hpredu/object$high.pct*10
predobject
}

predict.ddpcomp<-function(object,newdata,alpha=0.05,tpred=NULL,...){
if(is.null(tpred)){
tpred<-object$predtime
}
tpred<-tpred/object$high.pct*10

xpred<-model.matrix(object$Terms,newdata,xlev = object$xlevels)
xpred<-xpred[,-1]

xpred<-(xpred-matrix(rep(object$xmean,times=nrow(xpred)),nrow=nrow(xpred), byrow=TRUE))/matrix(rep(2*object$xsd,times=nrow(xpred)),nrow=nrow(xpred), byrow=TRUE)	
predobject<-.Call('DPWeibull_predcompreg', PACKAGE = 'DPWeibull',
object$alpharec1,object$lambdascaled1,object$betarec1,
object$alpharec2,object$lambdascaled2,object$betarec2,object$prec,
xpred,tpred,alpha)
class(predobject)<-"predddpcomp"
predobject$tpred<-object$predtime
predobject$alpha<-alpha
predobject$d<-predobject$d/object$high.pct*10
predobject$h<-predobject$h/object$high.pct*10
predobject$dpred<-predobject$dpred/object$high.pct*10
predobject$dpredl<-predobject$dpredl/object$high.pct*10
predobject$dpredu<-predobject$dpredu/object$high.pct*10
predobject$hpred<-predobject$hpred/object$high.pct*10
predobject$hpredl<-predobject$hpredl/object$high.pct*10
predobject$hpredu<-predobject$hpredu/object$high.pct*10
predobject
}

Try the DPWeibull package in your browser

Any scripts or data that you put into this service are public.

DPWeibull documentation built on Dec. 13, 2021, 1:07 a.m.