Displaying 20 results from an estimated 900 matches similar to: "Problem on flexmix when trying to apply signature developed in one model to a new sample"
2012 Oct 30
1
Data set BregFix in package flexmix
Dear list:
I would like to recreate how the artificial data set BregFix was generated
in package flexmix (thanks Bettina and Friedrich). The data set is
thoroughly described in Grun and Leisch's Computational Statistics & Data
Analysis 51(11) :5247-5252 but references to the appropriate seed number(s)
are missing (providing these details was certainly beyond the scope of the
authors'
2007 Feb 08
1
smartpred depends on fitted() in flexmix?
Hi,
I was going through the examples in smartpred. It seems there's an
unstated dependency on the fitted() function in package flexmix.
n = 20
set.seed(86)
x = sort(runif(n))
y = sort(runif(n))
library(splines)
fit = lm(y ~ ns(x, df=5))
plot(x, y)
lines(x, fitted(fit)) # won't work w/o prior loading of the flexmix
package.
newx = seq(0, 1, len=n)
points(newx, predict(fit,
2007 Apr 16
1
flexmix glm warning
Dear R-helpers & Friedrich Leisch,
I'm trying to fit a mixture of regression models on attached data set
by doing the following:
agl <- read.table("agl.txt")
mod1 <- flexmix(resp~gng+csc|Subject,data=agl,
model=FLXglm(family="binomial"),k=2)
The result is a (varying) number of the following warnings:
Warning messages:
1: non-integer #successes in a
2009 Nov 09
1
model based clustering with flexmix
Hello all,
I am trying to fit a truncated mixture model and I wrote a driver for
flexmix following the example in the vignette, but it doesn't work for
me: it assigns all data points to one component only, e.g.:
> > source('bugged.R')
>
> Call:
> flexmix(formula = x ~ 1, k = 2, model = truncatedmodel(lower = -4,
> upper = 4))
>
> prior size
2005 Apr 12
1
R Package: mmlcr and/or flexmix
Greetings
I'm a relatively new R user and I'm trying to build a latent class model.
I've used the 'R Site Search' and it appears there's not much dialogue on
these packages
On mmlcr, I've gotten it working, but not sure if I'm using it correctly.
On flexmix, I can only seem to get results for one class.
I'm attaching my code below - if anyone
2011 Sep 01
3
betareg question - keeping the mean fixed?
Hello,
I have a dataset with proportions that vary around a fixed mean, is it
possible to use betareg to look at variance in the dispersion parameter
while keeping the mean fixed?
I am very new to R but have tried the following:
svec<-c(qlogis(mean(data1$scaled)),0,0,0)
f<-betareg(scaled~-1 | expt_label + grouped_hpi, data=data1, link.phi="log",
2013 Apr 23
2
flexmix
Hi,
I am trying to understand how to use the flexmix package, I have read the Leisch paper but am very unclear what is needed for the M-step driver. I am just fitting a simple linear regression model. The documentation is far from clear what the FLXmclust function does, but, it in principle could do all I need, however, I do not get sentible results as if I try the following the result is poor:
2009 Nov 08
2
Simple 2-Way Anova issue in R
Hello, I'm new to R and have been following many guides including the two-way
anova (http://www.personality-project.org/r/r.anova.html). Using that
walkthrough including the supplied data I do str(data.ex2) and receive the
appropriate types of data as follows:
> str(data.ex2)
'data.frame': 16 obs. of 4 variables:
$ Observation: int 1 2 3 4 5 6 7 8 9 10 ...
$ Gender :
2009 Sep 22
3
problem using KLdiv - flexmix
I am receiving the following error while executing KLdiv method:
kl<-KLdiv(y)
Error in function (classes, fdef, mtable) :
unable to find an inherited method for function "KLdiv", for signature
"data.frame"
the code is as follows:
require(flexmix)
KLdiv(y)
the object y (snippet, actual object contains 188 entries) is as follows:
s1p5 s2p5
1 0.6 0.8
2 0.8
2011 Jan 12
1
how to change strip text of effect plot
Dear r heper,
How can I change the strip text, for example (16,23] in the following
example, to other more informative text such as "high level" on the
fly?
library(effects)
Cowles$ex2 <- cut(Cowles$extraversion,3)
mod.cowles <- glm(volunteer ~ sex+neuroticism*ex2,data=Cowles, family=binomial)
eff.cowles <- allEffects(mod.cowles)
plot(eff.cowles,
2009 Jan 16
1
installing mclust and flexmix on linux
I've been trying to install some R packages such as mclust and flexmix on linux but have had the following error messages.
> I've been trying to install mclust on my notebook which has linpus linux lite os and I have installed R as well as some packages all right. However, when I tried to install mclust, it gave me the following messages. Any suggestions?
Tim
++++++++++++++++
1997 Nov 03
1
R-alpha: expression(..) objects -- c(ex1, ex2) `fails'
When investigating the
legend(x,y, expression ( <math>, <math> ), ....)
[as suggested by Jim Lindsey],
I once again stumbled over the fact that expression objects behave somewhat
``wierdly'', IMHO.
We can have expressions of length(.) > 1, and access elements of these like
list elements.
However, other natural things don't work with expressions; some example
2011 Sep 12
1
findFreqTerms vs minDocFreq in Package 'tm'
I am using 'tm' package for text mining and facing an issue with finding the
frequently occuring terms. From the definition it appears that findFreqTerms
and minDocFreq are equivalent commands and both tries to identify the
documents with terms appearing more than a specified threshold. However, I
am getting drastically different results with both. I have given the results
from both the
2010 Aug 28
1
star models
Hi,
I am traying to implement an STAR model, but I have some problems.
I am following the instruction of the model, that they are in:
http://bm2.genes.nig.ac.jp/RGM2/R_current/library/tsDyn/man/star.html
that they are from:
http://bm2.genes.nig.ac.jp/RGM2/pkg.php?p=tsDyn
The model is:
star(x, m=2, noRegimes, d = 1, steps = d, series, rob = FALSE, mTh,
thDelay, thVar, sig=0.05, trace=TRUE,
2010 May 10
1
how to get p-value from ave
Hi there,
I checked google for aov. usually one uses summary to see whether the p-value is small.
but I want to put aov in my script. how can I get the p-value, (0.1115, 0.6665, 0.6665 in the following example)?
thanks
YU
> datafilename="http://personality-project.org/r/datasets/R.appendix2.data"
> data.example2=read.table(datafilename,header=T)
> aov.ex2 =
2010 Apr 03
2
Using ifelse and grep
Good Morning,
I am trying to create a new column of character strings based on the first
two letters in a string in another column. I believe that I need to use some
combination of ifelse and grep but I am not totally sure how to combine
them. I am not totally sure why the command below isn't working. Obviously
it isn't finding anything that matches my criteria but I am not sure why.
Any
2012 Jan 18
2
Table Intersection
I've got two tables....
first one(table1):
ID chrom start end
Ex1 2 152 180
Ex2 10 2000 2220
Ex3 15 3000 4000
second one ( table2):
chrom location name
2 160 Alv
2 190 GNN
2 100
2018 Feb 15
1
RV: Problem_graphic
Hi everyone,
I am beginner using R but I try to learn more.
I need this graphic or similar but instead of tropical and temperate are three ontogenetic states
Larva Met Juv
CTMAX
CTMIN
SP SP SP
This is my scrip:
Ex = subset(Expr, Outlayer=="N")
Ex2 = subset (Ex, S0 == 1)
Ex3
2010 Apr 14
2
GAMM : how to use a smoother for some levels of a variable, and a linear effect for other levels?
Hi,
I was reading the book on "Mixed Effects Models and Extensions in
Ecology with R" by Zuur et al.
In Section 6.2, an example is discussed where a gamm-model is fitted,
with a smoother for time, which differs for each value of ID (4
different bird species). In earlier versions of R, the following code
was used
BM2<-gamm(Birds~Rain+ID+
2007 Jan 04
1
problem with function 'optimise' (PR#9438)
Full_Name: Karsten Krug
Version: 2.4.0
OS: Open Suse 10.0, Windows XP
Submission from: (NULL) (88.134.13.50)
I found a problem in the 'optimise' function for one dimensional optimisation.
Example 1:
Try to find a maximum of the function below with the use of 'optimise' in the
interval [0,0.5]. The function follows a parabola and has two local maxima
located at the margins of