Как использовать сводную функцию внутри функции ldply()-summarise для извлечения p-значений?
Пример данных:
(предустановлен датафрейм "Пуромицин")
library(reshape2)
library(plyr)
Puromycin.m <- melt( Puromycin , id=c("state") )
Puro.models <- dlply( Puromycin.m , .(variable) , glm , formula = state ~ value ,
family = binomial )
Я могу создать этот фрейм данных с извлеченными результатами:
ldply( Puro.models , summarise , "n in each model" = length(fitted.values) ,
"Coefficients" = coefficients[2] )
Но я не могу извлечь p-значения таким же образом. Я думал, что это сработает, но это не так:
ldply( Puro.models , summarise ,
"n in each model" = length(fitted.values) ,
"Coefficients" = coefficients[2],
"P-value" = function(x) summary(x)$coef[2,4] )
Как я могу извлечь p-значения в этот фрейм данных :) Пожалуйста, помогите!