Run Code
|
API
|
Code Wall
|
Users
|
Misc
|
Feedback
|
Login
|
Theme
|
Privacy
|
Patreon
Bootstrap Method
set.seed(5610) x <- rnorm(300,10,15) bootstrap <- function(B,data){ sam = B dat = data n = length(data) mean.vector = c() for( i in 1:sam){ a = sample(dat,n,replace = T) mean.vector[i] = mean(a) } hist(mean.vector,col = 'lightblue',freq = F,border = 'white') abline(v = mean(mean.vector), col = 'red3', lty = 2, lwd = 2) cat("\n The average estimate is: ", mean(mean.vector)) cat("\n Standard error: ", sqrt(var(mean.vector))) } mean(x) bootstrap(1000,x)
run
|
edit
|
history
|
help
0
Please
log in
to post a comment.
Ggplot
ceiling.r
AjusteLinearizado
Peso hora
teste
Multiple barplot
First blood
Ex15-12-20
Delete True Values
16-09-2020ExVitorc
Please log in to post a comment.