[R] probability from different values
Matthias Kohl
Matthias.Kohl at stamats.de
Thu Dec 27 15:22:11 CET 2007
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.
>
More information about the R-help
mailing list