similar to: Estimating a Normal Mixture Distribution

Dear R user I follow the steps defined in Modern applied statistics page(453) to use optim. However, when I run the following code the parameters seems way off and the third parameter(p3) stayed as the initial value. below is the code: ## data da=c(418,401,416,360,411,425,537,379,484,388,486,380,394,363,405,383,392,363,398,526) ### initial values pars=c(392.25, 507.25, 0.80)

Hey, I am confused about how to generate the sample data from a mixture of Multivariate Gaussian ditribution. For example, there are 2 component Gaussian with prior probability of 0.4 and 0.6, the means and variances are u1=[1 1]', Cov1=[1 0;0 1] and u2=[-1 -1]', Cov2=[1 0;0 1] repectively. So how can I generate a sample of 500 data from the above mixture distribution? Thanks. Fred

MailMarshal has stopped the following message: Message: B0000edbce.00000001.mml From: samba@samba.org To: 0HMN001SJJEO7P@meandermc.nl; 0HMN002PGJ3XTO@meandermc.nl; 0HMN003MSRVTJ4@meandermc.nl; 0HMN004HCKM5BS@meandermc.nl; 0HMN00CLFJAHEI@meandermc.nl; 0HMU001QIAZXD7@meandermc.nl; 0HMU00DMEHZEM8@meandermc.nl; 0HN50004MEFYL8@meandermc.nl; 3Darjan@meandermc.nl;

Hi, in the code below, I am drawing 1000 samples from two beta distributions, each time using the same random number seed. Using set.seed(80) produces results I expect, in that the differences between the distributions are very small. Using set.seed(20) produces results I can't make sense of. Around half of the time, it behaves as with set.seed(80), but around half of the time, it behaves

Hi All I need some real help in using virt-install. I have a fedora8 o/s + Xen installed and wanted to have RHEL4 as the Guest Domain. I have mounted the CD1 of RHEL4 and wanted to install this as guest domain using virt-install and not able to do and seeking help. Please let me know if the parameters to virt-install are correct and if not please let me know the correct parameters. Contents of

Any help is greatly appreciated.   Thanks Chandra Sekhar --- On Sun, 11/16/08, Chandra Sekhar <sekhar_15@yahoo.com> wrote: From: Chandra Sekhar <sekhar_15@yahoo.com> Subject: [Xen-users] Error while running virt-install To: xen-users@lists.xensource.com Date: Sunday, November 16, 2008, 11:23 PM Hi All I need some real help in using virt-install. I have a fedora8 o/s + Xen

Dear all, I would like to solve the following problem, which can be done with optimal control theory or dynamic programming: max(x,y) a*u1+b*u2+c*f1(u2) s.t. 0<u1<x, 0<u2<f2(x,u2), x'=f3(u1,u2,x) which can be rewritten if optimal control theory should be applied as H=a*u1+b*u2+c*f1(u2)+lambda*(x') s.t. 0<u1<x, 0<u2<f2(x,u2) The maximum principle

Hi, In grid 4.0, adding two complex units by `sum()` seems to give wrong results. In the following example, `u1 + u2` gives the correct result, but `sum(u1, u2)` also `sum(unit.c(u1, u2))` give the wrong results. ``` library(grid) u1 = 0.4*sum(unit(1, "inch"), unit(1, "mm")) u2 = 0.1*sum(unit(1, "inch"), unit(1, "mm")) u1 # [1] 0.4*sum(1inches, 1mm)

Hello! I think I've found a bug in some (older) versions of rsync. The problem does not seem to have been reported at bugzilla. And of the 9 bugs found there by searching on "--delete-after" all have got the name wayned@samba.org attached. I ask you therefore to check if this problem is found in the newest version of rsync as well! I've read the warnings at rsync.samba.org

Hi, If I have an R object UUU, where the second element is U2, based on "g" column of my.table my.table of UUU is: mmm ggg gindex map Info aaa123 U1 1 1 1 aaa124 U1 1 2 1 bbb1378 U2 2 1 1 bbb8888 U2 2 2 0 bbb1389 U2 2 3

-----BEGIN PGP SIGNED MESSAGE----- Hash: SHA1 > Message: 2 > Date: Tue, 19 Aug 2003 14:05:40 +0200 > From: samba-list@naev.de (Peter Koch) > Subject: [Samba] Samba creates User-ACL's > To: samba@lists.samba.org > Message-ID: <200308191205.OAA02625@intranet.naev.de> > > Dear Readers: > > I'm using Samba 2.2.8a with ACL-support and noticed the >

This is a function I coded a few years ago to calculate a confidence interval for a coefficient of variation. The code is based on a paper by Mark Vangel in The American Statistician. I have not used the function much, but it could be useful for comparing cv's from different groups. Kevin Wright confint.cv <- function(x,alpha=.05, method="modmckay"){ # Calculate the

Hi everybody, I try to simulate random numbers from a trivariate nested Archimedean copula. My aim is to correlate two processes with, e.g. theta2, as the so called child pair and then to correlate these two processes with a third one with theta1 (parent). This "figure" tries to capture what I am explaining theta1 theta2

> On Mar 11, 2018, at 3:32 PM, Neha Aggarwal <aggarwalneha2000 at gmail.com> wrote: > > Hello All, > I am facing a unique problem and am unable to find any help in R help pages > or online. I will appreciate your help for the following problem: > I have 2 data-frames, samples below and there is an expected output > > R Dataframe1: > C1 C2

Hello everybody! I’m trying to define the optimal number of surveys to detect the highest number of species within a monitoring season/session. To do this I want to run all the possible combinations between a set of samples and to calculate the total number of species for each combination of 2, 3, 4 …n samples events, so that at the end I will be able to define which is the lowest number of

There are two users u1 & u2 in the cluster. Some files are created by u1, and they are read only for u2. Of course u2 can read these files. Later these files are renamed by u1. Then I switch to the user u2. I find that u2 can't list or access the renamed files. I see these errors in log: [2018-01-15 17:35:05.133711] I [MSGID: 109045] [dht-common.c:2393:dht_lookup_cbk] 25-data-dht:

Dear list! I am interested in learning about MLE and I wonder whether it is possible to use the examples for maximum likelihood estimation given in 8.7 in MASSv3 with R? AFAIU R does not have a direct replacement for S-PLUS's ms() which the examples use for the fitting, but optim() may be of help for me. However, I am not sure how I can convert the use of ms() using optim()--is there an

On 2016-12-07 4:08 am, truename wrote: > OS: CentosOS 7 > > I have installed samba + openldap + smbldap-tools + pam by: > > yum --enablerepo=extras install -y epel-release > yum install -y smbldap-tools > yum install -y samba openldap openldap-clients openldap-servers > migrationtools > yum install -y nss-pam* > > I know that smbldap-tools is a

OS: CentosOS 7 I have installed samba + openldap + smbldap-tools + pam. I create users and groups by: sudo smbldap-groupadd -a g1 sudo smbldap-groupadd -a g2 sudo smbldap-useradd -a -P -g 1001 u1 sudo smbldap-useradd -a -P -g 1002 u1 groups u1 u1 : g1 groups u2 u2 : g2 I create shared folder: cd / sudo mkdir data sudo chown u1 /data sudo chgrp g1 /data sudo chmod ugo+rwx /data I edit

Hi, Our department has detected a bug in the implementation of the Kinderman-Ramage generator for normal random variates in version 1.7.0, which can be seen from the below R session. (Consecutive calls for chisq.test(...) always gives p-values very close to 0.) We have already encountered this bug in version 1.6.2 The error is in file R-1.7.0/src/nmath/snorm.c Here is a patch for this file to