h**********6
发帖数: 66 | 1 如题。本人不是统计专业，只会一些简单R，看到paper上用bootstrap获得95%CI，但不
知道怎么编程，，，若有牛人发个简单的R script例子，不胜感激（仅做学习参考之用
）
Email：t****************[email protected] 非常非常感谢 ：） | j*******2
发帖数: 309 | 2 try calibrate function in {rms} package.
calibrate(fit,B=N) | p********r
发帖数: 1465 | 3 for (i in 1:1000){
...
write.table(a, file = "",quote=FALSE,append=TRUE,row.names=FALSE,col.names=
FALSE)
}
quantile(a,0.975)
quantile(a,0.025) | h**********6
发帖数: 66 | 4 关键是resampling部分怎么写呢？
【在 p********r 的大作中提到】
: for (i in 1:1000){
: ...
: write.table(a, file = "",quote=FALSE,append=TRUE,row.names=FALSE,col.names=
: FALSE)
: }
: quantile(a,0.975)
: quantile(a,0.025)
| p********r
发帖数: 1465 | 5 说出你要做的东西。
【在 h**********6 的大作中提到】
: 关键是resampling部分怎么写呢？
| h**********6
发帖数: 66 | 6 假设A，B follow一个函数关系，对于每一个固定的B（已知），有5个A值与之对应，这
5个A假设follow normal distribution，然后我要通过A的varibility去估计C的var，C
和A follow以下关系： C=∑A*B。要用bootstrap估计C的confidence interval。
谢谢~~
【在 p********r 的大作中提到】
: 说出你要做的东西。
| p********r
发帖数: 1465 | 7 write.table(c,"")
for (i in 1:1000){
a <- matrix(rnorm(5,0,1),nrow=1,ncol=5)
b <- matrix(1:5)
c <- a %*% b
write.table(c, file = "",quote=FALSE,append=TRUE,row.names=FALSE,col.names=
FALSE)
}
test <- read.table("", header=T)
c <- test$c
quantile(c,0.975)
quantile(c,0.025) | h**********6
发帖数: 66 | 8 这个是simulation？bootstrapping不是要从原有data（就是那些A）中做resampling么？
【在 p********r 的大作中提到】
: write.table(c,"")
: for (i in 1:1000){
: a <- matrix(rnorm(5,0,1),nrow=1,ncol=5)
: b <- matrix(1:5)
: c <- a %*% b
: write.table(c, file = "",quote=FALSE,append=TRUE,row.names=FALSE,col.names=
: FALSE)
: }
: test <- read.table("", header=T)
: c <- test$c
| p********r
发帖数: 1465 | 9 B是不是random的？
么？
【在 h**********6 的大作中提到】
: 这个是simulation？bootstrapping不是要从原有data（就是那些A）中做resampling么？
| g**********t
发帖数: 475 | | h**********6
发帖数: 66 | 11 B不是
【在 p********r 的大作中提到】
: B是不是random的？
:
: 么？
| p********r
发帖数: 1465 | 12 bs <- NULL
for (i in 1:5){
a <- matrix(rnorm(5,0,1),nrow=1,ncol=5)
b <- matrix(1:5)
c[i] <- a %*% b}
for (i in 1:1000){
bs.sample <- sample(c, length(c), replace=TRUE)
bs[i] <- median(bs.sample, na.rm=TRUE)
}
quantile(bs,0.975)
quantile(bs,0.025) | h**********6
发帖数: 66 | 13 Thank you so much!
【在 p********r 的大作中提到】
: bs <- NULL
: for (i in 1:5){
: a <- matrix(rnorm(5,0,1),nrow=1,ncol=5)
: b <- matrix(1:5)
: c[i] <- a %*% b}
: for (i in 1:1000){
: bs.sample <- sample(c, length(c), replace=TRUE)
: bs[i] <- median(bs.sample, na.rm=TRUE)
: }
: quantile(bs,0.975)
|
|
