search for: sigmasq

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...