使 table/matrix 的概率密度和相关的中断
Make table/matrix of probability densities and associated breaks
以下代码采用向量 V1，并对由 V1 组成的随机正态样本进行 10000 次自举，用结果创建一个包含 10000 列的矩阵。然后为该矩阵创建直方图。
V1 <- c(0.18, 0.2, 0.24, 0.35, -0.22, -0.17, 0.28, -0.28, -0.14, 0.03, 0.87, -0.2, 0.06, -0.1, -0.72, 0.18, 0.01, 0.31, -0.36, 0.61, -0.16, -0.07, -0.13, 0.01, -0.09, 0.26, -0.14, 0.08, -0.62, -0.2, 0.3, -0.21, -0.11, 0.05, 0.06, -0.28, -0.27, 0.17, 0.42, -0.05, -0.15, 0.05, -0.07, -0.22, -0.34, 0.16, 0.34, 0.1, -0.12, 0.24, 0.45, 0.37, 0.61, 0.9, -0.25, 0.02)
xxx <- round(sapply(1:10000, function(i) rnorm(length(V1), mean=sample
(V1, length(V1), replace=TRUE), sd(V1))),2)
h <- hist(xxx, plot=T)
我想以矩阵或 table 格式创建它的概率密度函数的打印输出，即得到一个矩阵 1.0, 0.9, 0.8, 0.7, 0.6, 0.5, 0.4, 0.3, 0.2, 0.1, 0, -0,1 -0.2, -0.3, -0.4, -0.5, -0.6, -0.7, -0.8, -0.9, -1.0 作为中断，相关的概率密度在下一列。
我的问题有两个。首先，指定我想要的休息时间失败了。其次，用 h$breaks 和 h$density 制作矩阵也失败了。任何见解将不胜感激。谢谢。
#This works, but I want to specify the breaks
h <- hist(xxx, plot=T, breaks=20)
#This gives error "some 'x' not counted; maybe 'breaks' do not span range of 'x'"
h <- hist(xxx, plot=T, breaks=c(1.0, 0.9, 0.8, 0.7, 0.6, 0.5, 0.4, 0.3, 0.2, 0.1, 0, -0,1 -0.2, -0.3, -0.4, -0.5, -0.6, -0.7, -0.8, -0.9, -1.0))
#This gives error number of rows of result is not a multiple of vector length
ddd<- cbind(h$breaks, h$density)
一开始你有一个类型错误，你混淆了“.”和“，”。
第二个休息点必须像这样跨越整个观察范围：
h <- hist(xxx, plot=F, breaks=c(max(xxx),1.0, 0.9, 0.8, 0.7, 0.6, 0.5, 0.4, 0.3, 0.2, 0.1, 0, -0.1, -0.2, -0.3, -0.4, -0.5, -0.6, -0.7, -0.8, -0.9, -1.0,min(xxx)))
最后，您想选择“中音”而不是间歇：
ddd<- cbind(h$mids, h$density)
以下代码采用向量 V1，并对由 V1 组成的随机正态样本进行 10000 次自举，用结果创建一个包含 10000 列的矩阵。然后为该矩阵创建直方图。
V1 <- c(0.18, 0.2, 0.24, 0.35, -0.22, -0.17, 0.28, -0.28, -0.14, 0.03, 0.87, -0.2, 0.06, -0.1, -0.72, 0.18, 0.01, 0.31, -0.36, 0.61, -0.16, -0.07, -0.13, 0.01, -0.09, 0.26, -0.14, 0.08, -0.62, -0.2, 0.3, -0.21, -0.11, 0.05, 0.06, -0.28, -0.27, 0.17, 0.42, -0.05, -0.15, 0.05, -0.07, -0.22, -0.34, 0.16, 0.34, 0.1, -0.12, 0.24, 0.45, 0.37, 0.61, 0.9, -0.25, 0.02)
xxx <- round(sapply(1:10000, function(i) rnorm(length(V1), mean=sample
(V1, length(V1), replace=TRUE), sd(V1))),2)
h <- hist(xxx, plot=T)
我想以矩阵或 table 格式创建它的概率密度函数的打印输出，即得到一个矩阵 1.0, 0.9, 0.8, 0.7, 0.6, 0.5, 0.4, 0.3, 0.2, 0.1, 0, -0,1 -0.2, -0.3, -0.4, -0.5, -0.6, -0.7, -0.8, -0.9, -1.0 作为中断，相关的概率密度在下一列。
我的问题有两个。首先，指定我想要的休息时间失败了。其次，用 h$breaks 和 h$density 制作矩阵也失败了。任何见解将不胜感激。谢谢。
#This works, but I want to specify the breaks
h <- hist(xxx, plot=T, breaks=20)
#This gives error "some 'x' not counted; maybe 'breaks' do not span range of 'x'"
h <- hist(xxx, plot=T, breaks=c(1.0, 0.9, 0.8, 0.7, 0.6, 0.5, 0.4, 0.3, 0.2, 0.1, 0, -0,1 -0.2, -0.3, -0.4, -0.5, -0.6, -0.7, -0.8, -0.9, -1.0))
#This gives error number of rows of result is not a multiple of vector length
ddd<- cbind(h$breaks, h$density)
一开始你有一个类型错误，你混淆了“.”和“，”。
第二个休息点必须像这样跨越整个观察范围：
h <- hist(xxx, plot=F, breaks=c(max(xxx),1.0, 0.9, 0.8, 0.7, 0.6, 0.5, 0.4, 0.3, 0.2, 0.1, 0, -0.1, -0.2, -0.3, -0.4, -0.5, -0.6, -0.7, -0.8, -0.9, -1.0,min(xxx)))
最后，您想选择“中音”而不是间歇：
ddd<- cbind(h$mids, h$density)