Erreur dans Si la déclaration dans R

Question

ci-dessous est le code. Il semble que je n'ai pas réussi à définir les vecteurs 'nobs', 'cor' correctement et l'instruction if essaye d'évaluer une valeur NA. Quelqu'un peut-il me signaler mon erreur? Merci, Tamir

corr<-function(directory, threshold=0){ 

allfiles<-list.files("specdata", full.names = TRUE) 
NumOfFiles<-length(allfiles) 

n<-1 
id<-vector() 
nobs<-vector() 
cor<-vector() 

for (i in NumOfFiles){ 
    data<-read.csv(allfiles[i]) 
    cleandata<-data[complete.cases(data),] 
    id[n]=i 
    nobs[n]<-sum(complete.cases(data)) #Number of completre cases 
    n<-n+1 


if (nobs[n]>threshold){ 
    cor[n]<-cor(cleandata[sulfate],cleandata[nitrate]) 
    } 
} 

return(cor) 

}

Answer 1

Le code était pas loin. J'ai changé la position de n <- n+1. Dans votre exemple, il incrémentait trop tôt dans la fonction. Lorsque la fonction est descendue à if(nobs[n]>threshold), elle vérifiait APRÈS avoir fait n + 1, donc c'était toujours NA. Deuxièmement, je mets des guillemets autour des noms de polluants.

corr<-function(directory, threshold=0){ 

allfiles<-list.files("specdata", full.names = TRUE, pattern="csv") 
NumOfFiles<-length(allfiles) 

n<-1 
id<-vector() 
nobs<-vector() 
cor<-vector() 

for (i in 1:NumOfFiles){ 
    data<-read.csv(allfiles[i]) 
    cleandata<-data[complete.cases(data),] 
    id[n]=i 
    nobs[n]<- sum(complete.cases(data)) #Number of completre cases 

    if (nobs[n]>threshold){ 
     cor[n]<-cor(cleandata[,"sulfate"],cleandata[,"nitrate"]) 
     n<-n+1 
     } 

} 

return(cor) 
} 

Nous pouvons tester avec:

summary(corr("specdata")) 
    Min. 1st Qu. Median  Mean 3rd Qu.  Max. 
-1.00000 -0.05205 0.10840 0.13800 0.27890 1.00000 

cr <- corr("specdata", 150) 
head(cr) 
[1] -0.01895754 -0.14051254 -0.04389737 -0.06815956 
[5] -0.12350667 -0.07588814 

qui correspondent à la fois la sortie désirée.