Model comparison Bayesian Analysis on visit length

\(Models\ comparison\ trough\ Widely\ applicable\ information\ criterion\ (WAIC),\ using\ the dataset\\ with\ the\ visit\ length\ at\ the\ feeder\ when\ the\ next\ visit\ was\ greater\ than\ or\ equal\ to\ 600sc,\ from\ 7\ trials,\\ in\ pig\ production\)

1. Model comparison visit length as the next visit is greater or equal than 600 seconds (6340 records)

1.1. WAIC criterion

print(m1.600s.waic)

## 
## Computed from 30000 by 6258 log-likelihood matrix
## 
##           Estimate    SE
## elpd_waic -21071.1 211.5
## p_waic       179.3  20.2
## waic       42142.1 423.1

## Warning: 33 (0.5%) p_waic estimates greater than 0.4. We recommend trying loo
## instead.

print(m2.600s.waic)

## 
## Computed from 30000 by 6258 log-likelihood matrix
## 
##           Estimate    SE
## elpd_waic -21069.4 210.6
## p_waic       184.0  20.5
## waic       42138.9 421.2

## Warning: 37 (0.6%) p_waic estimates greater than 0.4. We recommend trying loo
## instead.

print(m3.600s.waic)

## 
## Computed from 30000 by 6258 log-likelihood matrix
## 
##           Estimate    SE
## elpd_waic -21072.1 211.1
## p_waic       186.1  20.8
## waic       42144.2 422.2

## Warning: 37 (0.6%) p_waic estimates greater than 0.4. We recommend trying loo
## instead.

1.2. Model comparison

loo_compare(m1.600s.waic,m2.600s.waic, m3.600s.waic)

##        elpd_diff se_diff
## model2  0.0       0.0   
## model1 -1.6       1.4   
## model3 -2.7       1.4

print(compare(m1.600s.waic,m2.600s.waic, m3.600s.waic))

##              elpd_diff se_diff  elpd_waic p_waic   waic    
## m2.600s.waic      0.0       0.0 -21069.4     184.0  42138.9
## m1.600s.waic     -1.6       1.4 -21071.1     179.3  42142.1
## m3.600s.waic     -2.7       1.4 -21072.1     186.1  42144.2