tests/insem.R

library(parfm)
data(insem)
head(insem)

insem$TimeMonths <- insem$Time * 12 / 365.25
  
set.seed(1)
insem <- insem[sample(1:nrow(insem), 150),]
parfm(Surv(TimeMonths, Status)~Heifer, cluster="Herd", data=insem,
             dist="weibull", frailty="gamma")

Try the parfm package in your browser

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

parfm documentation built on May 2, 2019, 5 p.m.