Displaying 9 results from an estimated 9 matches for "sigmasq".

Did you mean:
sigmask

2005 Oct 21

1

finite mixture model (2-component gaussian): plotting component gaussian components?

...####
data(faithful)
##################
# (2) get model
##################
library(mclust)
MyMixtureModel<-summary(EMclust(faithful$eruptions),faithful$eruptions)
##################
# (3) plot mixture model
##################
attach(MyMixtureModel)
mclust1Dplot(data=faithful$eruptions,z=z,mu=mu,sigmasq=sigmasq,pro=pro,ask=FALSE,type=c("density"))
do.call("mclust1Dplot",c(list(data=faithful$eruptions,ask=FALSE,type=c("density")),MyMixtureModel))
##################
# (4) plot components
##################
???
Any information you might be able to shed on this would be...

2010 Sep 22

0

Help with mclust package

...lust package on a variable "Tuberculin indurations"
recorded as mm. The file has only one variable. When I run the package I get
NULL value for mu and sigma. Can anybody say why?
This is the program:
library("mclust")
mc<-Mclust(x.trab,G=1:9,warn=TRUE)
mc
mc$mu
sqrt(mc$sigmasq)
and the output I get is
> library("mclust")
> mc<-Mclust(x.trab,G=1:9,warn=TRUE)
There were 11 warnings (use warnings() to see them)
> mc
best model: unequal variance with 3 components
> mc$mu
NULL
> sqrt(mc$sigmasq)
Error in sqrt(mc$sigmasq) : Non-numeric argument...

2007 Oct 10

11

please help me

dear list
I am student M.S. statistics in department statistics . I am working in the function "nls" in the [R 2.3.1] with 246 data and want to fit the "exp" model to vectors( v and u ) but I have
a problem to use it
u
5.000000e-13 2.179057e+03 6.537171e+03 1.089529e+04 1.525340e+04
1.961151e+04 2.396963e+04 2.832774e+04 3.268586e+04 3.704397e+04
4.140209e+04

2005 Dec 21

0

Help with Krige.conv using linear models

...ond
that. I have been having problems using krige.conv() to get a decent
kriged map using the linear model. The code I am using from my data is
as follows:
>modeltest=variofit(variotest, weights=?cressie?, cov.model=?linear?,
ini.cov.pars=c(80,1))
The output parameters are tausq = 9.855, sigmasq = 0.0087, phi=1.0
>krig=krige.conv(data, krige=krige.control(type.krig=?ok?,
obj.model=modeltest), locations=pred.grid)
At this point, krig$predict values have little to no variability (1.897
+/- 0.004), where I would expect values between 0 ? 15.
By running the same data, except using...

2005 Oct 22

0

package mclust: cdens, EMclust?

...y
real dataset will not have peaks this well separated,
but I needed to find a small example.)
##################
data(faithful)
library(mclust)
MyMixtureModel<-summary(EMclust(faithful$eruptions),faithful$eruptions)
attach(MyMixtureModel)
mclust1Dplot(data=faithful$eruptions,z=z,mu=mu,sigmasq=sigmasq,pro=pro,ask=FALSE,type=c("density"))
do.call("mclust1Dplot",c(list(data=faithful$eruptions,ask=FALSE,type=c("density")),MyMixtureModel))
# plot components
???
##################
Any information you might be able to shed on this
would be very much apprec...

2012 Oct 04

1

geoRglm with factor variable as covariable

...n =
100) #trial error
outmcmc.5 = glsm.mcmc(geoData9093, model= model.5, mcmc.input = mcmc.5)
mcmcobj.5 = prepare.likfit.glsm(outmcmc.5)
lik.5 = likfit.glsm(mcmcobj.5, ini.phi = 0.3, fix.nugget.rel = F)/
And the summary of lik.5 is:
likfit.glsm: estimated model parameters:
beta sigmasq phi tausq.rel
"1.2781" "0.5193" "0.0977" "0.0069"
likfit.glsm : maximised log-likelihood = 43.62
I'm fairly new to geostatistics, but I thought using a factor variable as
covariable would give me 4 intercepts (beta) as I have 4 levels in my...

2006 Jul 25

1

HELP with NLME

...log(1+exp(eta1));
eta2 = a0 + a1*lnndsTs;
llnorm = -1/(2*sigsq)*(lnbldT - eta2)**2 - .5*log(sigsq);
ll = llbin + llnorm;
end;
else do;
eta2 = a0 + a1*lnndsTs;
eta1 = b0 + b1*eta2 + u;
llbin = anybc.cens.ind*eta1 - log(1+exp(eta1));
ll = llbin;
end;
sigma2 = sigmasq*b1**2 /*variance of random effect;
model anybc.cens.indic ~ general(ll);
random u ~ normal(0,sigma2) subject = CaseID; */;
run;
####################################################################
R Code with error message:
me.km.nlme <- nlme(model = anybc.cens.indic ~
vs.flag*((anybc....

2007 Mar 23

0

plotting dnorm() issued from mclust models

...ve a problem in fitting lines() of the normal distributions
identified with Mclust on a histogram or a mclust1Dplot. Here is some
sample code to explain :
set.seed(22)
foo <- c(rnorm(400, 10, 2), rnorm(500, 17, 4))
mcl <- Mclust(foo, G=2)
mcl.sd <- sqrt(mcl$parameters$variance$sigmasq)
mcl.size <- c(length(mcl$classification[mcl$classification==2]),
length(mcl$classification[mcl$classification==1]))
x <- pretty(c(0:44), 100)
#### my plot of histogram and lines of normal distributions
#### SEEMS OK (or am I wrong ?) using frequencies :
histA <- hist(foo, bre...

2010 Mar 14

3

the error in DoSimulateRF function

.../geoR
geoR version 1.6-27 (built on 2009-10-15) is now loaded
-------------------------------------------------------------
grf: generating grid 300 * 300 with 90000 points
grf: process with 1 covariance structure(s)
grf: nugget effect is: tausq= 0
grf: covariance model 1 is: spherical(sigmasq=1, phi=0.1)
grf: simulation using the function GaussRF from package RandomFields
Error in DoSimulateRF(n = n, reg = register, paired = paired) :
.Random.seed is not an integer vector but of type 'double'
I tried to trace back the source of the error but installations seem
correct to me.
r...