R Under development (unstable) (2025-07-22 r88445 ucrt) -- "Unsuffered Consequences"
Copyright (C) 2025 The R Foundation for Statistical Computing
Platform: x86_64-w64-mingw32/x64

R is free software and comes with ABSOLUTELY NO WARRANTY.
You are welcome to redistribute it under certain conditions.
Type 'license()' or 'licence()' for distribution details.

R is a collaborative project with many contributors.
Type 'contributors()' for more information and
'citation()' on how to cite R or R packages in publications.

Type 'demo()' for some demos, 'help()' for on-line help, or
'help.start()' for an HTML browser interface to help.
Type 'q()' to quit R.

> library(parfm)
Loading required package: survival
Loading required package: optimx
> data(culling)
> head(culling)
  Cowid Time Status Herd Timeassess LogSCC
1     1  331      0    1         13   4.09
2     2  312      0    1          6   4.83
3     3   96      1    1         11   3.93
4     4  314      0    1         13   3.76
5     5  181      1    1          7   6.29
6     6  331      0    1         11   3.33
> culling <- culling[culling$Time > 0,]
> culling$TimeMonths <- culling$Time * 12 / 365.25
> 
> set.seed(1)
> culling <- culling[sample(1:nrow(culling), 150),]
> parfm(Surv(TimeMonths, Status)~LogSCC, cluster="Herd", data=culling,
+                dist="exponential", frailty="gamma")

Frailty distribution: gamma 
Baseline hazard distribution: Exponential 
Loglikelihood: -180.345 

       ESTIMATE SE    p-val 
theta  0.000    0.538       
lambda 0.022    0.013       
LogSCC 0.043    0.124 0.728 
---
Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1

Kendall's Tau: 0 
Warning message:
In parfm(Surv(TimeMonths, Status) ~ LogSCC, cluster = "Herd", data = culling,  :
  optimisation procedure did not converge,
              conv = : see ?optimx for details
> 
> proc.time()
   user  system elapsed 
   5.85    0.43    6.23