ho più o meno questa funzione:Non capisco errore "oggetto non trovato" all'interno della funzione
plot_pca_models <- function(models, id) {
library(lattice)
splom(models, groups=id)
}
e sto chiamando in questo modo:
plot_pca_models(data.pca, log$id)
risultati wich in questo errore :
Error in eval(expr, envir, enclos) : object 'id' not found
quando lo chiamo io, senza la funzione di involucro:
splom(data.pca, groups=log$id)
solleva questo errore:
Error in log$id : object of type 'special' is not subsettable
ma quando faccio questo:
id <- log$id
splom(models, groups=id)
si comporta come previsto.
Per favore qualcuno può spiegare perché si comporta così e come correggerlo? Grazie.
btw: Sono consapevole di domande simili qui, ad esempio:
- Help understand the error in a function I defined in R
- Object not found error with ddply inside a function
- Object disappears from namespace in function
ma nessuno di loro mi ha aiutato.
modificare: Come richiesto, non v'è piena funzione "plot_pca_models":
plot_pca_models <- function(data, id, sel=c(1:4), comp=1) {
# 'data' ... princomp objects
# 'id' ... list of samples id (classes)
# 'sel' ... list of models to compare
# 'comp' ... which pca component to compare
library(lattice)
models <- c()
models.size <- 1:length(data)
for(model in models.size) {
models <- c(models, list(data[[model]]$scores[,comp]))
}
names(models) <- 1:length(data)
models <- do.call(cbind, models[sel])
splom(models, groups=id)
}
EDIT2: sono riuscito a rendere il problema riproducibile.
require(lattice)
my.data <- data.frame(pca1 = rnorm(100), pca2 = rnorm(100), pca3 = rnorm(100))
my.id <- data.frame(id = sample(letters[1:4], 100, replace = TRUE))
plot_pca_models2 <- function(x, ajdi) {
splom(x, group = ajdi)
}
plot_pca_models2(x = my.data, ajdi = my.id$id)
che producono lo stesso errore come sopra.
'log' è una funzione in base R. È buona norma non assegnare un nome agli oggetti dopo le funzioni ... può creare confusione. – Thomas
Non chiamare il 'log' degli oggetti. 'log' è una funzione del pacchetto' base'. Cambia il nome e dicci se il problema si verifica ancora. – shadow
Questo perché reticolo utilizza una valutazione non standard. In realtà non dici quello che vuoi, ma prova 'print (splom (models, groups = id))' – hadley