Динамические функции в петли с R - PullRequest
Новая
       13

Динамические функции в петли с R

0 голосов
/ 02 ноября 2018

Эта функция вычисляет R ^ 2 в линейную регрессию, соответствующую пакету dp2distance, используя базу данных благосостояния 

 library(p2distance)
 data <- as.matrix(welfare)
 calculoFactoresPonderacion = function(matriz) {
    m <- dim(matriz)[1]
    n <- dim(matriz)[2]
    vec.results <- numeric()
    for (i in 1:(n - 1)) {
        vec.results[i] <- summary(lm(matriz[, i + 1] ~ matriz[, 
            1:i]))$r.squared
     }
     vect1 <- matrix(1, m, 1)
     coefs <- matrix(vec.results, m, n - 1, byrow = TRUE)
     mR.restado <- 1 - coefs
     mFacPond <- cbind(vect1, mR.restado)
    return(mFacPond)
    }

Я хочу вычислить регресию, используя пакет сплайнов через сплайны, но функция " bs " не работает с "matriz [, 1: i]", если i> 1. Точнее в этой части: 

vec.results[i] <- summary(lm(matriz[, i + 1] ~ bs(matriz[, 
            1:i])))$r.squared

Вы можете мне помочь?

...