我也看过这个网页，但是具体的计算，怎么计算？

Here is an example:
# First estimate a model containing the intercept only

fit1 <- glm(Species ~.,family=binomial(link='logit'),data=ir_train)
summary(fit1)

#Next estimate a model containing the predictors, X1 and X2.

fit2 <- glm(Species ~X1 + X2,family=binomial(link='logit'),data=ir_train)
summary(fit2)

#Find LM0 from summary(fit1) and LMF from summary(fit2)
# Use the equation 1-(LM0/LMF)^2/N

谢谢，这个程序很清楚。
我是在十折交叉验证时出现的问题。
（1）100个数数据随机分成十份（train（90个数据）和validation（10个数据）；
（2）train（90个数据）做回归分析，newlm <- glm(model5,data=train,family=binomial())
（3）用validation（10个数据）做样本，计算估计值，newpred <- predict(newlm,newdata=validation,type ="response")
（4）然后在选择90个数据做回归分析，另外10个数据估计。
（5）重复10次，所有的数据都有预测值了。
（6）计算Cox-Snell R2 （不会计算了）。
代码：
见附件！

这个具体公式怎么算呀？急求呀、、、