年龄转换率R

Question

data=data.frame("id"=c(1,1,1,1,2,2,2,2,3,3,3,3,4,4,4,4),
                "grade"=c(11,11,12,13,11,12,12,11,13,14,NA,NA,12,11,13,14),
                "age"=c(20,21,22,23,26,27,28,29,19,20,NA,NA,22,23,24,25))

我们有关于学生何时从一个年级升入另一年级的数据，并希望使用纵向数据来估计各年级基于年龄的过渡概率。这对我来说是一个难题。就像马尔科夫一样，但我们希望使用与使用多项式模型可以完成的一样多的基数R。输出理想的是按年龄显示的过渡矩阵。