This data set contains hourly observations on progesterone and an activity index at hourly intervals since calving on seven dairy cows.
reprocows is a data frame containing 13040 rows.
id is the cow ID,
progesterone is a measurement of the
hormone in ng/L taken from a milk sample, activity is a relative measure
of activity calculated from a pedometer.
There are a large number of missing values as progesterone is measured only at milking time (and at a farm manager's discretion). Missing values in activity occur due to hardware problems can occur with pedometers.
Danish Cattle Research Centre
Peters, A. and Ball, P. (1995), "Reproduction in Cattle," 2nd ed.
1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28
data(reprocows) data(reproai) data(reproppa) tm = 1600 J <- 3 init <- c(1,0,0) trans <- matrix(c(0,0,0,1,0,1,0,1,0),nrow=J) emis <- list(mu=c(0,2.5,0),sigma=c(1,1,1)) N <- as.numeric(table(reprocows$id)) train <- list(x=reprocows$activity,N=N) class(train) <- "hsmm.data" tmp <- gammafit(reproppa * 24) M <- max(N) d <- cbind(dgamma(1:M,shape=tmp$shape,scale=tmp$scale), # ppa sojourn directly estimated from ppa data set dunif(1:M,4,30), # oestrus between 4 and 30 hours dunif(1:M,15*24,40*24)) #not-oestrus between 15 and 40 days startval <- hsmmspec(init,trans,parms.emission=emis,list(d=d,type='gamma'), dens.emission=dnorm.hsmm) #not run (takes some time) #h.activity <- hsmmfit(train,startval,mstep=mstep.norm,maxit=10,M=M,lock.transition=TRUE)
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.