RT，
代码如下：x1<-c(2.8,2.5,2.9,2.7,2.6,2.8,3.1,3.0,2.5,2.4,2.0,2.1,1.9,1.8,2.2,2.0,2.1,2.2,2.1,1.8,2.4,2.0,2.2,2.1,1.9,1.8,2.0,2.1,2.2,2.0,2.0,2.3,1.6,1.8,2.2,2.3,1.9,1.5,2.0,2.1,2.8,2.6,2.1,2.8,1.9,2.2,2.1,1.8,1.7,2.1,2.4,1.5,2.1,1.8,1.5,2.0,2.2,1.8,2.7,1.8,2.5,2.4,2.3,2.0,2.1,2.5,2.0,2.4,2.1,2.1,2.0,2.1,2.3,2.2,2.5,2.4,2.3,2.3,2.5,2.5,2.4,2.5,2.4,2.1,2.5,2.4,2.6,2.3,2.2,2.2,1.6,1.6,2.0,2.0,1.8,2.0,2.5,2.3,2.1,2.4,2.1,2.1,2.2,2.1,2.1,2.1,2.1,2.2,2.1,2.3,2.2,2.2,2.1,2.1,2.1,2.3,2.2,2.2,2.2,2.1)
x2<-c(10.5,10.5,10.3,10.3,10.1,10.7,10.6,11,10.6,10.2,11,10.8,10.3,10,10.5,10,10.6,9.3,10.8,9.5,10.5,9.2,9,8,7,6.5,8.5,9.5,9.5,12.5,8.8,10.8,10.8,12.8,10.6,11.6,11.6,10.4,10,10,10.4,10,10.2,11.5,10.5,11.7,11.7,10.2,10.2,10.2,13.2,8.1,12.1,8.1,10.9,8.9,7.9,7.5,6.5,7.5,8.5,10.5,9.5,10.5,8.5,10.5,9.3,10.3,8,9,12,9,8.8,9.8,8.8,11.8,6.9,6.9,7.9,8.9,7.9,7.9,8.2,8.2,7.4,7.4,7.7,7.7,8.4,7.8,10.7,11.7,11.7,10,10,10,8.3,8.3,8.3,10.3)
x3<-c(2.5,2.6,3.5,2.5,3.0,3.0,2.7,2.7,3.3,3.3,2.6,2.7,2.9,2.9,3.0,2.8,2.8,2.7,2.7,3.0,3.2,3.7,3.3,3.3,3.0,2.8,3.0,2.7,2.7,2.9,2.7,2.7,2.8,2.8,2.6,2.8,2.9,2.9,2.6,3.0,2.4,2.4,2.5,2.5,2.5,2.5,2.4,2.4,2.4,2.5,2.4,2.3,2.1,2.6,2.8,2.9,2.2,2.5,2.5,2.7,2.5,2.4,2.7,2.0,2.4,2,2,2.3,2.3,1.8,3.5,3.7,3.8,3.6,3.6,3.2,3.2,3.4,3.4,3.4,3,2.8,2.6,2.4,3,2.8,2.7,2.6,2.8,2.7,3.2,3.4,3.4,3.3,3.4,3.4,3.4,3.4,3.5,3.7)
x4<-c(12.5,11.2,11.2,8.5,8.5,11.5,10.5,9.5,8.5,8.5,10.5,10.5,10.5,10.5,10.7,10.3,10.3,10.3,10.2,10.2,11.2,10.4,10.4,10.4,10.4,9.6,10.6,11.6,10.6,10.6,10.8,11.8,11.8,10.8,10.7,12.7,10.7,10.7,10.8,10.8,10.5,10.5,10.3,10.3,9.3,10.3,10.2,10.2,10.2,10.2,10.1,10.1,10.1,10.1,10.2,10.2,10.2,10.3,10.3,10.4,10.4,10.4,10.4,10.6,10.6,10.6,10.7,10.7,10.7,9.1,9.9,9.9,10,10,9,9,9,9,10,10,10.9,10.9,10.9,10,12,12,11,10.8,10.8,9.8,10.8,11.3,11.3,10.3,11.3,11.3,10.3,11.1,11.1,11.1)
x<-c(x1,x2,x3,x4)
hist(x,prob=T,breaks=60,xlim = c(0,15), ylim=c(0,0.4),col="gray",main="Kernel Density of stiffnes Per cells",xlab="module Young")
lines(density(x,bw=0.8),col="red",lwd=3)
还想问下，核密度函数，有可以导出的参数么，如何导出，比如两个峰的中间值，分离度
|
发表于 2014-6-25 21:49:57