Displaying 16 results from an estimated 16 matches similar to: "Mclust problem with mclust1Dplot: Error in to - from : non-numeric argument to binary operator"
2005 Oct 21
1
finite mixture model (2-component gaussian): plotting component gaussian components?
Dear Knowledgeable R Community Members,
Please excuse my ignorance, I apologize in advance if this is an easy question, but I am a bit stumped and could use a little guidance.
I have a finite mixture modeling problem -- for example, a 2-component gaussian mixture -- where the components have a large overlap, and
I am trying to use the "mclust" package to solve this problem.
I need
2005 Oct 22
0
package mclust: cdens, EMclust?
Dear Knowledgeable R Community Members,
Please excuse my ignorance -- I apologize in advance
if this is an easy question, but I am a bit stumped
and could use a little guidance regarding parameters
in 2 functions in the "mclust" package.
--------------------
PROBLEM DESCRIPTION
--------------------
I have a finite mixture modeling problem -- for
example, a 2-component gaussian
2007 Mar 23
0
plotting dnorm() issued from mclust models
Dear all
I have 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]),
2011 Sep 03
2
mclust: modelNames("E") vs modelNames=("V")
Hi,
I'm trying to use the library mclust for gaussian mixture on a numeric
vector. The function Mclust(data,G=3) is working fine but the fitting is not
optimal and is using modelNames="E". When I'm trying
Mclust(data,G=3,modelName="V") I have the following message:
Error in if (Sumry$G > 1) ans[c(orderedNames, "z")] else ans[orderedNames] :
argument is
2007 Jan 11
3
batch job GLM calculations
Hello
I want to batch job the calculation of many GLM-models, extract some values and store them in a file. Almost everything in the script below works (read file, extract values and write them to file) except I fail in indexing the GLM with the modelstructure it should run. Running GLM's conventionally is no problem.
Conventionally a GLM is calculated as:
2003 Sep 26
0
Gaussian mixture models, modelNames: EEV, VEV mclust2002
Hello,
I am using the parameterized EM- Algorithm in clustering context.
I use generated datasets, containing 5000 points in 20 dimensions with 5
clusters.
The following poses a problem for me:
emEst<-me(modelName="EEV",data=data[,1:numDimensions],z=unmap(x))
emEst<-me(modelName="VEV",data=data[,1:numDimensions],z=unmap(x))
coordProj(main="Model
2007 Jul 18
2
EM unsupervised clustering
Hi All,
I have a n x m matrix. The n rows are individuals, the m columns are variables.
The matrix is in itself a collection of 1s (if a variable is observed for an
individual), and 0s (is there is no observation).
Something like:
[,1] [,2] [,3] [,4] [,5] [,6]
[1,] 1 0 1 1 0 0
[2,] 1 0 1 1 0 0
[3,] 1 0 1 1 0 0
[4,] 0 1 0
2008 Mar 26
0
out of colors in Mclust with 12 clusters
I'm running Mclust 3.0 in R-2.7.0 and have run into a situation where
the optimal number of clusters, 12, seems to be greater than the number
of colors available to Mclust. The code below, which demonstrates
errors, does not if the max no of clusters is set at 10. sessionInfo,
after a restart of R and loading of packages, follows the output.
Is this a known problem with Mclust? Can I
2011 Sep 04
2
mclust: modelName="E" vs modelName="V"
Hi,
I'm trying to use the library mclust for gaussian mixture on a numeric
vector. The function Mclust(data,G=3) is working fine but the fitting is not
optimal and is using modelNames="E". When I'm trying
Mclust(data,G=3,modelName="V") I have the following message:
Error in if (Sumry$G > 1) ans[c(orderedNames, "z")] else ans[orderedNames] :
argument is
2012 Jul 19
1
Change log(J) to log(J+1) to stop log(0) from occurring in harModel
I think the code is part of the RTAQ package but is not included in it, as I
obtained it from
https://r-forge.r-project.org/scm/viewvc.php/pkg/RTAQ/R/HAR_model.R?view=markup&root=blotter&sortby=author&pathrev=1028.
It is not my code and I make no claim to other's good work, and apologize if
I should even be posting it I am not sure, but in the transform function it
allows to
2010 Jan 06
1
positive log likelihood and BIC values from mCLUST analysis
My question is with respect to mCLUST and the values of BIC and log
likelihood. The relevant part of my R script is:
######################### BEGIN MDS ANALYSIS #########################
#load data
data <- read.table("Ecoli33_Barry.dis", header = TRUE, row.names = 1)
#perform MDS Scaling
mds <- metaMDS(data, k = Dimensions, trymax = 20, autotransform =TRUE,
noshare = 0.1,
2012 Jul 19
1
Switching log(J) to log(J+1) to avoid log(0) in HAR-RVJ model
I am working with xts dependent data, and my code is as follows (the problem
is explained throughout):
dat <- getdat("prices")
dat <- read.zoo(dat, sep = "",format="%d/%m/%Y %H:%M",
tz="", FUN=NULL, regular=TRUE,
header=TRUE, index.column=1, colClasses=c("character",
"numeric"))
dat <- as.xts(dat)
2010 Apr 19
1
What is mclust up to? Different clusters found if x and y interchanged
Hello All...
I gave a task to my students that involved using mclust to look for clusters
in some bivariate data of isotopes vs various mining locations. They
discovered something I didn?t expect; the data (called tur) is appended
below.
p <- qplot(x = dD, y = dCu65, data = tur, color = mine)
print(p) # simple bivariate plot of the data; looks fine
mod1 <- Mclust(tur[,2:3])
mod1$G
mod2
2006 May 30
0
data fixtures for non-activerecord models
I have a model that as an instance maintains a session on a website, and has
methods to return some scraped data. Instead of hitting a URL over and over
during testing, I would like to just save the source to a file and use that
like a fixture. A mock object doesn''t seem appropriate, because there is a
lot of data here. What I''m doing now is storing each separate
2012 Nov 23
2
[LLVMdev] [cfe-dev] costing optimisations
On 23.11.2012, at 15:12, john skaller <skaller at users.sourceforge.net> wrote:
>
> On 23/11/2012, at 5:46 PM, Sean Silva wrote:
>
>> Adding LLVMdev, since this is intimately related to the optimization passes.
>>
>>> I think this is roughly because some function level optimisations are
>>> worse than O(N) in the number of instructions.
>>
2010 Sep 24
7
help
Estimados
Escribo para consultar sobre el uso de modelos mixtos anidados. Los
datos que estoy analizando provienen de censos de malezas en cuatro
tipos de paisajes de la región pampeana, en los que seleccioné al azar
igual número de lotes agrícolas cultivados con tres cultivos (maíz,
soja y trigo-soja). En cada lote censé el número de especies de
malezas en tres posiciones: el alambrado, el borde