我有原始数据,我希望看到什么样的截止水平导致截止水平以上的观察百分比.这是模拟:
data<-rnorm(100,50,30)
prop.table(table(data>10))
prop.table(table(data>20))
prop.table(table(data>30))
prop.table(table(data>40))
prop.table(table(data>50))
prop.table(table(data>60))
prop.table(table(data>70))
prop.table(table(data>80))
prop.table(table(data>90))
Run Code Online (Sandbox Code Playgroud)
这是输出:
FALSE TRUE
0.1 0.9
FALSE TRUE
0.16 0.84
FALSE TRUE
0.29 0.71
FALSE TRUE
0.36 0.64
FALSE TRUE
0.51 0.49
FALSE TRUE
0.61 0.39
FALSE TRUE
0.75 0.25
FALSE TRUE
0.86 0.14
FALSE TRUE
0.91 0.09
Run Code Online (Sandbox Code Playgroud)
但这是一种粗略而低效的方式,因为你会同意.为无限地计算每个截止值的相应百分比的Instread,我想构建一个表示该关系的图,其中X轴表示所有可能的截止水平的范围,Y轴表示从0到100的百分比.类似于这样的:
请忽略绘图的轴标签等,这只是提供一般的例子.有什么建议?
我相信您正在寻找ecdf()创建经验累积分布函数的函数.
data<-rnorm(1000,50,30)
a = ecdf(data)
plot(a)
Run Code Online (Sandbox Code Playgroud)
| 归档时间: |
|
| 查看次数: |
288 次 |
| 最近记录: |