±±¾©Ê¯ÓÍ»¯¹¤Ñ§Ôº2026ÄêÑо¿ÉúÕÐÉú½ÓÊÕµ÷¼Á¹«¸æ
²é¿´: 647  |  »Ø¸´: 7
µ±Ç°Ö÷ÌâÒѾ­´æµµ¡£
µ±Ç°Ö»ÏÔʾÂú×ãÖ¸¶¨Ìõ¼þµÄ»ØÌû£¬µã»÷ÕâÀï²é¿´±¾»°ÌâµÄËùÓлØÌû

loappleve

½ð³æ (СÓÐÃûÆø)

[½»Á÷] ÇóÖú ÈçºÎʵÏÖleava one out

ÓöàÔªÏßÐԻع齨ģÒÔºó Èç¹û×öleave one out ¼ìÑé
ÓÐʲôÈí¼þ¿ÉÒÔʵÏÖ ×Ô¼ºÐ´µÄ»° ÐèÒªÊä³öʲô½á¹û

[ Last edited by csfn on 2008-12-29 at 20:16 ]
»Ø¸´´ËÂ¥

» ²ÂÄãϲ»¶

ÒÑÔÄ   »Ø¸´´ËÂ¥   ¹Ø×¢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µÄ»ØÌû
²é¿´È«²¿ 8 ¸ö»Ø´ð

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µÄ»ØÌû
×î¾ßÈËÆøÈÈÌûÍÆ¼ö [²é¿´È«²¿] ×÷Õß »Ø/¿´ ×îºó·¢±í
[¿¼ÑÐ] ²ÄÁϹ¤³Ì085601Êý¶þÓ¢Ò»335Çóµ÷¼Á +6 Ë«ÂíβƦÀϰå2 2026-03-31 6/300 2026-04-04 22:29 by hemengdong
[¿¼ÑÐ] µ÷¼Á +4 ÊÇ¿ÉÀÖ²»ÊÇ¿ÉÀÖ 2026-04-04 4/200 2026-04-04 19:41 by ÌÆãå¶ù
[¿¼ÑÐ] 085600£¬320·ÖÇóµ÷¼Á +14 ´ó²öС×Ó 2026-04-04 15/750 2026-04-04 16:27 by Î޼ʵIJÝÔ­
[¿¼ÑÐ] 334Çóµ÷¼Á +8 ÔøÑöÖ® 2026-04-03 8/400 2026-04-04 11:16 by w_xuqing
[¿¼ÑÐ] »¯Ñ§µ÷¼ÁÇóÖú +6 LULONG1 2026-04-03 6/300 2026-04-03 23:13 by qzxyhcsy
[¿¼ÑÐ] ¿¼Ñе÷¼Á +3 15615482637 2026-04-03 3/150 2026-04-03 22:50 by ms629
[¿¼ÑÐ] 350Ò»Ö¾Ô¸±±¾©º½¿Õº½Ìì´óѧ08500²ÄÁÏ¿ÆÑ§Ó빤³ÌÇóµ÷¼Á +5 kjnasfss 2026-04-03 5/250 2026-04-03 22:29 by Î޼ʵIJÝÔ­
[¿¼ÑÐ] 343Çóµ÷¼Á +8 °®î¿°í 2026-03-29 8/400 2026-04-03 22:19 by –¹Æ?
[¿¼ÑÐ] 281Çóµ÷¼Á +10 aaawhy 2026-04-03 10/500 2026-04-03 21:42 by lbsjt
[¿¼ÑÐ] 266Çóµ÷¼Á +18 ÑôÑôÍÛÈû 2026-04-01 18/900 2026-04-03 18:38 by zllcz
[¿¼ÑÐ] 0856£¬269·ÖÇóµ÷¼Á +15 ÓÐѧÉϾÍÐÐÇóÇóÁ 2026-03-30 18/900 2026-04-03 16:50 by melodiousnow
[¿¼ÑÐ] ¹¤¿Æ 267Çóµ÷¼Á +5 wanwan00 2026-04-02 7/350 2026-04-03 14:14 by zhangdingwa
[¿¼ÑÐ] 338Çóµ÷¼Á£¬Ò»Ö¾Ô¸ÄÜÔ´¶¯Á¦£¬ÍâÓïÊÇÈÕÓï203 +5 zzz£¬£¬r 2026-04-02 5/250 2026-04-03 09:45 by À¶ÔÆË¼Óê
[¿¼ÑÐ] Çóµ÷¼Á +7 Aniyaio 2026-04-02 7/350 2026-04-02 16:42 by zzsw+
[¿¼ÑÐ] 346Çóµ÷¼Á +5 Ö£³ÏÀÖ 2026-04-02 5/250 2026-04-02 16:38 by SZW_UJN
[¿¼ÑÐ] һ־Ըͬ¼Ã´óѧ323·Ö£¨080500£©Çóµ÷¼Á +6 yikeniu 2026-04-01 6/300 2026-04-02 14:19 by smileboy2006
[¿¼ÑÐ] 08¹¤¿ÆÇóµ÷¼Á290·Ö +5 1314Åõ»¨ 2026-04-02 8/400 2026-04-02 13:16 by ÇÇßÕßÕßÕ
[¿¼ÑÐ] 085600£¬320·ÖÇóµ÷¼Á +5 ´ó²öС×Ó 2026-04-01 6/300 2026-04-01 19:40 by ÌÆãå¶ù
[¿¼ÑÐ] 313Çóµ÷¼Á +6 Âô¸ö¹Ø×Ó°É 2026-03-31 6/300 2026-03-31 10:58 by Jaylen.
[¿¼ÑÐ] Ò»Ö¾Ô¸»ªÖÐʦ·¶»¯Ñ§332·ÖÇóµ÷¼Á +3 Lyy930824@ 2026-03-29 3/150 2026-03-30 20:15 by DHUSHUAI
ÐÅÏ¢Ìáʾ
ÇëÌî´¦ÀíÒâ¼û