Я пытаюсь экспортировать вывод моей регрессионной модели в файл .txt или .csv. Мне удалось экспортировать основную информацию из summary (), терминов (например, перехват), p-значений и т. Д. c. Однако экспортированный файл не дает мне стандартной информации об остатках, статистике R2 и F, которая отображается в нижней части вывода.
Вот код, который я использую model = lm(cs~group + sex, data = metadata) #adjusting for sex.
Вывод
Call: lm (formula = cs ~ group + sex , данные = метаданные)
Остатки: Мин. 1 квартал Медиана 3 квартал Макс. -16,917 -7,733 -2,506 6,934 32,494
Коэффициенты: Расчетная стандартная. Ошибка t значение Pr (> | t |) (Перехват) 14,051 3,767 3,731 0,000555 * groupHFD -3,987 4,853 -0,822 0,415880 groupHFD + Exe 1,571 4,853 0,324 0,747817 groupHFD + Exe + Gen 4,874 4,989 0,977 0,334131 groupHFD + Gen -1,189 4,853 -0,245 0,807610 полМужчина 9,633 3,104 3,103 0,003378
--- Значение коды: 0 '*' 0,001 '**' 0,01 '*' 0,05 '.' 0,1 '' 1
Остаточная стандартная ошибка: 10,85 на 43 степенях свободы Множественный R-квадрат: 0,2286, Скорректированный R-квадрат: 0,139 F-статистика c: 2,549 на 5 и 43 DF, p-значение : 0,04168
#Save the regression output as table in cvs table.
write.table(csregression, file = 'csregression.csv', sep = ",", row.names = F, col.names = T)
Файл .csv показывает следующее
Estimate Std. Error t value Pr(>|t|) (Intercept) 14.051 3.767 3.731 0.000555 groupHFD
-3,987 4,853 -0,822 0,415880 groupHFD + Exe 1,571 4,853 0,324 0,747817 groupHFD + Exe + Gen 4,874 4,989 0,977 0,334131 groupHFD + Gen -1,189 4,853 -0,245 0,807610 sexMale 9,633 3,104 3,103 0,003378
Есть ли у кого-нибудь такая же проблема? Как я мог экспортировать недостающую информацию, информацию F-статистики c и R-квадрат?