2017-03-08 3 views
1

저는 R의 매우 기본적인 사용자입니다. 그래서 질문의 단순함이나 공식이 부족한 경우 미리 사과드립니다.barplot에 t 값과 신뢰 구간 추가 R

나는 하나의 연속적인 수치 변수와 두 개의 레벨이 각각있는 두 개의 계수가있는 큰 데이터 세트를 가지고 있습니다.

ExampleBarPlot <- ggplot(new.table,aes(Lexicality,ErpMinAv,fill=Frequency)) + geom_bar(stat="identity",position="dodge") + xlab("Lexicality") + ylab("Microvolts") + labs(title = "Frequency effect for singular nouns and pseudoword controls") + scale_y_continuous("Microvolts",breaks = round(seq(0, -20, by = -0.5),1)) + guides(fill=guide_legend(title="Frequency"))+ scale_colour_manual(values = c("blue","red")) 

줄거리는 다음과 같습니다

wordhigh.mu <- -2 
    wordlow.mu <- -2.5 
    pswordhigh.mu <- -1.5 
    pswordlow.mu <- -1.5 
    sigma <- 0.3 
wordshigh <- rnorm(50,mean = wordhigh.mu,sd=sigma) 
wordslow <- rnorm(50,mean = wordlow.mu,sd=sigma) 
pswordshigh <- rnorm(50,mean = pswordhigh.mu,sd=sigma) 
pswordslow <- rnorm(50,mean = pswordlow.mu,sd=sigma) 
value <- c(wordshigh,wordslow,pswordshigh,pswordslow) 
LexicalitySample <- c(rep("Word",100),rep("Pseudoword",100)) 
FrequencySample <- c(rep("High",50),rep("Low",50),rep("High",50),rep("Low",50)) 
new.table <- data.frame(ErpMinAv=value,Lexicality=LexicalitySample,Frequency=FrequencySample) 

내가 ggplot를 사용하여 내 데이터를 음모 관리 :

이 (자세한 이하) 생성/인공 데이터를 기반으로 내 데이터의 재구성이다 이 :

Plot using simulated data

내가 뭘 N을하고 싶습니다 ow는 가명 간의 빈도 차이가 크지 않음을 보여 주지만 단어 간에는 유의미합니다. 이를 위해 의미 통계 (필자의 경우에는 t 값)와 신뢰 구간을 갖는 것이 좋습니다. 이것들을 계산하는 방법을 알고 있지만, 나는 그것을 barplot에 추가하는 방법을 모른다.

나는 인터넷에서 광범위하게 보았지만 내가보고 싶은 결과를 얻은 예를 찾을 수 없었다.

모든 도움을 주시면 감사하겠습니다.

답변

3

나는 평균과 신뢰 구간을 원한다고 가정 할 것입니다.

현재 누적 막대 그래프가 있기 때문에 합계를 플로팅합니다.

enter image description here

우리는 수단을 계산하는 stat_summary()을 사용할 수 있으며, 부트 스트랩 신뢰 구간 : 우리는 테두리 색상을 추가 할 때 볼 수있는 응답에 대한

ggplot(new.table, aes(Lexicality,ErpMinAv,fill=Frequency)) + 
    stat_summary(geom = 'bar', fun.y = mean, position = position_dodge(0.9)) + 
    stat_summary(geom = 'errorbar', fun.data = mean_cl_boot, position = position_dodge(0.9), 
       width = 0.5) + 
    scale_y_continuous("Microvolts",breaks = round(seq(0, -20, by = -0.5),1)) 

enter image description here

+0

감사합니다, 그리고에 대한 나를 stat_summary에 소개해 줬어. 굉장히 유용하다! – HernanLG