2015-01-26 4 views
-2

Мне было интересно, существует ли более полная функция summary() в R, которая, возможно, включает в себя больше модельных показателей, таких как доверительные интервалы вокруг оценок, возможно, лог-правдоподобие, AIC, такие вещи BIC. Я знаю, что довольно легко просто вызвать другие функции на объекте модели, чтобы получить эти показатели, но, черт возьми, было бы круто, если бы все было упаковано в итоговую функцию!Более полная функция summary() в R?

+0

Для какой-то модели? 'summary()' вызывает специфические для объекта методы (например, 'summary.lm()', 'summary.glm()', 'summary.merMod()' и т. д.). Вы также можете написать свою собственную функцию. – gung

+0

в основном lm и glm – moku

+0

Если у вас есть объект 'lm', тогда вызывается' summary() 'на нем –

ответ

1

Не обязательно ДОПОЛНИТЕЛЬНО, но я стараюсь писать свои собственные функции сводки. Затем я могу адаптироваться к конкретному приложению. Кажется, мне всегда нужно «n», и он не содержится в каких-либо доступных резюме.

sumstats<-function(x) { 
    n<-NROW(x) 
    mean<-mean(x) 
    max<-max(x) 
    min<-min (x) 
    SD<- sd(x) 
    hist<-hist(x) 
    summary<-list(n=n,mean=mean,max=max,min=min,SD=SD) 
    return(summary) 
} 
+0

Cool thxs Я сделал это больше вокруг модели, но вот моя функция: summaryRobust <- function (model) { 10 regOutput <- cbind (as.data.frame (summary (model) [4]), as.data.frame (confint (model))) colnames (regOutput) <- c ("Estimate", "Std. Error", "T-Value", "P-Value", "CI 2.5%", "CI 97.5%") метрики <- as.data.frame (cbind (AIC (модель), BIC (модель), logLik (модель))) colnames (метрики) <- c ("AIC", "BIC", "logLik") par (mfrow = c (2,2)) plot (model) list (Model_Output = regOutput, Other_Metrics = показатели) } – moku

0

Это больше того, что я надеялся в сводном функции:

summaryRobust <- function(model){ 
    regOutput <- cbind(as.data.frame(summary(model)[4]), as.data.frame(confint(model))) 
    colnames(regOutput) <- c("Estimate", "Std. Error", "T-Value", "P-Value", "C.I. 2.5%", "C.I. 97.5%") 

    rSq <- summary(model)$r.squared 
    Adj_rSq <- summary(model)$adj.r.squared 
    modelFit <- rbind(rSq, Adj_rSq) 
    rownames(modelFit) <- c("R-squared", "Adjusted R-squared") 
    colnames(modelFit) <- c("") 

    metrics <- as.data.frame(cbind(AIC(model), BIC(model), logLik(model))) 
    colnames(metrics) <- c("AIC", "BIC", "logLik") 
    rownames(metrics) <- c("") 

    par(mfrow=c(2,2)) 
    plot(model) 
    out <- list(Model_Output = regOutput, Model_Fit = modelFit, Other_Metrics = metrics) 
    return(out) 
} 
Смежные вопросы