Joyce Clemente
2024-05-01
#prices at time points for locales
pr1 <- t(byprov53[rn1, 6:8] )
pr2 <- t(byprov53[rn2, 6:8] )
#percent diff between time points
d12 <- round(((pr2 - pr1)/pr1) *100, 2)
#avg and %RSD over time period
avg <- round(apply(byprov53[rn1:rn2, 6:8], 2, mean), 2)
RSD <- round( (apply(byprov53[rn1:rn2, 6:8], 2, sd) / avg) * 100, 2)
#percent diff from national avg per time point
prCAN1 <- round( ( (pr1 - (byprov53[rn1, 6] ))/ (byprov53[rn1, 6]) ) *100, 2)
prCAN2 <- round( ( (pr2 - (byprov53[rn2, 6] ))/ (byprov53[rn2, 6]) ) *100, 2)