search for: modelnam

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 of length zero In addition: Warning message: In pickBIC(object[as.character(G), modelN...

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 of length zero In addition: Warning message: In pickBIC(object[as.character(G), model...

...The following example causes a crash of R: # univariate example require(mclust) X <- c(rnorm(100, mean=1), rnorm(100, mean=5)) # I also tried: X <- as.matrix(c(rnorm(100, mean=1), rnorm(100, mean=5))) Xmap <- cbind(c(rep(1, 100), rep(0, 100)), c(rep(0, 100), rep(1, 100))) Xm <- mstep(modelName="V", data=X, z=Xmap) # CRASH: em(modelName=Xm$modelName, data=X, mu=Xm$mu, sigmasq=Xm$sigmasq, pro=Xm$pro) For bivariate data the following example works: # bivariate example require(mclust) Y <- as.matrix(cbind(c(rnorm(100, mean=1), rnorm(100, mean=5)), c(rnorm(100, mean=1), rnorm...

...etty easy to answer. Thanks in advance for any help. Ken What are the pros/cons of starting a finite mixture model at the “m” step versus the “e” step (where “m” is the maximization step and “e” is the expectation step of the EM algorithm)? In particular, are there any reasons for using em(modelName=XXX) versus me(modelName=XXX). Other than MCLUST, I’ve not seen a finite mixture model “program” give such an option. Would it make sense to fit both models and take the one with the largest log likelihood? Rather than the hc() function performing cluster analysis for all of G possible clus...

It looks important and I''m too scared to remove it. The boys on IRC are saying it doesn''t do anything. Confirm/deny.

Dear list members, I am using Mclust in order to deconvolute a distribution that I believe is a sum of two gaussians. First I can make a model: > my.data.model = Mclust(my.data, modelNames=c("E"), warn=T, G=1:3) But then, when I try to plot the result, I get the following error: > mclust1Dplot(my.data.model, parameters = my.data.model$parameters, what = "density") Error in to - from : non-numeric argument to binary operator Also, I'd like to allow for...

...a coxph object using the Breslow estimator). Lets say I'm interested in a survival model with two (dichotomous) factors and a continuous covariate. Variable Possible Values Factor1 0 or 1 Factor2 0 or 1 Covariate 0 to 100 I fit my model: modelname <- coxph(Surv ~ Factor1 + Factor2 + Covariate, data = data) If I then ask for: baselineA <- basehaz(modelname, centered=FALSE) I am fairly certain that baselineA will provide me with the cumulative hazard evaluated at Factor1 = 0, Factor2 = 0, Covariate = 0. Yet, if I ask for: baselineB &l...

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 EEV",data,dimens=c(1,2),type="classification",ask=F,mu=emEst$mu,sigma=em Est$sigma,z=emEst$z) coordP...

I keep on receiving the message below after submitting the following line using the mclust package. m2 is a 99 X 1 column vector. * em(modelName = "E", m2, mu = c(25, 50), sigmasq=10, pro = c(0.4, 0.6)) Error in as.double.default(data) : (list) object cannot be coerced to double. Why do I receive this error? Thank, Brian C. Newquist Research Statistician Constella Health Sciences 2605 Meridian Parkway, Suite...

Hi.... Is there some way to configure Rails to not create a db/migrate/###_model_name.rb when I do a script/generate model ModelName? It''s a little annoying when I''m setting up a new application since I prefer creating all the first tables in a 001_inital_schema.rb file. Then when I generate models for those tables I get migration scripts for all of them again that I have to delete. And it''s a bit an...

How do I get scaffold to generate code for a table whos name is not plural of the modelname I give to it. I.E., If I have an existing database, and wish to generate scaffold code for existing tables, is this possible??? _______________________________________________ Rails mailing list Rails-1W37MKcQCpIf0INCOvqR/iCwEArCW2h5@public.gmane.org http://lists.rubyonrails.org/mailman/listinf...

...SIZE/Kandidatenmodelle_Berechnung/inputfiles/" t.tuti <- read.table(paste(t.url, "All_animals.txt", sep=""),header=T) collect.results <- function(x) { #resets vectors which will be filled i <- 0 AICA <- NA; #put models names hierarchically in vector modelnames <- c("1=global", "2=biotic1", "3=biotic2", "4=abiotic") #keep track of changes in model names and number for (i in 1:length(modelnames)) #model structure of the four models given for all models to run #global modelstructure <- c( "ZlogHRS ~...

...rectly into the table using mysql (for live it''ll be about .5mil rows) and configured a Rails search page along the exact same lines as ferret_demo. Unfortunately, search wouldn''t return any results; following the advice given to other posters, I deleted the index/development/[ModelName] folder. I expected then for the index to automatically get rebuilt, but unfortunately any search then resulted in: : Error occured at <fs_store.c>:318 Error: exception 2 not handled: Couldn''t open the file to read vendor/plugins/acts_as_ferret/lib/acts_as_ferret.rb:284:in `sear...

...s.data.frame(data)) Hit <Return> to see next plot: Hit <Return> to see next plot: Hit <Return> to see next plot: Error in 1:L : NA/NaN argument In addition: Warning messages: 1: best model occurs at the min or max # of components considered in: summary.mclustBIC(Bic, data, G = G, modelNames = modelNames) 2: optimal number of clusters occurs at min choice in: Mclust(as.data.frame(anc.st.mat)) 3: insufficient input for specified plot in: coordProj(data = data, parameters = x$parameters, z = x$z, what = "classification", That's puzzling because the example given by ?Mc...

...fine, the parameter estimates are reasonable # across 3 chains samplesBgr("beta") # crash samplesBgr("beta", plot=FALSE) # also crash 2. I'd like to run a bunch of models via BRugs, saving the state after each one, like so: for (m in models) { fit <- BRugsFit(paste(m$modelname, ".bug"), ...) save(fit, paste(m$modelname, ".RData)) } After I load the saved fitted objects in a new R process, BRugs functions like samplesStats() no longer work. Is it possible to somehow load the saved objects and use samplesStats and other BRugs functions?

...sets. R and Emacs freeze, my SystemMonitor in Ubuntu shows that the CPU grinds to a halt but the memory is exhausted. Shouldn't a process in R exit and return an error when memory can no longer be allocated? Wishful thinking???? Mark > exprs.clust <- Mclust(exprs(AOP.sig), G=9:12, modelNames="VEI") Warning messages: 1: In summary.mclustBIC(Bic, data, G = G, modelNames = modelNames) : best model occurs at the min or max # of components considered 2: In Mclust(exprs(AOP.sig), G = 9:12, modelNames = "VEI") : optimal number of clusters occurs at max choice >...

My code executes fully: <% for column in ModelName.content_columns %> <th><%= column.human_name %></th> <% end %> but it does not pick up any columns in my model''s table that end in "_id". Is this by design? Is there a helper method that will reference the related tables? Or will I need to...

...to different cpu types, research i should be able to mask one can anyone confirm this is doable with these two boxes, and if so which should be masked to which ?? xe host-cpu-info cpu_count : 8 vendor: GenuineIntel speed: 2493.790 modelname: Intel(R) Xeon(R) CPU L5420 @ 2.50GHz family: 6 model: 23 stepping: 10 flags: fpu de tsc msr pae mce cx8 apic sep mtrr mca cmov pat clflush acpi mmx fxsr sse sse2 ss ht nx constant_tsc aperfmperf pni vmx est ssse3 s...

...executed installation by a standard command emerge nut. All has passed successfully. I use UPS POWERCOM BNT-1500AP. It is connected to serial port 1, i.e. for my system it/dev/tts/0. In/etc/nut/ups.conf [powercom] driver = powercom port = /dev/tts/0 desc = "PowerCom BNT-1500AP on IP-PBX" modelname = "BNT-1500AP" linevoltage = 230 type = KIN1500AP also trying: [powercom] driver = powercom port = /dev/tts/0 desc = "PowerCom BNT-1500AP on IP-PBX" modelname = "BNT-1500AP" linevoltage = 230 type = Egys At start by a command/etc/init.d/upsdrv start gives out the...

...the asymmetric marginal distributions and then simulating from those distributions. For this I have tried to build a function as below "mclustSim" <- function(x){ library(mclust) nr=dim(x)[1] nc=dim(x)[2] NumberOfSim=100 for (i in 1:nc){ y<-Mclust(as.matrix(x[,i])) z<-sim(modelName=y$modelName,parameters=y$parameters,n=NumberOfSim) #print(c(colnames(x[i]),z[,2]),digits=3) } } z is the simulated return series. How can I say: "print z next to previous z and form a matrix" How could I get an matrix as a result..? So that I would have simulations from asset1 distribu...