similar to: about EM algorithm

Dear R experts, i have a vector z , i have to do something after z is sorted. how can i find the original index, i.e., before sorting, of a certain element in the sorted vector . thanks in advance regards ----------------------------------------------- Peng Jiang ?? ,Ph.D. Candidate Antai College of Economics & Management ???????? Department of Mathematics ??? Shanghai Jiaotong

Hi, Is there any package that I can use to simulate the Gaussian Mixture Model , which is a mixture modeling method that is widely used in statistical learning theory. I know there is a mclust, however, I think it is a little bit different from my problem. Thanks very much.. regards. -------------------------- Peng Jiang ?? Ph.D. Candidate Antai College of Economics &

Hi, dear R experts , I am new. I met this problem when I am trying to learn how to use the nlminb() function. I tried the example which the document provides ( as the following code ) and R gives no response . I don't know whether it is running or not and it takes a very long time but still output nothing so I just close the session window. my question is is there any method

Hi, I first report the installing error to the r-sig-mac mailing list but it seems nobody ever encounter this annoying problem and I got no replies. so I am trying to forward it to this list, hope it never bothers. Anyone had the same problem, or what does the warning message mean ? Thank you in advance ! Begin forwarded message: > From: Peng Jiang <jp021@sjtu.edu.cn> > Date:

Hi,Dear all R experts, I am trying to do the 2-way contingency table analysis by fitting the loglinear models. However, I found my table has several empty cells which are theoretically missing values.I have no idea of how to solve them coz we cannot compute the simulated p-value with zero marginals.Does someone have some suggestions? Please help me out, thanks a lot! Cheers, Yan

Hi all, I am a biological student and need your help in statistics. I have two sets of binomial distributed numbers: {a1, a2, ..., an} and {b1, b2, ..., bn}. How can I get the distribution of the ratios of the two sets of numbers {a1/b1, a2/b2, ..., an/bn}? Is there a formula to transform the distributions? Or where I can start if I want to learn necessary techniques? Thank you a lot! xiechao

Hi, Is the following function built in somewhere? concat = function(v) { res = "" for (i in 1:length(v)) res = paste(res,v[i],sep="") res } e.g. > concat(c("12","3","45")) [1] "12345" Cheers, Andy -- Andy Fugard, Postgraduate Research Student Psychology (Room F3), The University of Edinburgh, 7 George

I'm trying to do realize the following: I have 4 condtions. If all conditions are satisfied I will paste("PASS") If any of these is not satisfied I will paste("FAIL"). But I have to paste the corresponding failure. ifelse is a good solution but for a 2 conditions. Maybe switch or something like this. Does anyone have an idea how to do? Thanks in advance. Adel -- View

Hello there. Is there any function in R that can do cluster on a set of data that has both categorical and numerical variables? thanks. siangli

Dear R Users, I would like to plot three time series on the same graph, two axis on the left and one axis on the right. The time series that I am graphing on the left do not share a similar scale: one has a range of 1:100 and the other a range of 25000:70000. How can I display the tick marks for both on the left hand side without superimposing them, i.e. one set of tick marks in black right

Hello. Somebody knows how can I generate a set of n random vectors of a normal distribution of several variables? For example, I want to generate n=100 random vectors of two dimensions for a normal with mean c(0,1) and variance matrix: matrix(c(2,1,1,3),2,2). Thanks in advance, Arnau.

Dear R listers, I posted this problem several days ago but it seems nobody answered. I use nlm to optimize a given function ,but it always generates the following warnings " Error in nlm(foo, theta.start) : non-finite value supplied by 'nlm' " I don't know why ,can anybody give me some hints ?? thanks in advance. regards .

Dear list, hi ! I am a R beginner and I have a function to optimize . alpha = argmin{ f(x,alpha) } I want alpha to be in [0,1]. Is there any function that can work? I use nlm() but i can't fix the domain of alpha. thanks in advance _______________________ Jiang Peng, Ph.D. Candidate Department of Mathematics & Antai college of Economics and Management Shanghai Jiao

Dear R listers, I have a very annoying problem using nlm(). I want to find the minimizer of my target function, if written in \LaTeX is f(\mu1,\mu2,\sigma1,\sigma2) = \sum_i^n( w_ig_t(z_i) ), where g_t(z) is a pdf of bivariate normal distribution and z_i is my samples. I cannot get the estimation result generated by nlm(), and I got the following errors " Error in

Hello, I have installed samba in a Linux RedHat 8.0 PC in my lab as a file server. The clients of samba are PCs with MS Windows98. Now I meet some problems, I can see samba server name in Network Neighborhood, but when I click the icon, I get this error message after a long delay: Cannot link to Linuxserver(server name). And I used Ethereal to track the packets to find the problem, I found that

Hi, I am Ken.. I am trying to add dud(driver update disk) to centos5.6 using network method, I can successfully add dud to centos5.6 using http and ftp method, but fail when using nfs method. However those all three methods are officially supported in the centos5.6. As indicated in the centos5.6 website, we can see they not only support http, ftp, but also nfs.

After PrologEpilogCodeInserter I found that the instructions that restore callee saved registers S0,S1,LR are in the wrong location ,the instructions are: %S0<def> = LD %SP, 36; mem:LD4[FixedStack2] %S1<def> = LD %SP, 40; mem:LD4[FixedStack1] %LR<def> = LD %SP, 44; mem:LD4[FixedStack0] (LR is the Return address register) the whole code of print-machineinstrs are: # After

Hi, Guys I am starting looking the source code of Speex. Maybe a little silly, but I can not find the implementation of the function speex_encode(). Only something like below void speex_encode(void *state, float *in, SpeexBits *bits) { (*((SpeexMode**)state))->enc(state, in, bits); } Does that mean some of the detail implementation is embedded in lib without source code? I am new in

>Date: Wed, 23 Sep 1998 21:40:38 -0500 >From: "Daryl L. Biberdorf" <darylb@superserve.com> >To: samba@samba.anu.edu.au >Subject: Pages printed >Message-ID: <3.0.5.32.19980923214038.00816310@pop.a001.sprintmail.com> > >We have the need to track printer usage (pages used). >Right now, my Samba printer share defines the print command >as "print

Dear Prof Brian Ripley: I am very sorry for that I use the unappropriate words "It's my great pleasure" when I read this mail again! I give my apology sincerely to you ,unreservedly and with my heart! You are such a warm-hearted man.And you help me to solve my problems.But owing to my poor English, I found I had made such a foolish mistake. My brain must have been stupid