Displaying 2 results from an estimated 2 matches for "sigmat".
Did you mean:
sigma
2007 May 27
2
[Bioc-devel] promptClass
...mptClass.R
===================================================================
--- promptClass.R (revision 41719)
+++ promptClass.R (working copy)
@@ -165,7 +165,7 @@
if (nmeths > 0) {
.meths.body <- " \\describe{"
for (i in 1:nmeths) {
- .sigmat <- sigsList(methnms[i], where)
+ .sigmat <- sigsList(methnms[i], where=whereClass)
for (j in seq_along(.sigmat)) {
if (!all(is.na(match(.sigmat[[j]],clName)))) {
methn.i <- escape(methnms[i])
Index: RMethodUtils.R
==================...
2008 May 16
0
How to determine sensible values for 'fnscale' and 'parscale' in optim
...se when finding,
for example, the MLEs of a bivariate normal distribution using optim.
Here's code for this example:
-----------------------------------------------
library(MASS) # needed mvrnorm
library(mvtnorm) # needed for dmvnorm
set.sed(20080516)
n=1000
mu1=3
mu2=5
sig1=7
sig2=20
rho=.5
sigmat=matrix(c(sig1^2,sig1*sig2*rho,sig1*sig2*rho,sig2^2),2)
xy=mvrnorm(n,c(mu1,mu2),sigmat) # n = 1000 observations from this
# distribution.
obj=function(par,xy) # The function to maximize.
{
mu=par[1:2]
sigmat=matrix(c(par[3]^2,par[3]*par[4]*par[5],par[3]*par[4]*par[5]...