similar to: Problem in optimization of Gaussian Mixture model

Please keep the list in cc. Sorry, it didn't work as expected. Maybe someone else have an appropriate solution. Best, Ulrik On Sa., 26. Aug. 2017, 12:57 niharika singhal <niharikasinghal1990 at gmail.com> wrote: > Hi > > Thanks for you mail, > I really appreciate your time on my problem > > I have posted this problem on > > >

Hi, Thanks for your mail, and time It is not working for some arguments, when mean value is like >6. case mc0 <- c(0.08844446,0.1744455,0.1379778,0.1209769,0.1573065,0. 1134463,0.2074027) rv <-UnivarMixingDistribution(Norm(486.4255, 53.24133), Norm(664.0713, 3.674773), Norm(669.0484, 4.101381),

Hello, I have some parameters from Mclust function. The parameters are in the form *parametersDf * * mu_1 mu_2 var_mc1 var_mc2 c1 c2 * *2 1.357283 2.962736 0.466154 0.1320129 0.5258975 0.4741025 * *21 8.357283 9.962736 0.466154 0.1320129 0.5258975 0.4741025 * Each row in the above data frame

this is probably not really a R specific question, if so apologies for off-topic posting: I'm interested in the probability density function of the maximum values from repeated samples of size N from a normal distribution: smp <- rnorm(N, meanval, stdev) with some mean 'meanval' and standard deviation 'stdev'. I would like to know what is the frequency distribution of

Dear all, I want to use the deltavar() function from emdbook. I can use it directly from the command terminal but within a function it behaves weird. Working example: ---------------------- library("emdbook") fn <- function() { browser() y <- 2 print(deltavar(y*b2, meanval=c(b2=3), Sigma=1) ) } x <- 2 print(deltavar(x*b1, meanval=c(b1=3), Sigma=1) ) y<-3 fn()

Hello, I have got a question. I've got a matrix (mail end) with the colnames x, y, z. In this matrix are different measurements. x and y are risign coordinates. My question. Always, if the "x" AND "y" coordinates are the same, I want to get the mean of their z values. e.q. " x" AND "y" in line1 and line8 are identical: 29 4.5 --> mean of

I haven't been an R lister for a bit, but I hope to enlist someone's help here. I think this is a simple question, so I hope the answer is not much trouble. Can you please respond directly to this email address in addition to the list (if responding to the list is warranted)? I have a histogram and I want to see if the data fit a Poisson distribution. How do I do this? It is

I have a time series with interval of 2.5 minutes, or 24 observations per hour. I am trying to find a 1 hr moving average, looking backward, so that moving average at n = mean(n-23 : n) The time series has about 1.5 million rows, with occasional gaps due to poor data quality. I only want to take a 1 hour moving average for those periods that are complete, i.e. have 24 observations in the

I have a Gaussian mixture model with some parameters mean=(506.8644,672.8448,829.902) sigma=(61.02859,9.149168,74.84682) c=(0.1241933, 0.6329082, 0.2428986) And the plot look something like below.[image: enter image description here] <https://i.stack.imgur.com/4uUQ9.png> Also, if I change my parameters to mean=(2.15,2.0,2.9) sigma=(0.1,0.1,0.1) c=(1/3,1/3,1/3) Then plot would change

I have data that looks like O.lengthO.age 176 1 179 1 182 1 ... 493 5 494 5 514 5 606 5 462 6 491 6 537 6 553 6 432 7 522 7 625 8 661 8 687 10 704 10 615 12 (truncated) with a simple VonB growth model from within nls(): plot(O.length~O.age, data=OS) Oto = nls(O.length~Linf*(1-exp(-k*(O.age-t0))), data=OS, start=list(Linf=1000, k=0.1, t0=0.1), trace=TRUE) mod <- seq(0, 12)

Hi, I would like to have 3d plot, and I found that there is a command "scatterplot3d". x is V, y=sigmaV, and z=TSE (note: After some calculation, I update TSE). How could I do to have 3d plot I want? V<-seq(1279,1280,,100) sigmaV<-seq(0.28,0.29,,100) TSE<-matrix(0, length(V),length(sigmaV)) ps: TSE[i,j] corresponds to V[i] and sigmaV[j]. Thanks, Kate [[alternative HTML

Using R-2.12.1 and latticeExtra-0.6-14, I would like to understand why a rootogram displaying samples from the Poisson distribution looks like I expected it, whereas a rootogram using the normal distribution does not: library(latticeExtra) rootogram(~rpois(1000, lambda = 50), dfun = function(x) dpois(x, lambda = 50)) rootogram(~rnorm(1000), dfun = function(x) dnorm(x,mean(x),sd(x))) I

> On 26 Aug 2017, at 16:39, niharika singhal <niharikasinghal1990 at gmail.com> wrote: > > Hi, > > Thanks for your mail, and time > > It is not working for some arguments, when mean value is like >6. > > > case > > mc0 <- c(0.08844446,0.1744455,0.1379778,0.1209769,0.1573065,0. > 1134463,0.2074027) > > rv

Hello all, Can someone tell me why the following mixture of two log-normal distributions does not get truncated? What puzzles me is that the function works almost always, but for certain combinations (like the one below), it does not. # R code example library(distr) mix<-UnivarMixingDistribution(Lnorm(3.2,0.5),Lnorm(5.4,0.6),mixCoeff=c(0.3,0.7))

Hello, I am having some problem in retrieving the list elements Code Mean_of_mc is a list as described below mean_mc_1 mean_mc_2 mean_mc_3 mean_mc 0.9577246 2.000516 2.962521 mean_mc 7.957725 9.000516 9.962521 optim(par=c(lVal=2, mu=Mean_of_mc), fn=xnorm, method = "BFGS", control=list(fnscale=-1)) xnorm= function(par) { print(par)

------------------------------------------------------------------------------- New versions released for the distr family of package ------------------------------------------------------------------------------- We would like to announce the availability on CRAN (with possibly a minor delay until on every mirror) of new versions of our packages in the "distr"-family (version 2.1),

------------------------------------------------------------------------------- New versions released for the distr family of package ------------------------------------------------------------------------------- We would like to announce the availability on CRAN (with possibly a minor delay until on every mirror) of new versions of our packages in the "distr"-family (version 2.1),

---------- Forwarded message ---------- From: niharika singhal <niharikasinghal1990 at gmail.com> Date: Sun, Aug 27, 2017 at 11:57 AM Subject: Re: Find maxima of a function To: "David Winsemius [via R]" <ml+s789695n4745009h56 at n4.nabble.com>, "Ismail SEZEN [via R]" <ml+s789695n4744993h60 at n4.nabble.com>, Ulrik Stervbo <ulrik.stervbo at gmail.com>

Hi, gfortran trunk and 9-branch now have an option to automatically generate C prototypes for old-style F77 procedures. I just did for a in *.f; do gfortran -fsyntax-only -fc-prototypes-external $a > ${a%.f}.h; done in the src/modules/lapack directory. This generates header files which contain prototypes like int ilaenv_ (int *ispec, char *name, char *opts, int *n1, int *n2, int *n3,

Hello list, I have a problem with the value of the result size ($msize in omegascript) when collapsing results. The index contains 151452 documents. I'm using Omega 1.0.18 on FreeBSD (I tried both the version in ports and the latest one from xapian.org). This is my indexscript: uniqueid: boolean=Q unique=Q field=uniqueid objectid: field=objectid boolean=XID value=0 objecttype: field=type