- 阅读权限
- 255
- 威望
- 0 级
- 论坛币
- 45705 个
- 通用积分
- 2921.8361
- 学术水平
- 408 点
- 热心指数
- 417 点
- 信用等级
- 407 点
- 经验
- 18062 点
- 帖子
- 470
- 精华
- 0
- 在线时间
- 344 小时
- 注册时间
- 2017-9-10
- 最后登录
- 2022-12-5
已卖:55份资源
讲师
还不是VIP/贵宾
TA的文库 其他... XA0小金库
充实*课堂
- 威望
- 0 级
- 论坛币
 - 45705 个
- 通用积分
- 2921.8361
- 学术水平
- 408 点
- 热心指数
- 417 点
- 信用等级
- 407 点
- 经验
- 18062 点
- 帖子
- 470
- 精华
- 0
- 在线时间
- 344 小时
- 注册时间
- 2017-9-10
- 最后登录
- 2022-12-5
 | 开心 2022-9-4 00:30:59 |
|---|
签到天数: 91 天 连续签到: 1 天 [LV.6]常住居民II
|
88论坛币
|
求助,在R中运行到最下面出现error mapping` must be created by `aes(),这是哪里出了问题,明明以前也是这样写ggplot 里的
画图,怎么这次就出错了?完全无头绪!!
~~~~~~~~~
a=c(1,2,3,1,2,3)
b=c(1,1,1,2,2,2)
y=c(27,21,33,15,6,11)
k=c(295,416,308,474,540,498)
aa=as.factor(a)
bb=as.factor(b)
options(contrasts=c("contr.sum","contr.sum"))
collator.out=glm(y~aa+bb,family = poisson,offset = log(k))
summary(collator.out)
collator.fit=predict.glm(collator.out,type = 'response',se.fit=TRUE)
y.fit=collator.fit$fit
y.sef=collator.fit$se.fit
yfL=y.fit-2*y.sef
yfH=y.fit+2*y.sef
data.jam=data.frame(a,aa,b,bb,y,y.fit,y.sef,yfL,yfH)
library(ggplot2)
w=ggplot(data.jam,mapping = aes_(a,y))
w + geom_point(data.jam, aes_(a,y,color=b))
|
|