On May 16, 2010, at 6:16 PM, Agustín Muñoz M. (AMFOR) wrote:
Dear, I want to make an application to calculate quantile within a
For()
I tried the following without success:
ej.
date
p_val <- matrix(sample(10, 1000, replace=TRUE), 200,5)
test 1
rr <- paste("p_val$",names(p_val[1]), sep="")
p_val does not have any names>
names(matrix(1:4, ncol=2))
NULL
quant <- quantile(rr, probs = c(0, 10, 20, 30, 40, 50, 60, 70, 80, 90,
100)/100, na.rm=FALSE, type=1)
test 2
rr <- noquote(paste("p_val$",names(p_val[1]), sep=""))
quant <- quantile(rr, probs = c(0, 10, 20, 30, 40, 50, 60, 70, 80, 90,
100)/100, na.rm=FALSE, type=1)
test 3
quant <- quantile(p_val[1], probs = c(0, 10, 20, 30, 40, 50, 60, 70,
80,
90, 100)/100, na.rm=FALSE, type=1)
Your argument is a single number:
> matrix(1:4, ncol=2)[1]
[1] 1
The only thing that works for me is:
quant <- quantile(p_val$1, probs = c(0, 10, 20, 30, 40, 50, 60, 70,
80,
90, 100)/100, na.rm=FALSE, type=1)
but the evil, he will not let me do this for a dynamic For()
I think more details are needed. You have provide neither full code,
not an error message for you failed attempts. You probably want to
work with:
p_val[, i]
and should test it with
p_val <- matrix(sample(10, 1000, replace=TRUE), 200,5)
> quant <- quantile(p_val[,1] , probs = c(0, 10, 20, 30, 40, 50, 60,
70, 80,
+ 90, 100)/100, na.rm=FALSE, type=1)
>
> quant
0% 10% 20% 30% 40% 50% 60% 70% 80% 90% 100%
1 1 2 3 4 6 7 8 9 10 10
The moral: Avoid using the "$" operator for looping tasks. Use the
indexing "[" operator using the [,i] or [,1] forms instead.
I hope I can help, as they always do.
from now, thank you very much.
Atte.
Agustín
[[alternative HTML version deleted]]
______________________________________________
R-help@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
West Hartford, CT
______________________________________________
R-help@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.
