Displaying 3 results from an estimated 3 matches for "sseed".
Did you mean:
speed
2006 Aug 30
2
density() with from, to or cut and comparrison of density()
Hi
the function density() does normally integrate to one - I've checked it
and it works and I also read the previous threads.
But I realised that it does not integrate to one if I use from, to or cut.
My scenario: I simulated densities of a plants originating from an sseed
source at distance zero. Therefore the density of the plants will be
highest close to zero. Is there anything I can do to have this pattern?
If I use 'from' or 'cut', the resulting densities do not integrate to
one which I need as I want to compare different density curves.
Ny seco...
1997 Nov 24
0
R-alpha: random number generator -- S-plus's
...4 5 6 7 8 9 10 12 12 ]
###--> only [4:6] really change !!
###- They have 6 + 6 + 2 = 14 bits each (the ranges are 0:63, 0:63, and 0:3),
###- i.e., 2^14 = 16384 different possibilities.
###- ~~~~ ~~~~~
###===> Period of set.seed(.) =: P | 2^14
X
## test this period
tst.sseed <- function(i)
{
X ## Purpose: test set.seed(.)
X ## -------------------------------------------------------------------------
X ## Author: Martin Maechler, Date: 9 May 94, 11:54
X set.seed(i)
X .Random.seed
}
X
for (i in c(0:7, 11, 17, 20, 23, 40, 53, 59, 67, 76, 97)) {
X cat("i=&...
2009 Feb 12
3
proposed simulate.glm method
...eed <- attr(S1, "seed")
identical(S1, simulate(mod1, nsim = 4))
S2 <- simulate(mod1, nsim = 200, seed = 101)
rowMeans(S2)/10 # after correcting for binomial sample size, should be
about
fitted(mod1)
plot(rowMeans(S2)/10)
lines(fitted(mod1))
## repeat identically:
(sseed <- attr(S2, "seed")) # seed; RNGkind as attribute
stopifnot(identical(S2, simulate(mod1, nsim = 200, seed = sseed)))
}
--
Ben Bolker
Associate professor, Biology Dep't, Univ. of Florida
bolker at ufl.edu / www.zoology.ufl.edu/bolker
GPG key: www.zoology.ufl.edu/bolker/benbolke...