search for: binomin

Dear All I have run the following GLM binominal model on a dataset composed by the following variables: TRAN_DURING_CAMP_FLG enviados bono_recibido 0 1 benchmark 0 1 benchmark 0 1 benchmark 0 1 benchmark 0 1...

Hi I have two vectors, prob and size, and I want to add the random deviates of these two, i.e. sum( sapply( 1:length(prob), function(n){ rbinom(1, size(n), prob(n) } ) ) My problem is that I have to do this for a large number of value combinations. Is there a faster way of doing this? Rainer

...; -----Oprindelig meddelelse----- >> Fra: r-help-bounces at r-project.org >> [mailto:r-help-bounces at r-project.org] P? vegne af Rainer M Krug >> Sendt: 24. oktober 2007 09:11 >> Til: Charles C. Berry >> Cc: r-help >> Emne: Re: [R] analytical solution to Sum of binominal >> distributed random numbers? >> >> Hi Charles >> >> thanks for the pointing out that size and prob can be vectors >> as well - I tried it out but used 1 as the number of >> observations, assuming that and it only gave me one randon >> mumbewr...

Hi, this is my first time using R. I want to simulate the following process: "in a population of size N, there are i individuals bearing genotype A, the number of those bearing A is j in the next generation, which following a binominal distribution (choose j from 2*N, the p is i/2*N), to plot the probability of the next generations, my script is as follows. It cannot run successfully, declaring that the "ylim should be limited. " I wonder where the bug is. Thanks very much! generation<-function(i,N) { m<-1;...

...Recently I conduct a project about expected return that involves using R for calculations. I encounter three questions as follow: Q1. How to use R to demonstrate bootstrapping and Monte Carlo simulation under 10,000 times of random selection? Q2. How to do a loop in R? Q3. How to use R to do binominal tree calculation under multiple periods? I'm perplexed by the R commands, so any hints or examples will be highly appreciated. Andy _________________________________________________________________

Hello, I like to fit data against a negative binominal distribution x2<-c(rep(10,14),rep(9,8),rep(8,13),rep(7,11),rep(6,6),rep(5,18),rep(4,7),re p(3,21),rep(2,33),rep(1,55),rep(0,225)) f2<-fitdist(x2,"nbinom",method="mle") plot(f2) summary(f2) gofstat(f2) I receive the following result: Fitting of the distrib...

...4 1 BONO6EUROS 626 1793 What I want achive is if there is a redemption or purchasing rate significatively different between each group, and see between which group there is difference Now, I have the following doubts: a. I understand I should run a multiple comparison test, like maybe a GLM with binominal distribution, but I'm not sure it is a correct procedure, considering that some of the groups (for instance the fourth one, with n=83) are quite smaller than the main other groups. Is the model I choose the correct one for this kind of analysis, and I should exclude the smaller groups? b. I...

R 2.3.0 Linux, SuSE 10.0 Hi I have two problems with mle - probably I am using it the wrong way so please let me know. I want to fit different distributions to an observed count of seeds and in the next step use AIC or BIC to identify the best distribution. But when I run the script below (which is part of my original script), I get one error message for the first call of mle: Error in