probability from different values
Matthias Kohl wrote:
one addition ...
Hello Ricardo,
another solution could be using package distr:
library(distr)
A <- c(18,18,18,19,20,21,22,23,24,25,26,27,28)
DA <- DiscreteDistribution(A)
# maybe
# support(DA)
# plot (DA)
B <- c(82,83,84,85,85,86,87,88,89,90,91,91,92)
DB <- DiscreteDistribution(B)
# support(DB)
# plot(DB)
DC <- DB - DA# convolution with of DB with (-DA)
# support(DC)
# plot(DC)
to get the probabilities for the C-values: p(DC)(support(DC))
sorry, I meant d(DC)(support(DC))
hth,
Matthias
Richard.Cotton at hsl.gov.uk wrote:
i have two vectors (A and B) and i need create another vector (C)
from the subtraction of A's values with the B's values. How can i
estimate the probability of C's values if i have differents values
combinations of A and B that can result in the same value? like 90
-20 = 70 and 91 - 21 = 70.
example:
B = {18 18 18 19 20 21 22 23 24
25 26 27 28....}
A = {82 83 84 85 85 86 87 88 89
90 91 91 92 ....}
Try this:
A = c(18,18,18,19,20,21,22,23,24,25,26,27,28)
B = c(82,83,84,85,85,86,87,88,89,90,91,91,92)
lenA = length(A)
lenB = length(B)
AA = rep(A, lenB)
BB = rep(B, each=lenA)
table(AA-BB)/(lenA*lenB)
Regards,
Richie.
Mathematical Sciences Unit
HSL
------------------------------------------------------------------------
ATTENTION:
This message contains privileged and confidential inform...{{dropped:20}}
______________________________________________ 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.
______________________________________________ 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.
______________________________________________ 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.