24小时热门版块排行榜    

北京石油化工学院2026年研究生招生接收调剂公告
查看: 650  |  回复: 7
当前主题已经存档。

loappleve

金虫 (小有名气)

[交流] 求助 如何实现leava one out

用多元线性回归建模以后 如果做leave one out 检验
有什么软件可以实现 自己写的话 需要输出什么结果

[ Last edited by csfn on 2008-12-29 at 20:16 ]
回复此楼
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

snoopyzhao

至尊木虫 (职业作家)


spur(金币+1,VIP+0): 感谢参与!欢迎常来!
我以前在这里贴过一个用 R 语言写的 leave-one-out 的小程序,可以找找看,呵呵
2楼2008-01-11 14:19:04
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

2006Jerry

金虫 (正式写手)

没有搜到啊!
能否请snoopyzhao发一份啊!
谢谢!
gzg_123@tom.com
3楼2008-01-11 14:43:02
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

snoopyzhao

至尊木虫 (职业作家)

★ ★ ★ ★ ★
zzgyb(金币+5,VIP+0):谢谢您的参与,欢迎您常来计算模拟版。
两年前写的,仅供参考。
CODE:
cv.lm <- function(obj, loo = TRUE) {
#
# The Leave-One-Out (LOO) and/or Leave-Group-Out (LGO) Cross-Validation in R for (Multiple) Linear Regression.
#
# Input:
#   obj: the model of MLR
#   loo: logic, if FALSE, do LGO CV
#
# Note:
#   the following two varible should be assigned according to the number of observation in the data set.
#   m: the number of loop when doing validation.
#   n: the number of removed observations.
#
# Output:
#   q.squared: cross-validation relation coefficient.
#   SDEP: Standard Deviation of Error of Predictions
#   newsq: variance in Y explained only for LOO CV
#
# Usage:
# loo <- cv.lm(obj)
# lgo <- cv.lm(obj, loo = FALSE)
#
# Copyright (C) 2005
#
        data <- data.frame(y = obj$model[,1], x = obj$model[, -1])
        N <- nrow(data)
        if (loo == TRUE) {
                ytest <- data[, 1]
                newrsq <- numeric(N)
                ypred <- numeric(N)
                for(i in 1:N){
                        newtrain <- data[-i, ]
                        xtest <- data[i, -1]
                        newfm <- lm(y ~., data = newtrain)
                        newrsq[i] <- summary(newfm)$r.squared
                        ypred[i] <- predict(newfm, xtest)
                }
        }
        else {
                m <- 40
                n <- round(N / 10)
                newrsq <- numeric(m)
                ytest <- numeric(m*n)
                ypred <- numeric(m*n)
                for (i in 1:m) {
                        j <- sample(N, n, replace = FALSE)
                        newtrain <- data[-j, ]
                        xtest <- data[j, -1]
                        ytest[((i-1) * n + 1):(i * n)] <- data[j, 1]
                        newfm <- lm(y ~., data = newtrain)
                        ypred[((i-1) * n + 1):(i * n)] <- predict(newfm, xtest)
                }
        }
        q.squared <- 1 - sum((ytest - ypred)^2) / sum((ytest - mean(ytest))^2)
        if (loo == TRUE) {
                SDEP <- sqrt(sum((ytest - ypred)^2) / N)
        }
        else {
                SDEP <- sqrt(sum((ytest - ypred)^2) / (m * n))
        }
        if (loo == TRUE) {
                return(list(q.squared = q.squared, SDEP = SDEP, newrsq = newrsq))
        }
        else {
                return(list(q.squared = q.squared, SDEP = SDEP))
        }
}

4楼2008-01-12 10:09:19
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

loappleve

金虫 (小有名气)

clc
clear
[filename, pathname]= uigetfile('*.xls');
file=[pathname filename];
a=xlsread(file);
[m,n]=size(a);
b=zeros(m-1,n);
c=zeros(m,n-1);
d=zeros(m,1);
k=1;
for i=1:m
    k=1;
    for j=1:m
        if j~=i
            b(k,=a(j,;
            k=k+1;
        end
    end
    c(i,=(b(:,1:n-1)\b(:,n))';
    d(i)=a(i,1:n-1)*(c(i,');
end
5楼2008-01-12 12:12:21
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

2006Jerry

金虫 (正式写手)

楼上的:你发的这个程序是干什么的呢?
6楼2008-01-14 18:20:35
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

juju550

金虫 (小有名气)

请高手们发一个可用的留一法Matlab程序啥,自己编的好象有些问题

着急着想用啊,先谢过了!

顺便把调用格式也说明一下

太感谢了!
7楼2008-01-16 23:59:55
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

匿名

用户注销 (小有名气)


本帖仅楼主可见
8楼2008-01-31 11:15:20
已阅   申请仿真EPI   回复此楼   编辑   查看我的主页
相关版块跳转 我要订阅楼主 loappleve 的主题更新
普通表情 高级回复 (可上传附件)
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 求调剂 +3 111623 2026-04-04 3/150 2026-04-05 01:35 by wxiongid
[考研] 材料调剂 +11 一样YWY 2026-04-02 13/650 2026-04-04 23:10 by 无际的草原
[考研] 化学357分,考研调剂 +8 .Starry. 2026-04-04 9/450 2026-04-04 22:52 by 理想在实践之中
[考研] 求生物学学硕调剂——364分 +7 云朵遛弯指南 2026-04-04 7/350 2026-04-04 22:49 by zhyzzh
[考研] 一志愿南航,数一英一学硕317求调剂!! +3 Acaciad 2026-04-04 3/150 2026-04-04 22:29 by lbsjt
[考研] 085600调剂 +4 1amJJ 2026-04-02 4/200 2026-04-04 21:53 by hemengdong
[考研] 341求调剂 +3 洛多罗 2026-04-02 4/200 2026-04-04 21:36 by 智能智慧
[考研] 化学调剂 +12 艾志恒 2026-04-03 13/650 2026-04-04 20:11 by xhai2011
[考研] 081700化学工程与技术 一志愿中海洋 323 求调剂学校 +15 披星河 2026-04-03 15/750 2026-04-04 16:25 by dongzh2009
[考研] 311求调剂 +20 zchqwer 2026-04-01 22/1100 2026-04-03 22:09 by lglzsd
[考研] 311求调剂 +11 勇敢的小吴 2026-04-02 11/550 2026-04-03 21:46 by qlm5820
[考研] 293求调剂 +5 末未mm 2026-04-02 6/300 2026-04-03 15:20 by 王保杰33
[考研] 英一数一408,总分284,二战真诚求调剂 +13 12.27 2026-03-30 15/750 2026-04-03 14:41 by 氮气气气
[考研] 工科341分调剂 +3 洛多罗 2026-04-03 3/150 2026-04-03 14:20 by 1753564080
[考研] 085600专硕材料与化工348分求调剂 +10 上学啦! 2026-04-01 11/550 2026-04-03 14:13 by 百灵童888
[考研] 一志愿陕西师范大学生物学317分 +5 1563日。 2026-04-02 5/250 2026-04-03 06:58 by ilovexiaobin
[考研] 279求调剂 +5 傅文秋 2026-04-02 5/250 2026-04-02 18:10 by 笔落锦州
[考研] 08生物与医药专硕初试346找调剂 +6 dianeeee 2026-04-01 7/350 2026-04-02 08:23 by guoweigw
[考研] 求调剂0703 +5 周嘉尧 2026-03-31 8/400 2026-04-01 20:32 by ltltkkk
[考研] 085600,320分求调剂 +5 大馋小子 2026-04-01 6/300 2026-04-01 19:40 by 唐沐儿
信息提示
请填处理意见