Re: [R] Random assignment

2010-10-15 Thread Michael Bedward
Hello again John, I was going to suggest that you just use qbinom to generate the expected number of extinctions. For example, for the family with 80 spp the central 95% expectation is: qbinom(c(0.025, 0.975), 80, 0.0748) which gives 2 - 11 spp. If you wanted to do look across a large number of

Re: [R] Random assignment

2010-10-15 Thread Michael Bedward
Hi John, I haven't read that particular paper but in answer to your question... > So if i do this for all the families it will be the same as doing the > simulation experiment > outline in the method above? Yes :) Michael On 15 October 2010 23:18, John Haart wrote: > Hi Michael, > > Thanks

Re: [R] Random assignment

2010-10-15 Thread John Haart
Hi Michael, Thanks for this - the reason i am following this approach is that it appeared in a paper i was reading, and i thought it was a interesting angle to take The paper is Vamosi & Wilson, 2008. Nonrandom extinction leads to elevated loss of angiosperm evolutionary history. Ecology Let

Re: [R] Random assignment

2010-10-15 Thread Michael Bedward
Hi John, The word "species" attracted my attention :) Like Dennis, I'm not sure I understand your idea properly. In particular, I don't see what you need the simulation for. If family F has Fn species, your random expectation is that p * Fn of them will be at risk (p = 0.0748). The variance on t

Re: [R] Random assignment

2010-10-15 Thread John Haart
Hi Denis and list Thanks for this , and sorry for not providing enough information First let me put the study into a bit more context : - I know the number of species at risk in each family, what i am asking is "Is risk random according to family or do certain families have a disproportionate

Re: [R] Random assignment

2010-10-15 Thread Dennis Murphy
Hi: I don't believe you've provided quite enough information just yet... On Fri, Oct 15, 2010 at 2:22 AM, John Haart wrote: > Dear List, > > I am doing some simulation in R and need basic help! > > I have a list of animal families for which i know the number of species in > each family. > > I a

[R] Random assignment

2010-10-15 Thread John Haart
Dear List, I am doing some simulation in R and need basic help! I have a list of animal families for which i know the number of species in each family. I am working under the assumption that a species has a 7.48% chance of being at risk. I want to simulate the number of species expected to