男性と女性の両方の被験者について、2 日間 (0 日目、1 日目) の 2 つのレベルのデータ (高、低) をプロットしたいと考えています。私は日ごと、レベルごとに面取りに成功しています。しかし、性別を組み合わせて特定することはできません。男性/女性を 0 日目と 1 日目に一緒に表示したいと思います。以下は、私が作成しようとしている特定のコードです。前もって感謝します
data <- function(ids,time_vec) {
obs.data <-
data.frame(expand.grid(ids,time_vec),DOSE=0,Conc=rnorm(13,10,2),Day=0)
names(obs.data) <- c("ID","TIME","DOSE","Conc","Day")
obs.data<-obs.data[order(obs.data$ID),]
return(obs.data)
}
test<-data(ids=1:4, time_vec= seq(0,120,10))
test$Gender<-ifelse(test$ID==1|test$ID ==3,"Male","Female")
test$Day<-ifelse(test$ID==1|test$ID==2,"Day 1","Day 2")
test$DoseLevel<-ifelse(test$ID==1|test$ID==2,"Low","High")
gf1<-ggplot(test,aes(x=TIME, y=(Conc), group=interaction(DoseLevel,Day,
Gender)))+
geom_line(size=1.25)+
facet_grid(DoseLevel~.,as.table=FALSE)
gf1
gf2<-gf1+ geom_point(aes(shape=factor(Day), fill=factor(Day),
colour=factor(Day)),size=4,show_guide=TRUE)+
scale_shape_manual(values=c(21, 21))+
scale_fill_manual(values=c("black","white"))+
scale_colour_manual(values=c("black","black"))
gf2