замените NA на укороченные нормальные значения распределения в r

Question

Я пытаюсь заменить NAs на укороченные нормальные значения распределения. Сначала я использовал sample следующим образом и функция работала:

v.new <- replace(vector,v, sample(8,length(v),replace =FALSE)) 

Однако, когда я пытаюсь использовать rtnorm, кажется, не работает. Я получил сообщения об ошибках, и для замены НС требуются возрастные периоды на нужный интервал. Любое предложение сделать эту работу?

library(msm) 

# Some data 
data("airquality") 
airquality$Ozone 

# My function 
add.trunc.to.NAvector <- function(vector){ 
v <- NULL 
for(i in 1:length(vector)){ 
if(is.na(vector[i])==TRUE) 
    v <- append(v, i) 
} 
mean.val <- mean(vector) 
sd.val <- sd(vector) 
min.val <- mean.val - 4 * sd.val 
max.val <- mean.val + 4 * sd.val 
v.new <- replace(vector,v, rtnorm(length(v), lower = min.val, upper = max.val)) 


return(v.new) 
} 

Answer 1

Не должно ли это работать?

v <- airquality$Ozone 
v.new <- v 
indices <- which(is.na(v)) 
m <- mean(v[-indices]) 
s <- sd(v[-indices]) 
v.new[indices] <- rtnorm(length(indices), lower = m-4*s, upper = m+4*s)