[R] NA command in a 'for' loop

Rui Barradas ru|pb@rr@d@@ @end|ng |rom @@po@pt
Mon Apr 20 08:05:37 CEST 2020


Hello,

Instead of

d[d == 0] <- NA

try

d[] <- lapply(d, function(x) {is.na(x) <- x == 0; x})


Also, in the first for loop

paste(i, sep = "")

does nothing, it's the same as i.
And the same for

(d2$V4 == 1) == TRUE

Since (d2$V4 == 1)  already is FALSE/TRUE there is no need for

(.) == TRUE


Hope this helps,

Rui Barradas



Às 20:52 de 19/04/20, Helen Sawaya escreveu:
> Dear R experts,
> 
> I am using a 'for' loop to apply commands to multiple datasets (each file is one participant). The only one not working is the command that identifies zeros in my datasets and changes them to NAs. But when I look at the output, zeros ("0") are still present. Surprisingly, the functions work fine when I apply them to a single dataset (outside the loop). I've tried:
> 
> all.files <- list.files(".")
> txt.files <- grep("threat.txt",all.files,value=T)
> 
> for(i in txt.files){
>    d <- read.table(paste(i,sep=""),header=F)
>    d[d==0] <- NA #replace zeros with NA
>    write.table(d, paste0(i,".tlbs.txt"), quote=FALSE, row.names=TRUE)}
>    d<-d[ ,-c(10,11)]
>    d2<-d[complete.cases(d), ]
>    d2$V4<-as.numeric(d2$V4)
>    congruent <- (d2$V4 == 1) == TRUE
>    x <- get_tlbs(d2$V14, congruent, prior_weights = NULL, method = "weighted", fill_gaps = FALSE)
>    write.table(x, paste0(i,".tlbs.txt"), quote=FALSE, row.names=TRUE)}
> 
> I've also tried:
> 
> for(i in txt.files){
>    d <- read.table(paste(i,sep=""),header=F)
>    if (0 %in% d)
>    {replace_with_na(d,replace = list(x = 0))} # replace zeros with NA
>    d<-d[ ,-c(10,11)]
>    d2<-d[complete.cases(d), ]
>    d2$V4<-as.numeric(d2$V4)
>    congruent <- (d2$V4 == 1) == TRUE
>    x <- get_tlbs(d2$V14, congruent, prior_weights = NULL, method = "weighted", fill_gaps = FALSE)
>    write.table(x, paste0(i,".summaryoutput.txt"), quote=FALSE, row.names=TRUE)}
> 
> Thank you for your help.
> Sincerely
> Helen
> 
> 	[[alternative HTML version deleted]]
> 
> ______________________________________________
> R-help using r-project.org mailing list -- To UNSUBSCRIBE and more, see
> https://stat.ethz.ch/mailman/listinfo/r-help
> PLEASE do read the posting guide http://www.R-project.org/posting-guide.html
> and provide commented, minimal, self-contained, reproducible code.
>



More information about the R-help mailing list