similar to: Mixture of Univariate Normals

Hi, I would like apply different link functions from Aranda-Ordaz (1981) family to large binary dataset (n = 2000). The existing links in glm for binomial data (logit, probit, cloglog) are not adequate for my data, and I need to test some other transformations. Is it possible to do this in R? And how? Thank you for your help, /Sharon

Hi R-help, Is there such a thing as a function in R for fitting a GLM where the response is distributed as a Beta distribution? In my case, the response variable is a percentage ([0,1] and continuous). The current glm() function in R doesn't include the Beta distribution. Thank you for any help on this topic. Sincerely, Sharon K?hlmann

Full_Name: Sharon Kuhlmann Version: 0.99 OS: solaris2.6 Submission from: (NULL) (129.16.167.231) When I use the function > eigen(x), it gives me a "Process R bus error", and the program quits. I'm running the following R version: platform sparc-sun-solaris2.6 arch sparc os solaris2.6 system sparc, solaris2.6 status Patched

Hi, I have some trouble with the result of a fortran function calling the R binomial RNG: The C function rbinom is wrapped as in the file attached. My main fortran program starts as call rndstart() and ends as call rndend() I happen to call the binomial RNG within a loop as b = ggrbinom(1.d0,0.5d0) write(*,*) 'b=',b In certain cases, after a few iterations in the loop, b get

We are currently analyzing data on children clustered in day care-centers (DCC). We have tried to use geepack and gee libraries to estimate an overall incidence rate for absences (=number of absences/risk time) by specifying geese(number.absences ~ offset(log(risktime)), id=day.care.id, family=poisson("log"), data=dcc, corstr="exch",

Hi, I am trying to adjust the position of the boxplots when using bwplot(). In boxplot() there is the alternative to modify the position of the boxes by using boxplot(y~x, at=...). However, I don't seem to find a comparable option in bwplot(). I would like in some instances to add more space between the boxplots, e.g. that I could specify at=c(1, 2, 5, 6) so that the first two boxplots are

Hi, I am trying to create a graph using Arial font (as required by PLoS One). I have read probably all posts in R-help on this topic, as wells as R-news 2006. The code I have been trying is following: Arial <- Type1Font(family="Arial", metrics=c("ArialMT.afm", "arial-BoldMT.afm", "Arial-ItalicMT.afm",

All: I am trying to fit a mixture of 2 normals with > 110 million observations. I am running R 2.5.1 on a box with 1gb RAM running 32-bit windows and I continue to run out of memory. Does anyone have any suggestions. Thanks so much, Tim [[alternative HTML version deleted]]

An embedded and charset-unspecified text was scrubbed... Name: not available URL: <https://stat.ethz.ch/pipermail/r-help/attachments/20091125/6a78f78b/attachment-0001.pl>

hi i'm looking for a gibbs sampling algorithm for R for the case of mixture of K normals, and in particular for the case of bivariate normals. i'd be grateful if anyone could send its own R-routine, at least for the univariate case. thank you in advance matteo

Hello, I am trying to find out if R can do the following: I have a mixture of normals say f = 0.2*Normal(2, 5) + 0.8*Normal(3,2) How do I find the difference in the densities at any particular point of f and at Normal(2,5)? -- Thanks, Jim. [[alternative HTML version deleted]]

Hi, I have a problem of estimating a mixture of two normal distributions. I need to find the starting points automatically, since this is a part of a larger piece of image processing code. I found the mix2normal1 function in VGAM package that mentions a method of finding starting values for mu1 and mu2 but refers the reader to a book by Everitt and Hand. Unfortunately, I do not have an easy

Can anyone send me the R code for fitting a mixture of normal distributions. Or in fact any'combination' of mixtures of distributions. I looke dta mix but it is still not clear to me. Jim [[alternative HTML version deleted]]

Hi, I would like to know if it is possible to have a "R code" to estimate the parameters of a mixture of bivariate (or multivariate) normals via EM Algorithm. I tried to write it, but in the estimation of the matrix of variance and covariance, i have some problems. I generate two bidimensional vectors both from different distribution with their own vector means and variance and

Hello, I have a concrete statistical question: I have a sample of an univariate mixture of an unknown number (k) of normal distributions, each time with an unknown mean `m_i' and a standard deviation `k * m_i', where k is known factor constant for all the normal distributions. (The `i' is a subscript.) Is there a function in R that can estimate the number of normal distributions k

Yes. I'm running The 2.3.99.3 win32 version. And if you can to provide the " mime.types " and " crossdomain.xml " files by yourself it will be great. TNX! -----Original Message----- From: "Thomas B. R?cker" [mailto:thomas at ruecker.fi] Sent: Tuesday, October 08, 2013 5:58 PM To: icecast-dev at xiph.org Cc: Yaniv Sharon Subject: Re: [Icecast-dev]

I have an "interesting" problem that I didn't think would ever happen... Using rsync over ssh to mirror a local directory on a remote host. On the remote host the disk becomes full during the rsync run. What I find then is an ssh process on the remote host with no child processes, i.e. the remote rsync has correctly terminated, but the shell has not. On the lost host, rsync is

Hi, I am currently working on fitting a mixture density to financial data. I have the following data: http://s000.tinyupload.com/?file_id=00083355432555420222 I want to fit a mixture density of two normal distributions. I have the formula: f(l)=πϕ(l;μ1,σ21)+(1−π)ϕ(l;μ2,σ22) my R code is: normalmix<-normalmixEM(dat,k=2,fast=TRUE) pi<-normalmix$lambda[1] mu1<-normalmix$mu[1]

So there is no usage for the " mime.types " in ice win32? -----Original Message----- From: icecast-dev-bounces at xiph.org [mailto:icecast-dev-bounces at xiph.org] On Behalf Of Thomas B. R?cker Sent: Tuesday, October 08, 2013 5:09 PM To: icecast-dev at xiph.org Subject: Re: [Icecast-dev] icecast-2.3.99.3 Yaniv Sharon <yaniv.sharon at gmail.com> wrote: >Hi, > >Where to

https://bugzilla.samba.org/show_bug.cgi?id=10092 Summary: Hang when remote dest's disk is full Product: rsync Version: 3.0.9 Platform: All OS/Version: All Status: NEW Severity: normal Priority: P5 Component: core AssignedTo: wayned at samba.org ReportedBy: list0570 at paradise.net.nz