There are still a lot of things that you did not specify, but maybe this
will get you started:
n.groups <- sample(3:6, 1)
n.species <- sample(2:10, n.groups, replace=TRUE)
n.animals <- sample( 10:25, sum(n.species), replace=TRUE )
mu.g <- rnorm(n.groups, 100, 50)
mu.species <- rnorm( sum(n.species
On Fri, Nov 9, 2012 at 3:12 PM, Andrea Goijman
wrote:
> Thanks Sarah,
>
> Ok, I'll start over... forget the functions that I wrote before and lets
> start from scratch.
> Suppose I want to simulate data; for example n species that belong to g
> groups. Lets say 10 species, where 4 belong to group
Thanks Sarah,
Ok, I'll start over... forget the functions that I wrote before and lets
start from scratch.
Suppose I want to simulate data; for example n species that belong to g
groups. Lets say 10 species, where 4 belong to group 1, 3 to group 2, and
the rest to group 3.
I want to simulate those
Andrea,
I don't understand what you want. Why don't you start over, and
explain what the inputs are, and what your desired outputs are.
Using your previous code:
# pre-specified number of groups
groups <- 3
# what is species? It needs to be defined before
# we can assign a value to a component o
Yes, I tried your suggestion, but it didn't help. It just creates a
tri-dimentional array for p... and that is not what I want p per species
(within groups) to be an unbalanced array...
I need to have a single value of p per species, and the total number of
elements in p would be whatever number r
Andrea,
I simply meant that I couldn't run your code assigning a value to
species$groups
because the code didn't include any information about creating the R
object species. Thus, I changed the name of that R object to
speciesgroups and altered your code so that it runs by removing the
erroneous c
Hi Sarah,
Maybe I expressed myself wrong, but so far, I don't have the species; I'm
just wanting to generate simulated data. For example, creating an
unbalanced (and random) number of species per group, and then run the
for-loops
for example:
Group1: 3 species
Group2: 5 species
Group3: 8 species
Hi Andrea,
Without knowing what species is, I can't run your code as is, but try this:
groups<-3
speciesgroups <- as.integer(runif(groups,1,10))
# I'd use sample(1:10, groups, replace=TRUE)
p<-array(NA,dim=speciesgroups)
Mostly you're trying to use c() on something that's already a vector.
Sa
I know this seems like a very easy question (and maye it is) but I've been
trying to simulate nested data and been unsucessful so far..
I want to simulate a varying number of species within a group; and then
create an array to store the results of my for-loop. For example:
groups<-3
species$group
9 matches
Mail list logo
