2012-11-19 5 views
4

내일 예정인 과제에 대한 신뢰 구간을 가진 상호 작용 플롯을 그려야합니다. 이것은 내 데이터 세트입니다 :신뢰 구간이있는 상호 작용 플롯

Snail0 

    water lightness total_migration 
    control dark 9.5 
    control dark 20 
    control dark 22 
    control dark 24 
    control dark 10 
    control dark 18 
    control dark 19 
    control dark 19.5 
    control dark 23 
    control dark 2 
    control dark 10.5 
    control dark 22.5 
    control dark 26 
    control dark 32 
    control dark 2 
    control dark 12.5 
    control dark 13.5 
    control dark 16.5 
    control dark 19 
    control dark 20 
    control dark 21 
    control dark 22 
    control dark 24 
    control dark 25 
    control light 11 
    control light 20.5 
    control light 22.5 
    control light 24 
    control light 21 
    control light 23 
    control light 23.5 
    control light 28 
    control light 30 
    control light 18.5 
    control light 19.5 
    control light 19.5 
    control light 23.5 
    control light 24 
    control light 2 
    control light 20 
    control light 23 
    control light 28.5 
    control light 29 
    control light 16 
    control light 17.5 
    control light 26.5 
    control light 28.1 
    control light 30.5 
    crab light 6 
    crab light 14 
    crab light 15 
    crab light 17 
    crab light 19 
    crab light 7 
    crab light 14 
    crab light 17 
    crab light 22 
    crab light 26 
    crab light 4.5 
    crab light 10 
    crab light 11 
    crab light 21 
    crab light 24 
    crab light 14 
    crab light 20 
    crab light 22.5 
    crab light 23 
    crab light 14 
    crab light 16 
    crab light 18.5 
    crab light 22 
    crab light 23 
    crab dark 14 
    crab dark 19 
    crab dark 20 
    crab dark 23 
    crab dark 25 
    crab dark 9 
    crab dark 18 
    crab dark 18 
    crab dark 20 
    crab dark 21 
    crab dark 15 
    crab dark 15 
    crab dark 15 
    crab dark 19 
    crab dark 23 
    crab dark 12 
    crab dark 14.5 
    crab dark 21 
    crab dark 21.5 
    crab dark 22 
    crab dark 14.5 
    crab dark 16 
    crab dark 19 
    crab dark 20 
    crab dark 22.5 
    snail dark 6 
    snail dark 15 
    snail dark 20 
    snail dark 22 
    snail dark 15 
    snail dark 15 
    snail dark 16 
    snail dark 17 
    snail dark 21 
    snail dark 15 
    snail dark 20 
    snail dark 15 
    snail dark 16 
    snail dark 19.5 
    snail dark 23 
    snail dark 24 
    snail dark 4 
    snail dark 18.5 
    snail dark 20 
    snail dark 24 
    snail dark 25 
    snail light 13 
    snail light 13 
    snail light 13 
    snail light 18 
    snail light 15 
    snail light 16 
    snail light 17 
    snail light 17 
    snail light 19 
    snail light 14 
    snail light 14 
    snail light 19 
    snail light 7 
    snail light 13 
    snail light 13 
    snail light 13 
    snail light 16 
    snail light 3 
    snail light 9 
    snail light 15 
    snail light 22 
    snail light 23 

interaction.plot(Snail0$water,Snail0$lightness,Snail0$total_migration) 

내가 여기 ToothGrowth data과 하단에 완료 어떤 식으로 뭔가를 시도하고있다,하지만 난 분명히 내가 오류 메시지 만 얻을으로 뭔가 잘못하고있는 중이 야 .

Snail0$water <- factor(Snail0$water) 


df<-with(Snail0, aggregate(total_migration, list(lightness=lightness, water=water), mean)) 


df$se<-with(Snail0, aggregate(total_migration, list(lightness=lightness, water=water), function(x) sd(x)/sqrt(10)))[,3] 


gp<-ggplot(df, aes(x=water, y=total_migration, colour=lightness, group=lightness)) 

gp+geom_line(aes(linetype=lightness)) + geom_point(aes(shape=lightness)) + geom_errorbar(aes(ymax=x+se, ymin=x-se)) 

오류 메시지 :

Error: Aesthetics must either be length one, or the same length as the dataProblems:lightness, water, lightness, lightness 

분명히, 내가 뭘하는지 아무 생각이 없다, 그래서 사람이 내 코드를 수정 도와주세요 수 있습니까?

+1

는 http://stats.stackexchange.com/questions/9477/how-to- 참조 : 당신이 ggplot 명령에 xtotal_migration를 교체 할 경우

따라서, 귀하의 코드가 작동 그림 그리기 - 상호 작용 - 작의 - 신뢰 구간 함께 –

답변

3

aggregate 명령으로 개체 df을 만들면 종속 변수를 포함하는 열은 더 이상 total_migration이 아니라 x이됩니다.

gp <- ggplot(df, aes(x=water, y=x, colour=lightness, group=lightness)) 
gp + geom_line(aes(linetype=lightness)) + 
    geom_point(aes(shape=lightness)) + 
    geom_errorbar(aes(ymax=x+se, ymin=x-se)) 

enter image description here

+0

정말 고마워요! 너는 내 하루를 보냈다. :) – SunRoad