cut2 error
Hi Charles,
Sorry, now I remember that I was working with another cut2() function (from Rui's solution) when I tried your data.
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
cut2 <- function(x, g = 0){
??? cut(x, breaks = c(-Inf, seq(min(x), max(x), length.out = g)))
}
dat2<-dat1[!is.na(dat1)]
cut2(dat2,g=4)
A.K.
From: Charles Determan Jr <deter088 at umn.edu>
To: David Winsemius <dwinsemius at comcast.net>
Cc: arun <smartpink111 at yahoo.com>; R help <r-help at r-project.org>
Sent: Wednesday, October 17, 2012 5:04 PM
Subject: Re: [R] cut2 error
To: David Winsemius <dwinsemius at comcast.net>
Cc: arun <smartpink111 at yahoo.com>; R help <r-help at r-project.org>
Sent: Wednesday, October 17, 2012 5:04 PM
Subject: Re: [R] cut2 error
David, I am pursuing this because this is just one variable analyzed within a loop where most have the 4 defined quartiles to assign 4 groups.? When I see that error, I wanted to understand it to see if there is something wrong.? That is why I isolated it and tried the particular variable and still got the error.? Is this something to just ignore and add a component in the loop saying there isn't four groups with that variable? Regards On Wed, Oct 17, 2012 at 3:58 PM, David Winsemius <dwinsemius at comcast.net> wrote: >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 > >