我得到了一个数据集，并进行如下编码以估计方差分量。在数据集中，“ y”是响应，在线性混合模型中，“燃料”已固定，而“天”具有随机影响。由于“天”，我可以估计随机效应，并且没有选择具有对应于随机效应估计的p值。如何估算p值？还有其他包装吗？任何答案将不胜感激。
day<-c(1, 1, 2, 2, 3, 3)
fuel<-c(1,1,1,2,2,2)
y<-c(32.32048, 32.28583, 31.97912, 33.61024, 33.33915, 33.63342)
dataTB<-data.frame(day,fuel,y)
library(nlme)
contrastTB.lme<-lme(y~fuel, random= ~1 | day,data=dataTB)
summary(contrastTB.lme)