Предыдущий пользователь спросил Как мне добавить доверительные интервалы к отношениям шансов в таблице звездочета? и наметили четкое решение проблемы.
В настоящее время я печатаю свои таблицы вручную, что отнимает очень много времени. пример моей напечатанной таблицы. Вот ссылка на используемый файл .txt.
Моя модель имеет размер как зависимую переменную (категориальную) и пол (категориальный), возраст (непрерывный) и год (непрерывный) как независимые переменные. Я использую mlogit
для моделирования взаимосвязи между переменными.
Код, который я использовал для модели, выглядит следующим образом:
tattoo <- read.table("https://ndownloader.figshare.com/files/6920972",
header=TRUE, na.strings=c("unk", "NA"))
library(mlogit)
Tat<-mlogit.data(tattoo, varying=NULL, shape="wide", choice="size", id.var="date")
ml.Tat<-mlogit(size~1|age+sex+yy, Tat, reflevel="small", id.var="date")
library(stargazer)
OR.vector<-exp(ml.Tat$coef)
CI.vector<-exp(confint(ml.Tat))
p.values<-summary(ml.Tat)$CoefTable[,4]
#table with odds ratios and confidence intervals
stargazer(ml.Tat, coef=list(OR.vector), ci=TRUE, ci.custom=list(CI.vector), single.row=T, type="text", star.cutoffs=c(0.05,0.01,0.001), out="table1.txt", digits=4)
#table with coefficients and standard errors
stargazer(ml.Tat, type="text", single.row=TRUE, star.cutoffs=c(0.05,0.01,0.001), out="table1.txt", digits=4)
Код stargazer
, который я пробовал, показан ниже для небольшой части моих данных:
library(stargazer)
OR.vector<-exp(ml.Tat$coef)
CI.vector<-exp(confint(ml.Tat))
p.values<-summary(ml.Tat)$CoefTable[,4] #incorrect # of dimensions, unsure how to determine dimensions
stargazer(ml.Tat, coef=list(OR.vector), ci=TRUE, ci.custom=list(CI.vector), single.row=T, type="text", star.cutoffs=c(0.05,0.01,0.001), out="table1.txt", digits=4) #gives odds ratio (2.5%CI, 97.5%CI)
Выходные данные отношения шансов и доверительного интервала:
stargazer(ml.Tat, type="text", single.row=TRUE, star.cutoffs=c(0.05,0.01,0.001), out="table1.txt", digits=4) #gives coeff (SE)`
Я могу комбинировать отношения шансов с доверительными интервалами или стандартными ошибками или коэффициенты с доверительными интервалами и стандартными ошибками, но когда я пишу все три вместе, функция ci=TRUE
, кажется, перезаписывает значение SE по умолчанию.
Для моей диссертации мне нужны таблицы, чтобы показать коэффициенты, стандартные ошибки, доверительные интервалы и отношения шансов (и p-значения в некотором формате). Есть ли способ у звездочета включить все четыре вещи? Возможно, в двух разных столбцах? Я могу экспортировать таблицы в Excel, однако без всех четырех вещей в одной и той же таблице звездочета я застрял вручную, соединяя две приведенные выше таблицы вместе. Это не имеет большого значения для 1 таблицы, но я работаю с 36 моделями, всем которым нужны таблицы (для моей диссертации).
Как с помощью звездочета показать все четыре вещи? (отношение шансов, доверительные интервалы, коэффициенты и стандартные ошибки)
glm
. Я не думаю, что он может создать такую модель. У меня есть загруженные данные, и фрейм данных показывает 18 171 строку с переменной размера, имеющей 4 уровня, поэтому я предполагаю, что вы исключили 9 случаев unk. Если вы хотите получить конкретный ответ на эту проблему, вам необходимо включить код обработки и моделирования данных. - person IRTFM   schedule 20.11.2016stargazer
принимает несколько моделей и включает их в отдельные столбцы. Поэтому создайте вторую модель (например, ml.TatOR) и замените коэффициенты отношениями шансов. Затем в Stargazer включите обе модели и соответствующие аргументы:stargazer(ml.Tat, ml.TatOR, type="text",ci = c(F,T),column.labels=c("coefficients","odds ratio"), single.row=TRUE)
- person paqmo   schedule 20.11.2016ml.Tat
, чтобы вstargazer
выводе были показаны четыре вещи, которые мне нужны в моей таблице? Моя текущая модель выглядит какml.Tat<-mlogit(size~1|age+sex+yy, Tat, reflevel="small", id.var="date")
, и я ничего не уточняю о выходе модели. Как перейти от вывода, который показывает коэффициенты и стандартные ошибки, к отношениям шансов и доверительным интервалам? - person Blundering Ecologist   schedule 20.11.2016ml.TatOR <-mlogit(size~1|age+sex+yy, Tat, reflevel="small", id.var="date")
. Заменить коэффициенты отношениями шансов:ml.TatOR$coef <- exp(ml.TatOR$coef)
. Затем поместите обе модели в свойstargazer
вызов, как указано выше. Аргументci = c(F,T)
подавляет доверительный интервал в первом столбце, но не во втором. - person paqmo   schedule 20.11.2016