rm(list =ls())
setwd("C:\\Users\\17657\\Desktop\\ska3")
options(repos='http://cran.rstudio.com/')
library(xlsx)
baseline<- read.xlsx('临床信息.xlsx',sheetIndex = 1,header = T,encoding = 'UTF-8')
baseline<-baseline[1:1096,]
base<-subset(baseline,select = c(age,T,M,N,stage,ER,PR,HER2,subtype))
Freq<-lapply(base[,1:9],table)
prop<-lapply(Freq[1:9],prop.table)
prop
character<-c(names(Freq[1]),names(Freq[[1]]))
Noc<-c(NA,paste0(Freq[[1]],'(',prop[[1]],')'))
char<-NULL
characteristics<-data.frame('characteristics'=character,' number of class'=Noc)
for (i in 1:9) {character<-c(names(Freq[i]),names(Freq[[i]]))
Noc<-c(NA,paste0(Freq[[i]],'(',prop[[i]]*100,')'))
characteristics<-data.frame('characteristics'=character,' number of class(%)'=Noc)
char<-rbind(char,characteristics)
}
write.xlsx(char,'基线表.xlsx',col.names=T,row.names=F,showNA=F)
好啦可以啦~
q1, 首先要确定是barplot还是hist，如果是barplot的话，应该不存在breaks的问题，因为barplot的传入参数是个矩阵；我假设你要画的是个hist，我偶遇过这个问题，我的理解是hist的breaks的值要能被范围整除才行；比如x=1:200，break=7的话，就只能画出4个柱来，但如果breaks=10就没问题；基本上是这样的，偶尔也有例外；比如break=5就不行....奇怪得很
最后，没办法的办法，就只能用barplot代替hist了，barplot肯定不会有这个问题，统计下hist参数中的分布情况，转换成矩阵，用barplot吧；
q2, 貌似一般都用一组因素把这些类别区分开，我用abcde，表示你的小学，中学...了，比如这样：
a=1:7b=8:10c=c(9,10,11)d=c(40,55)e=100:110f=factor(c(rep(1,sum(length(a),length(b),length(c))),rep(2,sum(length(d),length(e)))))#先用c()生成数组，在转换成factor，其实数组也ok的，不过plot()中两个数组和factor不一样 x=c(a,b,c,d,e)plot(x~f)q3, 就我所知不行；yes或no一定也要是能映射到x，y范围内的点才行；你是想表示分类结果吗？如果是的话，通常用颜色，或者在点旁边的text表示。
q4, 举个例子吧
x=-50:50y=x^2+x+1z=10*abs(x)+1 plot(x,y,type='l')lines(x,z,lty=3)legend(c('type1','type2'), x=-20,y=2500, col=c('black','red'), lty=c(1,3))legend的x和y是legend的左上角，匿名参数是类型名称，col,lty,pch 是对应的颜色，线类型，和点类型。
最后，我现在多用ggplot2，如果不抵触的话可以看看，和R的基础作图包思路不是很一样，但是图很清新的；
如果还有问题，建议把数据集data.frame粘贴几行上来，我也试试；
你先写一个函数,然后在R的界面加载新建一个脚本:
fun<-function(x){
if(x>=0 &x<=3) y=3*x+2
else y=2*x-0.5*x^2
y
}
保存为fun.R
然后在R界面里
>source("fun.R")
>x<-seq(0,6,0.01) 每隔0.01取一个点
>y<-fun(x)
>plot(x,y)
