R语言里怎么做十折交叉验证？

SATjustice

在R里怎么才能自动做十折交叉验证啊？

难道只能自己做十遍然后计算十次误差再取平均么？

nan.xiao

很多方法都内置了工具

boot::cv.glm<br />
lars::cv.lars<br />
elasticnet::cv.enet<br />
GAMens::GAMens.cv

</p>

可以看到大概也就是那样不知有没有通用性很强的工具自己的方法自定义一个也好

可以统一到一些框架下来做比如caret 提供了多种做验证的方法在从很多方法中选择的时候就很方便了就是指定一个参数的事情

fisheat

ten fold crossvalidation 好像不是做十遍的意思啊

而是应该是随机分成十组，9/10 development，1/10 validation。

nan.xiao

LS正解

楼主可能想说repeated K-fold CV吧。。。

SATjustice

回复第4楼的 nan.xiao：我就是4楼的意思，我初学R语言，想问问该怎么实现啊？比如我用lm来作模型，还需要手动设置训练集和测试集么？还是说电脑自动帮我搞定了？谢谢啦

nan.xiao

现存的方法很多都会自带现成的工具

线性模型boot::cv.glm应该可以直接做

用caret的话

require(caret)</p>
<p>fitControl = trainControl(<br />
    method = "repeatedcv",<br />
    number = 10,<br />
	repeats = 5,<br />
	returnResamp = "all")</p>
<p>x = jitter(10:30)<br />
y = jitter(40:60)</p>
<p>glmFit = train(y, x,<br />
    method = "glm",<br />
    trControl = fitControl)</p>
<p>print(glmFit)

</p>

修改train()中的参数method 就可以比较几十种不同的分类/回归方法

修改fitControl()中的参数method即可选择不同的验证方法

目前有boot, boot632, cv, repeatedcv, LOOCV, LGOCV, oob这几种

修改number指定折数或者重抽样迭代次数

修改repeats指定repeatedcv做多少次重复

不过以上都是浮云你依然完全拥有自己去实现的自由。。。

flyhi

顶LS

SATjustice

回复第6楼的 nan.xiao：谢谢啦，我去试试！

juvenchan

回复第1楼的 SATjustice：我想問你有code嗎?

我要做十次十折, 但我真的不會寫

vz33

回复第6楼的 nan.xiao：应该是

glmFit = train(x, y,

method = "glm",

trControl = fitControl)

vz33

回复第9楼的 juvenchan：

<br />
require(caret)</p>
<p>fitControl = trainControl(<br />
    method = "repeatedcv",<br />
    number = 10,<br />
    repeats = 10,<br />
    returnResamp = "all")</p>
<p>x = jitter(10:30) # 自己修改这两行，用你的data<br />
y = jitter(40:60)#</p>
<p>glmFit = train(x, y,<br />
    method = "glm",   # glm 可以换很多其他方法<br />
    trControl = fitControl)</p>
<p>summary (glmFit)<br />

</p>

alma2004

又学习了[s:13]

xiaoqianxixi

回复第11楼的 vz33：您好，我想问一下，怎么把我的数据跟这个交叉验证联系起来呢？语句要怎么写呢？谢谢！

xiaoqianxixi

回复第6楼的肖楠：您好，如果我用randomForest，语句如下：

library（randomForest）

x=read.table("1.txt")

set.seed(150)

x.rf<-randomForest(V22~.,data=x,importance=TRUE,proximity=TRUE)

print(x.rf)

我想加入5-折交叉验证，那么我需要在哪加语句呢？

x = jitter(10:30)

y = jitter(40:60)

以上这两行是什么意思呢？

cr7madrid

回复第14楼的 xiaoqianxixi：在trainControl里修改method=“cv”和number=5

HarryYu

统计真是深奥阿，我还是一心一意的搞好金融和算法，统计上向大家学习

neversarah

学习了

nan.xiao

发现最近经常手写并行可重现的重复n次k折交叉验证以保持手感。[s:11]

这个东西非常容易写错，经典程度可以作为 statistical programmer 的面试题。

platinum_

require(caret)
fitControl = trainControl(
method = "oob",
number = 10,
repeats = 5,
returnResamp = "all")

glmFit = train(Train80$FICO,Train80$FlagD60, #Train80$FlagD60 这个是个分类变量，
method = "rf", #但是我换成连续变量还是报同样的错orz
trControl = fitControl)
print(glmFit)

报错：Error in if (p <= len) { : argument is of length zero
不知道是什么原因，求解答，谢谢！

cjsun0000

[未知用户]
楼主，你交叉验证的编程你会了吗