similar to: Help MLE

I have a data with 236 observations. After plotting the histogram, I found that it looks like non-central t distribution. I would like to get MLE for mu and df. I found an example to find MLE for gamma distribution from "fitting distributions with R": library(stats4) ## loading package stats4 ll<-function(lambda,alfa) {n<-200 x<-x.gam

I am just trying to teach myself how to use the mle function in R because it is much better than what is provided in MATLAB. I am following tutorial material from the internet, however, it gives the following errors, does anybody know what is happening to cause such errors, or does anybody know any better tutorial material on this particular subject. >

Hi, How can we set the environment for the minuslog function in mle()? The call in this code fails because the "ll" function cannot find the object 'y'. Modifying from the example in ?mle: library(stats4) ll <- function(ymax=15, xhalf=6) { -sum(stats::dpois(y, lambda=ymax/(1+x/xhalf), log=TRUE)) } fit.mle <- function(FUN, x, y) { loglik.fun <- match.fun(FUN)

Dears, I'm trying to find the parameters (a,b, ... l) that optimize the function (Model) described below. 1) How can I set some constraints with MLE2 function? I want to set p1>0, p2>0, p3>0, p1>p3. 2) The code is giving the following warning. Warning: optimization did not converge (code 1) How can I solve this problem? Can someone help me? M <- 14 Y = c(0, 1, 0, 0, 0,

Hi, I have written out the log-likelihood function to fit some data I have (called ONES20) to the non-central chi-squared distribution. >library(stats4) >ll<-function(lambda,k){x<-ONES20; 25573*0.5*lambda-25573*log(2)-sum(-x/2)-log((x/lambda)^(0.25*k-0.5))-log(besselI(sqrt(lambda*x),0.5*k-1,expon.scaled=FALSE))} > est<-mle(minuslog=ll,start=list(lambda=0.05,k=0.006))

There is a thread on this topic already: http://finzi.psych.upenn.edu/Rhelp10/2010-August/250934.html I'm rather mystified by a similar problem and wondering whether I've overlooked something obvious. I'm running with latest versions of R and all packages, and latest version of JAGS running under Windows 7. Here's the problem. I have some source code. It's given below -

Hello, I have successfully installed SAMBA4 Alpha12 as PDC. I added two machines. Ubuntu 11.04 Linux and Windows 7. They both can authenticate and log in to the machines. WHich is great. The problem comes when users trying to change their password. This is not working! It's constantly complaining for the complexity. However I used strong and absolutely different to any previous password.

I am trying to run WinBUGS 1.4 from the Ubuntu 6.10 Linux distribution. I am using the R2WinBUGS packages with the source file listed below. WinBUGS appears to run properly, but I get the following message after WinBUGS starts in WINE. Does anyone know what may be causing this error and what the correction may be? Thanks ERROR MESSAGE: fixme:ole:GetHGlobalFromILockBytes cbSize is 13824

Dear R-help readers, I am writing to you in order to ask you a few questions about distribution fitting in R. I am trying to find out whether the set of event interarrival times that I am currently analyzing is distributed with a Gamma or General Pareto distribution. The event arrival granularity is in minutes and interarrival times are in seconds, so the values I have are 0, 60, 120, 180, and

my problem actually arised with fitting the data to the weibulldistribution, where it is hard to see, if the proposed parameterestimates make sense. data1:2743;4678;21427;6194;10286;1505;12811;2161;6853;2625;14542;694;11491; ?? ?? ?? ?? ?? 14924;28640;17097;2136;5308;3477;91301;11488;3860;64114;14334 how am I supposed to know what starting values i have to take? i get different

Hello! I've set up samba 4 with great success since alpha12, but i've always used a locally installed version of BIND. Is it possible to use a seperate BIND server instead? I'd like to not run BIND on my file server. I have not seen any instructions to this effect, so perhaps someone would please point me in the correct direction? Thank You Very Much!! ++AMARU

Hi All, I am using Samba4 Alpha12, and i am trying to edit the password policy (which I believe is derived from the PDC's policy), but I am trying to figure out whether Samba4 can respond to GPOs. If It cannot, how would I define the password policy for the domain? Thanks.

Dear all, after switching from Samba4-alpha10 to -alpha12, it is no longer possible to manage user accounts via the RSAT on Windows 7 Professional. The program says that name information can't be retrieved, because the domain controller isn't configured properly (or the like, I read it in German). However, joining the domain and logging in with administrator or simple user account is

I am trying to see if Hyper-V 2008 R2 Failover Clustering will work with Samba 4 as the Active Directory server. Have installed Samba 4 (version alpha12) on CentOS 5.4. Anybody else tried this before? Please share your experiences here. Thanks, Ajay

Hallo, I'm working with the mle function and I would like to ask you a couple of questions. My goal is to construct the historical value of v1(t), v2(t) and v3(t) using the maximum likelihood estimation. So, I need to optimize the following log-likelihood: sum(E1_f[t,]*(v1*teta1[] + v2*teta2[] + v3*teta3[]) - E_f[t,]*log(1 + exp(v1*teta1[] + v2*teta2[] + v3*teta3[]))) (E_f and E1_f

Dear all I've got the dataset data:2743;4678;21427;6194;10286;1505;12811;2161;6853;2625;14542;694;11491; ?? ?? ?? ?? ?? 14924;28640;17097;2136;5308;3477;91301;11488;3860;64114;14334 I know from other testing that it should be possible to fit the data with the exponentialdistribution. I tried to get parameterestimates for the exponentialdistribution with R, but as the values of the parameter

-----BEGIN PGP SIGNED MESSAGE----- Hash: SHA1 Hi there, I'm test the Account Lockout Threshold policy (Computer Setting) in samba4-alpha12 but dont work. I the user and Computer to OU and nothing. Work fine this policy? I need a password expiry warning too... any suggestion Thanks a lot -- Ivan Martinez -----BEGIN PGP SIGNATURE----- Version: GnuPG v1.4.9 (GNU/Linux) Comment: Using GnuPG

Hello all! I am new to the mailing lists, as well as Samba. I was looking for an Active Directory alternative for linux, and finally landed on Samba4. My installation is alpha12, installed via apt on Ubuntu 10.10. I am trying to create roaming profiles, however, if "path = /UserData" then the profiles are viewable from every user account by simple accessing "

Hi, Actually I am facing a similar problem. I would like to fit both an ordinary (symmetric) and a non-central t distribution to my (one-dimensional) data (quite some values.. > 1 mio.). For the symmetric one, fitdistr or funInfoFun (using fitdistr) from the qAnalyst package should do the job, and for the non-central one.. am I right to use gamlss(x ~ 1, family=GT()) ? Anyway, I am a little

I've run into persistent problems with OpenBUGS crashing when using BRugs .53 and .71, and am hoping someone has suggestions. There is obviously something unusual going on in my environment, but I'm at a loss as to where to begin to try to solve it. In a nutshell, what happens is that, as soon as I call "modelCheck()" in BRugs, it gets an error or crashes ... but only some of