Проблемы в вязании в PDF в R - PullRequest
Новая
       9

Проблемы в вязании в PDF в R

0 голосов
/ 15 февраля 2020

У меня есть следующий код, который запускается при выполнении блока кода: 

library(boot)
data("FishMercury")
# part a 
hist(FishMercury$Mercury)
boxplot(FishMercury$Mercury)

#part b

boot.fish = numeric(10^4)
for(i in 1: 10^4)
{
  boot.fish[i] = mean(sample(FishMercury$Mercury, 10, replace = TRUE))
}
boot.fishMean = mean(boot.fish)
boot.fishSD = sd(boot.fish)

quantile(boot.fish, .025, .0975)


#part c
FishMercury2 <- FishMercury[FishMercury < max(FishMercury)]
for (i in 1:n){boot.mean[i] <- mean(sample(FishMercury2, length(FishMercury2),replace = TRUE))}
mean(boot.mean)
sd(boot.mean)
quantile(boot.mean, prob= c(.025,.975))

Однако, когда я пытаюсь связать в PDF, я получаю ошибку: 

Error in eval(expr, envir, enclos) : object 'boot.mean' not found
Calls: <Anonymous> ... handle -> withCallingHandlers -> withVisible -> eval -> eval
Execution halted

Может кто-нибудь помочь мне решить эту проблему ??

1 Ответ

0 голосов
/ 16 февраля 2020

В #part c вашего кода отсутствуют некоторые вещи, которые приводят к ошибкам. 


#part c

FishMercury2 <- FishMercury[ FishMercury < max( FishMercury)]

# first error, you use an uninitiated object n
n <- length( FishMercury2)
# second error: you are replacing elements in boot.mean, but never initiated boot.mean.
# so, for example, create one:
boot.mean <- numeric( length = n)

for (i in 1: n){
  boot.mean[i] <- mean(sample(FishMercury2, length(FishMercury2),replace = TRUE))
  }

mean(boot.mean)
sd(boot.mean)
quantile(boot.mean, prob= c(.025,.975))
...