24小时热门版块排行榜    

北京石油化工学院2026年研究生招生接收调剂公告
查看: 644  |  回复: 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 的主题更新
普通表情 高级回复 (可上传附件)
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 320分人工智能调剂 +7 振—TZ 2026-04-03 7/350 2026-04-05 00:42 by chongya
[考研] 0817化学工程与技术求调剂,一志愿中海洋319 +8 lv945 2026-04-04 8/400 2026-04-04 23:41 by lqwchd
[考研] 290求调剂 +7 luoziheng 2026-04-04 7/350 2026-04-04 23:17 by lqwchd
[考研] 322求调剂:一志愿湖南大学 材料与化工(085600),已过六级。 +15 XX小邓 2026-03-29 17/850 2026-04-04 22:28 by hemengdong
[考研] 323求调剂 +8 李佳乐1 2026-04-04 8/400 2026-04-04 22:26 by hemengdong
[考研] 材料调剂 +12 一样YWY 2026-04-03 12/600 2026-04-04 21:46 by hemengdong
[考研] 278求调剂 +3 依旧! 2026-04-02 4/200 2026-04-04 20:27 by 蓝云思雨
[考研] 一志愿华南师范361分,化学求调剂 +7 Nicole88888 2026-04-01 7/350 2026-04-04 18:28 by macy2011
[考研] 22408求调剂 354分 可跨专业 +3 hannnnnnn 2026-04-04 3/150 2026-04-04 14:35 by 土木硕士招生
[考研] 311求调剂 +20 zchqwer 2026-04-01 22/1100 2026-04-03 22:09 by lglzsd
[考研] 求调剂机会 +5 意染ivy 2026-04-03 5/250 2026-04-03 15:13 by qoooooo614
[考研] 325分化学调剂 +5 15771691647 2026-04-02 5/250 2026-04-03 09:58 by ChemPharm
[考研] 262求调剂 +6 励志一定发文章 2026-04-02 7/350 2026-04-03 09:54 by linyelide
[考研] 调剂 +3 osbbx 2026-04-02 3/150 2026-04-03 07:47 by cc8418
[考研] 330分求调剂 +15 qzenlc 2026-03-29 15/750 2026-04-02 21:41 by 百灵童888
[考研] 求调剂,一志愿 南京航空航天大学 ,080500材料科学与工程学硕,总分289分 +11 @taotao 2026-03-29 11/550 2026-04-02 10:04 by realme321
[考研] 一志愿北京科技大学085601材料工程英一数二初试总分335求调剂 +5 双马尾痞老板2 2026-03-31 5/250 2026-04-01 09:04 by oooqiao
[考研] 一志愿西交大080500材料学硕349 +6 jqx1258 2026-03-31 7/350 2026-03-31 21:08 by yuq
[基金申请] 面上5B能上会吗? +8 redcom 2026-03-29 8/400 2026-03-31 15:53 by niuailing
[考研] 一志愿中海洋320化学工程与技术学硕求调剂 +8 披星河 2026-03-30 8/400 2026-03-31 08:53 by lbsjt
信息提示
请填处理意见