search for: timetyp

...e left part of the plot are the time values before this date (scale_x_datetime major = 1 year), and the right part of the plot are the time values after this date (scale_x_datetime major=1 day). Hereby also the R-code (simplified): ggplot(data_plot.melt,aes(timevalue,ID)) + geom_point(aes(groups=timetype,colour=timetype,shape=timetype)) + facet_grid(type ~pos,scales="free",space="free") + xlab(NULL) + ylab(NULL) The scales of y has to be free, because the number of ID's per type differ. The scales of x has to be free, so the scales differ in the left and right part of th...

...~ dnorm( 0.0, 1.0E-4)# definition of sigmas2<-1/taus <-sqrt(s2)# calculation of the sample variancefor (i in 1:n){ c.time[i]<-time[i]-mean(time[]) }sy2 <- inprod( c.time[], c.time[] )/(n-1)# calculation of Bayesian version R squaredR2B <- 1 - s2/sy2# Expected y for a typical delivery timetypical.y <- beta0 + beta1 * mean(cases[]) + beta2 * mean(distance[])}INITSlist( tau=1, beta0=1, beta1=0, beta2=0 )DATA (LIST)list( n=25,time = c(16.68, 11.5, 12.03, 14.88, 13.75, 18.11, 8, 17.83,79.24, 21.5, 40.33, 21, 13.5, 19.75, 24, 29, 15.35,19, 9.5, 35.1, 17.9, 52.32, 18.75, 19.83, 10.75),dist...