Estoy usando la funcionprcomp Para calcular los dos primeros componentes principales. Sin embargo, mis datos tienen algunos valores de NA y, por lo tanto, la función produce un error. El na.action definido parece no funcionar a pesar de que se menciona en el archivo de ayuda?prcomp

Aquí está mi ejemplo:

d <- data.frame(V1 = sample(1:100, 10), V2 = sample(1:100, 10))

prcomp(d, center = TRUE, scale = TRUE, na.action = na.omit)

d$V1[5] <- NA
d$V2[7] <- NA

prcomp(d, center = TRUE, scale = TRUE, na.action = na.omit)

Estoy usando la versión R más nueva 2.15.1 para Mac OS X.

¿Alguien puede ver la razón mientras?prcomp falla

Aquí está mi nuevo ejemplo:

d <- data.frame(V1 = sample(1:100, 10), V2 = sample(1:100, 10))

result <- prcomp(d, center = TRUE, scale = TRUE, na.action = na.omit)

result$x

d$V1[5] <- NA

result <- prcomp(~V1+V2, data=d, center = TRUE, scale = TRUE, na.action = na.omit)

result$x

¿Es posible retener la fila 5 en PC1 y PC2? Por supuesto, en mi conjunto de datos real tengo más de dos columnas de variables y solo faltan algunas de ellas y no quiero perder la información restante oculta en los otros valores.