0
我是R新手 我想將描述性統計量和直方圖都添加到pdf中。R:將文本和圖形輸出到pdf
以下代碼使用ggplot2成功生成兩個直方圖。但是描述(來自心理包)函數不會出現在pdf中。 我該如何包含兩者?
library(psych)
library(foreign)
library(nnet)
library(ggplot2)
library(reshape2)
# direct output to a file
sink("C:\\Users\\jake\\Dropbox\\__iKoda\\datafiles\\OutputR.txt", append=FALSE, split=TRUE)
gc()
memory.limit()
options(max.print=1000000)
results <- read.csv("C:\\Users\\jake\\Dropbox\\__iKoda\\datafiles\\results.csv")
pdf(file="C:\\Users\\jake\\Dropbox\\__iKoda\\datafiles\\plots.pdf")
timesTrimmedComplete=processITStimes(results,"TSICompleted")
print(describe(timesTrimmedComplete$totaltimemins))
freq=generateQplot(timesTrimmedComplete$totaltimemins,"histogram", 1)
print(freq)
timesTrimmedINComplete=processITStimes(results,"_TSIIncomplete")
print(describe(timesTrimmedINComplete$totaltimemins))
freq1=generateQplot(timesTrimmedINComplete$totaltimemins,"histogram", 1)
print(freq1)
dev.off()
########################################################################################
generateQplot<-function(dataVector, plotType, binWidthValue)
{
freq=qplot(dataVector,geom=plotType, binwidth=binWidthValue)
return(freq)
}
processITStimes<-function(resultsData, statusCode)
{
completeResults <- resultsData[grep(statusCode, resultsData$Final_Status), ]
times <- completeResults[, grep("*duration*", colnames(completeResults))]
times[is.na(times)] <- 0
times$totaltime <- rowSums(times[,2:ncol(times)])
times$totaltimemins <-round(times$totaltime/60, digits=0)
times$rowId<-completeResults$RowId
print(statusCode);
print(describe(times$totaltimemins))
timesTrimmed<-times[times$totaltimemins<60,]
return(timesTrimmed)
}
sink()
謝謝,但我正在尋找包括兩套描述統計和二一套pdf的直方圖 – Jake
哦,我明白了。如果你想要製作一個包含多個對象的PDF文件,最好的辦法是創建一個降價文件並將其編輯爲pdf。這樣你就可以得到你的情節,你也可以從describe函數中打印出來。 對不起,我沒有更多的幫助。 –
這很有幫助! – Jake