[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