Hi Z,
I guess it should work.
The apply(geno1,1,...) applies on rows 1:nrow(geno1). In the example dataset,
there were 5 rows, the result vector (tmp) is also of length 5.
A.K.
________________________________
From: JiangZhengyu <zhyjiang2...@hotmail.com>
To: smartpink...@yahoo.com
Sent: Friday, September 28, 2012 10:55 PM
Subject: RE: [R] Errors in if statement
Thanks, Arun. The matrix format & put "NA" with 0 help. but how do I do all
the rows. Is the loop still OK?
Z
> Date: Fri, 28 Sep 2012 18:13:39 -0700
> From: smartpink...@yahoo.com
> Subject: Re: [R] Errors in if statement
> To: zhyjiang2...@hotmail.com
> CC: r-help@r-project.org
>
> Hi,
> Try this:
> geno<-read.table(text="
> P1 P2 P3 P4 P5
> 1 2 2 3 2
> 2 2 2 1 1
> 1 2 1 2 NA
> NA 2 3 4 5
> 1 1 3 1 3
> ",sep="",header=TRUE,stringsAsFactors=FALSE)
> geno1<-as.matrix(geno)
> geno1[is.na(geno1)]<-0
> tmp<-apply(geno1,1,function(x) ifelse((sum(x!=2)>3) & (sum(x==1)>=1) &
> (sum(x==3)>=1), 1,0) )
> tmp
> #[1] 0 0 0 0 1
> A.K.
>
>
>
>
> ----- Original Message -----
> From: JiangZhengyu <zhyjiang2...@hotmail.com>
> To:
> Cc: r-help@r-project.org
> Sent: Friday, September 28, 2012 4:16 PM
> Subject: [R] Errors in if statement
>
>
> Hi guys, I have many rows (>1000) and columns (>30) of "geno" matrix. I use
> the following loop and condition statement (adapted from someone else code).
> I always have an error below. I was wondering if anyone knows what's the
> problem & how to fix it.
> Thanks,Zhengyu ########### geno matrix P1 P2 P3 P4
> 1 2 2 3 2
> 2 2 2 1 1
> 1 2 1 2 NANA 2 3 4 5 ###########
> for(i in 1:4) {
> cat(i,"")
> if(sum(geno[i,]!=2)>3 && sum(geno[i,]==1)>=1 && sum(geno[i,]==3)>=1){
> tmp = 1
> }
> } ########### 1 2 Error in if (sum(geno[i, ] != 2) > 3 && sum(geno[i, ] == 1)
> >= 1 && sum(geno[i, :
> missing value where TRUE/FALSE needed
>
> [[alternative HTML version deleted]]
>
> ______________________________________________
> R-help@r-project.org mailing list
> 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.
>
______________________________________________
R-help@r-project.org mailing list
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.
