R 绘制趋势性检验结果图 实例操作

biostatistic

data1和data2是趋势性检验所得结果。
根据结果画图说明

1. 
   
2. data1<-sample(1:50,18)
   
3. data2<-sample(51:10,18)
   
4. 
   
5. data1<- matrix(data1,nrow=3)
   
6. data2<- matrix(data2,nrow=3)
   
7. 
   
8. datab<-round(sapply(1:6,function(x) data1[,x]/sum(data1[,x])),digits=2)*100
   
9. ldatab<-round(sapply(1:6,function(x) data2[,x]/sum(data2[,x])),digits=2)*100
   
10. 
    
11. colnames(datab)<-rep(c("非胖子","胖子"),3)
    
12. colnames(ldatab)<-rep(c("非胖子","胖子"),3)
    
13. 
    
14. 
    
15. par(mfrow=c(2, 3))
    
16. barplot(datab[,1:2],sub ="P=0.085 for trend",space=1,main="腰围",axes=F);
    
17. axis(1,labels=F,at=c(0,9))
    
18. axis(2,at=seq(from=0,to=100,by=20),labels=c("0%","20%","40%","60%","80%","100%"))
    
19. 
    
20. barplot(datab[,3:4],sub ="P=0.81 for trend",space=1,main="臀围",axes=F);
    
21. axis(1,labels=F,at=c(0,9))
    
22. axis(2,at=seq(from=0,to=100,by=20),labels=c("0%","20%","40%","60%","80%","100%"))
    
23. 
    
24. barplot(datab[,5:6],sub ="P=0.0321 for trend",space=1,main="体重",axes=F);
    
25. axis(1,labels=F,at=c(0,9))
    
26. axis(2,at=seq(from=0,to=100,by=20),labels=c("0%","20%","40%","60%","80%","100%"))
    
27. 
    
28. barplot(ldatab[,1:2],sub ="P=0.0091 for trend",space=1,main="腰围",axes=F);
    
29. axis(1,labels=F,at=c(0,9))
    
30. axis(2,at=seq(from=0,to=100,by=20),labels=c("0%","20%","40%","60%","80%","100%"))
    
31. 
    
32. barplot(ldatab[,3:4],sub ="P=0.0501 for trend",space=1,main="臀围",axes=F);
    
33. axis(1,labels=F,at=c(0,9))
    
34. axis(2,at=seq(from=0,to=100,by=20),labels=c("0%","20%","40%","60%","80%","100%"))
    
35. 
    
36. barplot(ldatab[,5:6],sub ="P=0.0101 for trend",space=1,main="体重",axes=F);
    
37. axis(1,labels=F,at=c(0,9))
    
38. axis(2,at=seq(from=0,to=100,by=20),labels=c("0%","20%","40%","60%","80%","100%"))
    
39. 
    

复制代码

以上数据纯属随机所得，只做实践例子参考。不代表真是结果。

停在八楼的35路

学习了，谢谢b版的实例，俺们这些新手就喜欢看这个，嘎嘎，先MARK一个，慢慢看

zcvb6

data1和data2是趋势性检验所得结果这个怎么来理解啊