similar to: Quartile regression question

Hi, I am wondering if there is a package for doing conditional logistic regression for nested case-control study as described in "Estimation of absolute risk from nested case-control data" by Langholz and Borgan (1997) where Horvitz-Thompson sampling weight (log of (number in the risk set divided by the number sampled)) is used with regression. In SAS Proc Phreg, this is implemented

hello. I'm a student that work on R.(version 1.7.0) I need to some data frames such as "possum","Beams"'"kiwishade","dolphins",... . I typed help.search("possum") and so on,but the response was "No help files found with alias or title matching 'possum' using fuzzy matching.". I dont know what is "fuzzy matching"

I am using the describe function in prettyR. I would like to add the 25% 75% quartiles to the summary table how do I do this I have tried describe(x.f, num.desc=c("mean", "median", "sd", "min", "max", "skewness", "quantile(x.f, na.rm=T, probs=seq(0.25, 0.75))", "valid.n")) help -- Let's not spend our time

Why am I getting a wrong result for quartiles? here is my code: > cbiomass = c(910, 1058, 929, 1103, 1056, 1022, 1255, 1121, 1111, 1192, > 1074, 1415) > summary(cbiomass) > IQR(cbiomass) The result R gives me is: For the summary > Min. 1st Qu. Median Mean 3rd Qu. Max. 910 1048 1088 1104 1139 1415 For IQR > 91.25 ********* The true Q1 is 1039

Hello, I have a doubt, but it is more statistic than just about R: How the people deal usually with negative percentile and quartile? In my concrete case, I want to know the distribution of an error, so the nearest it is to 0 the better (I think it would be optimun to have the nearest values to 0 in the lower percentiles). The best result would be making the percentile of the absolute value

Greetings R users, My goal is to generate quartile groups of each variable in my data set. I would like each experiment to have its designated group added as a subsequent column. I can accomplish this individually with the following code: brks <- with(data_variables, cut2(var2, g=4)) #I don't want the actual numbers, I need a numbered group data$test1=factor(brks,

Hi, all. I have a list: process <- c( 5 , 7 , 4 , 1 , 4 , 1) and I'd like to get each half (or each third or each quartile) of the list: process.firsthalf would be (5, 7, 4) and process.secondhalf would be (1, 4, 1). note that I'm not interested in the numeric quartiles (then I could use quantile or several other functions). what is the best way to get this kind of thing? for

R users: I am trying to replicate the boxplot output I achieve with Minitab in R. I realize that R gives the user many more options on the algorithm used to calculate the IQR than Minitab, so I concentrated on type=6 when using the quantile() function in R. The problem I am having is setting the upper and lower limit of the whisker based on the nearest actual data that should be included. If

Hello, I'm currently using a simple plot to visualize some mean values. I'm having ~200 datapoints on the x-axis, each has 10 records. I'm currently plotting only the mean value of each of the datapoints. What I need is a way to visualize the quartiles/error/whatever of these points. I thought about boxplots, but I have to many points on the xaxis - it would be impossible to see

When data contains both factor and numeric variables, how to get quartiles for all numeric variables? n <- 100 x1 <- runif(n) x2 <- runif(n) x3 <- x1 + x2 + runif(n)/10 x4 <- x1 + x2 + x3 + runif(n)/10 x5 <- factor(sample(c('a','b','c'),n,replace=TRUE)) x6 <- factor(1*(x5=='a' | x5=='c')) data1 <- cbind(x1,x2,x3,x4,x5,x6) data

For boxplot(), is it possible to pass in a parameter to change the default way that the 1st and 3rd quartiles are computed? (specifically, I'd like to use type 6 described in the quantile function). Also, what are the options for how outliers are computed, and how can one change them? Thank you [[alternative HTML version deleted]]