cut2 error
On Oct 17, 2012, at 1:52 PM, Charles Determan Jr wrote:
Hi A.K. I tried your code exactly as you presented but I only get the following output, did I somehow miss something? Was there something else loaded? Thanks, [1] [96, 99) [96, 99) [96, 99) [96, 99) [96, 99) 99 [96, 99) [96, 99) [9] 99 [96, 99) 99 [96, 99) [96, 99) [96, 99) [96, 99) [96, 99) [17] 99 [96, 99) [96, 99) [96, 99) [96, 99) 99 [96, 99) [96, 99) [25] 99 99 [96, 99) 99 [96, 99) 99 99 99 [33] 99 [96, 99) 99 [96, 99) 99 [96, 99) [96, 99) 99 [41] [96, 99) [96, 99) 99 99 [96, 99) 99 99 [96, 99) [49] [96, 99) [96, 99) 99 99 [96, 99) [96, 99) [96, 99) [96, 99) [57] [96, 99) [96, 99) [96, 99) 99 99 [96, 99) 99 99 [65] [96, 99) [96, 99) 99 99 [96, 99) [96, 99) [96, 99) [96, 99) [73] 99 [96, 99) [96, 99) [96, 99) 99 [96, 99) [96, 99) Levels: [96, 99) 99 Warning message: In min(xx[xx > upper]) : no non-missing arguments to min; returning Inf
For the quantiles requested with your data you get: > quantile(dat[1], prob=(1:4)/4, na.rm=TRUE) 25% 50% 75% 100% 98 98 99 99 So you only have 2 groups in your answer. Why are you persisting in trying to do this operation when you only had 4 values to start with?
David. > > On Wed, Oct 17, 2012 at 3:30 PM, arun <smartpink111 at yahoo.com> wrote: > >> Hi, >> Try this: >> dat1<-structure(list(var1 = c(97, 97, 98, 98, 97, 99, 97, >> 98, 99, 98, 99, 98, 98, 97, 97, 98, 99, 98, 96, 98, 98, 99, 98, >> 98, 99, 99, 98, 99, 98, 99, 99, 99, 99, 98, 99, 96, 99, 98, 98, >> 99, 97, 98, 99, 99, 97, 99, 99, 98, 98, 98, 99, NA, 99, 98, 98, >> 98, 98, 98, 98, 98, 99, 99, 98, 99, 99, 98, 98, 99, 99, 97, 98, >> 98, 98, 99, 98, 98, 98, 99, 98, 98)), .Names = "var1", row.names = >> c(NA, >> 80L), class = "data.frame") >> unique(dat1[,1]) >> #[1] 97 98 99 96 NA >> dat2<-dat1[!is.na(dat1)] >> cut2(dat2,g=4) >> # [1] (96,97] (96,97] (97,98] (97,98] (96,97] (98,99] >> (96,97] >> #[8] (97,98] (98,99] (97,98] (98,99] (97,98] (97,98] >> (96,97] >> #[15] (96,97] (97,98] (98,99] (97,98] (-Inf,96] (97,98] >> (97,98] >> #[22] (98,99] (97,98] (97,98] (98,99] (98,99] (97,98] >> (98,99] >> #[29] (97,98] (98,99] (98,99] (98,99] (98,99] (97,98] >> (98,99] >> #[36] (-Inf,96] (98,99] (97,98] (97,98] (98,99] (96,97] >> (97,98] >> #[43] (98,99] (98,99] (96,97] (98,99] (98,99] (97,98] >> (97,98] >> #[50] (97,98] (98,99] (98,99] (97,98] (97,98] (97,98] >> (97,98] >> #[57] (97,98] (97,98] (97,98] (98,99] (98,99] (97,98] >> (98,99] >> #[64] (98,99] (97,98] (97,98] (98,99] (98,99] (96,97] >> (97,98] >> #[71] (97,98] (97,98] (98,99] (97,98] (97,98] (97,98] >> (98,99] >> #[78] (97,98] (97,98] >> #Levels: (-Inf,96] (96,97] (97,98] (98,99] >> A.K. >> >> >> >> >> ----- Original Message ----- >> From: Charles Determan Jr <deter088 at umn.edu> >> To: r-help at r-project.org >> Cc: >> Sent: Wednesday, October 17, 2012 3:42 PM >> Subject: [R] cut2 error >> >> To R users, >> >> I am trying to use cut2 function from the 'Hmisc' library. >> However, when I >> try and run the function on the following variable, I get an error >> message >> (displayed below). I suspect it is because of the NA but I have no >> idea >> how to address the error. Many thanks to any insights. >> >> structure(list(var1 = c(97, 97, 98, 98, 97, 99, 97, >> 98, 99, 98, 99, 98, 98, 97, 97, 98, 99, 98, 96, 98, 98, 99, 98, >> 98, 99, 99, 98, 99, 98, 99, 99, 99, 99, 98, 99, 96, 99, 98, 98, >> 99, 97, 98, 99, 99, 97, 99, 99, 98, 98, 98, 99, NA, 99, 98, 98, >> 98, 98, 98, 98, 98, 99, 99, 98, 99, 99, 98, 98, 99, 99, 97, 98, >> 98, 98, 99, 98, 98, 98, 99, 98, 98)), .Names = "var1", row.names = >> c(NA, >> 80L), class = "data.frame") >> >> cut2(dat[,1], g=4) >> >> Warning message: >> In min(xx[xx > upper]) : no non-missing arguments to min; returning >> Inf >> >> Regards, >> Charles >> >> [[alternative HTML version deleted]] >> >> ______________________________________________ >> R-help at 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. >> >> > > [[alternative HTML version deleted]] > > ______________________________________________ > R-help at 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. David Winsemius, MD Alameda, CA, USA