R语言结果输出到word-stargazer
前言
?????關(guān)于R結(jié)果的導(dǎo)出,有不少的工具可以使用,包括:stargzer、R2word、Sweave R、Knitr、texreg等工具(來源于人大經(jīng)濟(jì)論壇討論貼:如何將R運(yùn)行結(jié)果自動(dòng)保存?)此處僅介紹使用stargazer導(dǎo)出結(jié)果到word.使用的數(shù)據(jù)集為R自帶的attitude數(shù)據(jù)集。
查看數(shù)據(jù)集
head(attitude) #查看前5行 str(attitude) #查看數(shù)據(jù)結(jié)構(gòu)可以看到,該數(shù)據(jù)集包含8個(gè)變量,30條數(shù)據(jù)。
描述性統(tǒng)計(jì)
library('stargazer') #在工作空間展示 stargazer(attitude, title = "Tab了e 1. descriptive statistic",type = "text",no.space=TRUE,summary.stat=c('median','mean','sd','min','max')) #導(dǎo)出到相應(yīng)路徑,結(jié)果為doc的表格 stargazer(attitude, title = "Table 1. descriptive statistic",type ="html", out='C:/Users/xxx/Desktop/descriptive statistic')???注:1)stargazer的summary.stat參數(shù)可以選擇要展示的統(tǒng)計(jì)量,統(tǒng)計(jì)量的簡寫可見stargazer參數(shù)解釋文件P23;
??????????2)summary參數(shù)=FLASE,結(jié)果為查看整個(gè)數(shù)據(jù)集或數(shù)據(jù)集的子集,默認(rèn)為TRUE;
??????????3)rownames設(shè)置為FALSE來去掉行號;
??????????4)其他見文件stargazer參數(shù)解釋文件。
相關(guān)性分析
cor_matr=cor(attitude) #cor_matr=cor(df1,method = 'spearman') #在工作空間展示 stargazer(cor_matr, title = "Table 2. variable correlations",type = "text") #導(dǎo)出到相應(yīng)路徑,結(jié)果為doc的表格 stargazer(cor_matr, title = "Table2. variable correlations",type ="html", out='C:/Users/xxx/Desktop/correlation.doc') #計(jì)算相關(guān)系數(shù)的顯著性 library('psych') corr.test(attitude, use = "complete",method = "pearson", adjust = "none") #可視化相關(guān)系數(shù)矩陣 library('corrplot') corrplot(corr = cor_matr,order = "AOE",type="upper",tl.pos = "d") corrplot(corr = cor_matr,add=TRUE, type="lower", method="number",order="AOE",diag=FALSE,tl.pos="n", cl.pos="n")
?????注:計(jì)算兩個(gè)變量之間的相關(guān)性顯著性可以直接使用cor.test函數(shù)(無需調(diào)包),method默認(rèn)是pearson,alternative默認(rèn)是雙尾(即總體相關(guān)系數(shù)不等于0)。關(guān)于相關(guān)性分析細(xì)節(jié)可見生信筆記文章R語言相關(guān)性分析。
回歸分析
linear.1 <- lm(rating ~ complaints + privileges + learning + raises + critical, data = attitude) linear.2 <- lm(rating ~ complaints + privileges + learning, data = attitude) #結(jié)果導(dǎo)出到工作空間 stargazer(linear.1,linear.2,title="Table 3. Regression Results",type='text',no.space=TRUE) #結(jié)果導(dǎo)出參數(shù)解釋 stargazer(linear.1,linear.2, title="Regression Results", dep.var.labels=c("Overall Rating"), #dv名稱修改 covariate.labels=c("Handling of Complaints","No Special Privileges", "Opportunity to Learn","Performance-Based Raises","Too Critical","Advancement"),#iv名稱修改 omit.stat=c("LL","ser","f"),#忽略的一些統(tǒng)計(jì)量 keep.stat="n"#刪除一些統(tǒng)計(jì)量,保留obersvation,只用一個(gè)參數(shù)即可 order=c("learning", "privileges"), #展示變量的順序,此處未寫全,實(shí)際中需要寫全 ci=TRUE, ci.level=0.90, #報(bào)告90%置信區(qū)間,而非標(biāo)準(zhǔn)誤 single.row=TRUE) #系數(shù)和標(biāo)準(zhǔn)誤在同一行
?????注:有時(shí)候需要隱藏某個(gè)變量的結(jié)果,僅用YES代表加入了模型,如year fixed effect。具體操作是回歸方程中加入factor(year),展示時(shí)候加上omit = c(“year”),omit.labels = c(‘year’)。
總結(jié)
?????1)stargazer輸出到word,需要改變參數(shù)type=‘html’,并加上輸出的路徑;
????? 2)no.space=TRUE解決輸出表格空一行的問題,試試不加這一項(xiàng)就知道差別了;
????? 3)描述性統(tǒng)計(jì)輸出想要的統(tǒng)計(jì)量需要統(tǒng)計(jì)量的簡寫,回歸結(jié)果導(dǎo)出刪除一些不想要的統(tǒng)計(jì)量,都可通過查看文件stargazer參數(shù)解釋解決;
????? 4)有時(shí)候需要隱藏某個(gè)變量的結(jié)果,僅用YES代表加入了模型,如year fixed effect。具體操作是回歸方程中加入factor(year),展示時(shí)候加上omit = c(“year”),omit.labels = c(‘year’)
????? 5)對導(dǎo)出結(jié)果不滿意的,可以在word中手動(dòng)修改。
?????另:附上官方示例文件-官方示例文件
?????最后:總結(jié)不易,轉(zhuǎn)載時(shí)請寫明出處,也可以打賞一丟丟哦!
總結(jié)
以上是生活随笔為你收集整理的R语言结果输出到word-stargazer的全部內(nèi)容,希望文章能夠幫你解決所遇到的問題。
- 上一篇: python 变量大小,进程和内存信息
- 下一篇: jenkins配置自动部署java程序