Hi Ana, >From your desired result, it looks like those two NA values in PLASER are the ones you want to drop. If so, try this:
b$pheno<-ifelse(b$PLASER==2 | b$FLASER==2 | is.na(b$PLASER) & b$FLASER == 2,2,1) and if I have it the wrong way round, swap FLASER and PLASER in the bit I have added. Jim On Sat, Jun 13, 2020 at 10:46 AM Ana Marija <sokovic.anamar...@gmail.com> wrote: > > Hello > > I have a data frame like this: > > > head(b) > FID IID FLASER PLASER > 1: fam1000 G1000 1 1 > 2: fam1001 G1001 1 1 > 3: fam1003 G1003 1 2 > 4: fam1005 G1005 1 1 > 5: fam1009 G1009 1 1 > 6: fam1052 G1052 1 1 > ... > > > table(b$PLASER,b$FLASER, exclude = NULL) > > 1 2 3 <NA> > 1 836 14 0 0 > 2 691 70 43 2 > 3 2 7 21 0 > <NA> 4 1 0 7 > > I am trying to make a new column "pheno" so that I reduce the number of NAs > > right now I am doing: > > > b$pheno=ifelse(b$PLASER==2 | b$FLASER==2,2,1) > > table(b$pheno, exclude = NULL) > > 1 2 <NA> > 859 828 11 > > I would like to reduce this number of NAs to be 7 > so I would like to have in "pheno column" > 7 NAs > 825 2s (825=691+14+70+7+43) > and the rest would be 1s (866=1698-7-825) > > How can I change the above command to get these numbers? > > Thanks > Ana > > ______________________________________________ > R-help@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. ______________________________________________ R-help@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.
