Displaying 20 results from an estimated 420 matches for "quantile".

2006 Mar 16

5

problem for wtd.quantile()

Dear R-users,
I don''t know if there is a problem in wtd.quantile (from library "Hmisc"):
--------------------------------
x <- c(1,2,3,4,5)
w <- c(0.5,0.4,0.3,0.2,0.1)
wtd.quantile(x,weights=w)
-------------------------------
The output is:
0% 25% 50% 75% 100%
3.00 3.25 3.50 3.75 4.00
The version of R I am using is: 2.1.0
Best,Jing

2012 Jul 14

1

Quantile Regression - Testing for Non-causalities in quantiles

Dear all,
I am searching for a way to compute a test comparable to Chuang et al.
("Causality in Quantiles and Dynamic Stock
Return-Volume Relations"). The aim of this test is to check wheter the
coefficient of a quantile regression granger-causes Y in a quantile range. I
have nearly computed everything but I am searching for an estimator of the
density of the distribution at several points of th...

2008 Sep 27

3

quantile / centile

I''m wondering if there is a simple way to assign a quantile to a vector in a
data frame, much like one could in Stata using centile. Let''s say I want 100
slices in my assignation. I can easily see what the limits of each slice by
using quantile:
quantile(my.df$my.var, probs=seq(0, 1, 0.01))
But how do I assign the appropriate value to each row/...

2009 Jul 21

5

how to change the "quantile" method in bwplot

Hi, everyone,
Since quantile calculation has nine different methods in R, I wonder how I
specify a method when calling the bwplot() in lattice. I couldn''t find any
information in the documentation. Thanks.
--
Jun Shen PhD
PK/PD Scientist
BioPharma Services
Millipore Corporation
15 Research Park Dr.
St Charles, M...

2018 Apr 19

2

create multiple categorical variables in a data frame using a loop

Hi All,
I want to create a categorical variable, cat.pfoa, in the file of pfas.pheno (a data frame) based on log2pfoa values. I can do it using the following code.
pfas.pheno <-within(pfas.pheno, {cat.pfoa<-NA
cat.pfoa[pfas.pheno$log2pfoa <=quantile(pfas.pheno$log2pfoa,0.25, na.rm =T)]<-0
cat.pfoa[pfas.pheno$log2pfoa >=quantile(pfas.pheno$log2pfoa,0.75, na.rm =T)]<-2
cat.pfoa[pfas.pheno$log2pfoa >=quantile(pfas.pheno$log2pfoa,0.25, na.rm =T)
&pfas.pheno$log2pfoa <=quantile(pfas.pheno$log2pfoa,0.75, na.rm =T)]&...

2017 Jun 15

1

"reverse" quantile function

David,
thanks for the response. In your response the quantile function (if I see correctly) runs on the columns versus I need to run it on the rows, which is an easy fix, but that is not exactly what I had in mind... essentially we can remove t() from my original code to make "res" look like this:
res<-apply(z, 1, quantile, probs=c(0.3))
but...

2017 Jun 16

2

"reverse" quantile function

It would depend on which one of the 9 quantile definitions you are using. The discontinuous ones aren't invertible, and the continuous ones won't be either, if there are ties in the data.
This said, it should just be a matter of setting up the inverse of a piecewise linear function. To set ideas, try
x <- rnorm(5)
curve(quantile(...

2016 Apr 19

2

Interquartile Range

...#39;s
function had a typo in it - was that the problem? Or did you not
change the call to ddply to use that function. Here is something
that might "work" for you:
library(plyr)
data <- data.frame(groupColumn=rep(1:5,1:5), col1=2^(0:14))
myIqr <- function(x) {
paste(round(quantile(x,0.25),0),round(quantile(x,0.75),0),sep="-")
}
ddply(data, ~groupColumn, summarise, col1_myIqr=myIqr(col1),
col1_IQR=stats::IQR(col1))
# groupColumn col1_myIqr col1_IQR
#1 1 1-1 0
#2 2 2-4 1
#3 3 12-24 12
#4...

2017 Jun 15

1

"reverse" quantile function

...as via R-help <r-help at r-project.org> wrote:
>
> Dear All,
>
> we have:
>
> t<-seq(0,24,1)
> a<-10*exp(-0.05*t)
> b<-10*exp(-0.07*t)
> c<-10*exp(-0.1*t)
> d<-10*exp(-0.03*t)
> z<-data.frame(a,b,c,d)
>
> res<-t(apply(z, 1, quantile, probs=c(0.3)))
>
>
>
> my goal is to do a 'reverse" of the function here that produces "res" on a data frame, ie: to get the answer 0.3 back for the percentile location when I have "res" available to me... For a single vector this would be done using ec...

2012 Oct 30

12

standard error for quantile

Dear all
I have a question about quantiles standard error, partly practical
partly theoretical. I know that
x<-rlnorm(100000, log(200), log(2))
quantile(x, c(.10,.5,.99))
computes quantiles but I would like to know if there is any function to
find standard error (or any dispersion measure) of these estimated
values.
And here is a the...

2016 Apr 19

1

Interquartile Range

To be precise:
paste(round(quantile(x,0.25),0),round(quantile(x,0.75),0),sep="-")
is an expression that evaluates to a character string:
"round(quantile(x,.25),0) - round(quantile(x,0.75),0)"
no matter what the argument of your function, x. Hence
return(paste(...)) will return this exact character string and ne...

2016 Apr 19

1

Interquartile Range

...ith having an open mind is that people keep coming along
and sticking things into it."
-- Opus (aka Berkeley Breathed in his "Bloom County" comic strip )
On Tue, Apr 19, 2016 at 10:25 AM, Bert Gunter <bgunter.4567 at gmail.com> wrote:
> To be precise:
>
> paste(round(quantile(x,0.25),0),round(quantile(x,0.75),0),sep="-")
>
> is an expression that evaluates to a character string:
> "round(quantile(x,.25),0) - round(quantile(x,0.75),0)"
>
> no matter what the argument of your function, x. Hence
>
> return(paste(...)) will return th...

2009 Jun 09

5

Non-linear regression/Quantile regression

Hi,
I''m relatively new to R and need to do a quantile regression. Linear
quantile regression works, but for my data I need some quadratic function.
So I guess, I have to use a nonlinear quantile regression. I tried the
example on the help page for nlrq with my data and it worked. But the
example there was with a SSlogis model. Trying to write
dat.nl...

1998 Feb 12

1

R-beta: Quantile function

Is the following behaviour of the quantile function what one would expect?
> a <- 1:100
> quantile(a,.6)
60%
60.4
Philippe
-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-
r-help mailing list -- Read http://www.ci.tuwien.ac.at/~hornik/R/R-FAQ.html
Send "info", "...

2011 Mar 28

2

quantile function -> I need only the quantile value itself

Hi,
I am using the quantile function currently and I have just bumped into a little problem.
I have a very small data frame something like this:
small_df <- c(7,3,4,7,1,10,12,1,12,4,4,8,6,11,9,10,4,13,3,9,6,5,2,10,7,14,2,7,10,10,7,8,2,11,3,10,11,3,11,14,12,7,6,11)
small_df
Now in the next step I would like to calculate...

2005 Jul 09

2

Quantile normalization and NA

Hi,
I am new to R,
I am doing quantile normalization with a dat matix of
384X124 and I find that while computing the quantile
normailzation it introduces ''NA'' into some of the
cells, can someone help me to overcome this problem ?
This is the command that goes like upto g62 for 124
colomns
>g1 <- normalize.quan...

2018 Apr 19

1

create multiple categorical variables in a data frame using a loop

Hello,
When programming it is better to use dat[["variable"]] than dat$variable.
So your code could be
pfas.pheno[[cat.var]] <- NA
pfas.pheno[[cat.var]][pfas.pheno[,i] <= quantile(pfas.pheno[,i],0.25,
na.rm =T)] <- 0
etc.
Untested.
Hope this helps,
Rui Barradas
On 4/19/2018 7:20 PM, Ding, Yuan Chun wrote:
> Hi All,
>
> I want to create a categorical variable, cat.pfoa, in the file of pfas.pheno (a data frame) based on log2pfoa values. I can do it u...

2016 Apr 19

1

Interquartile Range

...ple keep coming along
> and sticking things into it."
> -- Opus (aka Berkeley Breathed in his "Bloom County" comic strip )
>
>
> On Tue, Apr 19, 2016 at 10:25 AM, Bert Gunter <bgunter.4567 at gmail.com> wrote:
>> To be precise:
>>
>> paste(round(quantile(x,0.25),0),round(quantile(x,0.75),0),sep="-")
>>
>> is an expression that evaluates to a character string:
>> "round(quantile(x,.25),0) - round(quantile(x,0.75),0)"
>>
>> no matter what the argument of your function, x. Hence
>>
>> retur...

2002 May 14

2

quantile() and boxplot.stats()

Hello,
I faced something I can''t understand. When I use boxplot.stats(1:10) and
quantiles(1:10) the results are different for 25% and 75%:
> boxplot.stats(1:10)
$stats
[1] 1.0 3.0 5.5 8.0 10.0
> quantile(1:10)
0% 25% 50% 75% 100%
1.00 3.25 5.50 7.75 10.00
Actually, I expected the value 3 for 25% and 8 for 75% as results of
quantile(1:10). Can you please expla...

2016 Apr 19

1

Interquartile Range

...gt; -- Opus (aka Berkeley Breathed in his "Bloom County" comic strip )
>> >
>> >
>> > On Tue, Apr 19, 2016 at 10:25 AM, Bert Gunter <bgunter.4567 at gmail.com>
>> wrote:
>> >> To be precise:
>> >>
>> >> paste(round(quantile(x,0.25),0),round(quantile(x,0.75),0),sep="-")
>> >>
>> >> is an expression that evaluates to a character string:
>> >> "round(quantile(x,.25),0) - round(quantile(x,0.75),0)"
>> >>
>> >> no matter what the argument of yo...