[R] rlaplace using rmutil - HELP

2010-02-15 Thread vaibhav dua
Hi, I'm using rlaplace distribution to draw non normal random deviates. I'm using rmutil package for this. I know rlaplace needs location and dispersion parameters but I want to pass diag matrix as dispersion parameters and I'm not getting the correct values. Here is my code snippet: disp

Re: [R] rlaplace using rmutil - HELP

2010-02-15 Thread vaibhav dua
ictor From: Peter Ehlers Cc: r-help@r-project.org Sent: Mon, February 15, 2010 10:52:33 AM Subject: Re: [R] rlaplace using rmutil - HELP On 2010-02-15 8:26, vaibhav dua wrote: > Hi, > > I'm using rlaplace distribution to draw non normal random deviates. I'm usi

[R] HELP on Non-Linera Mixed-Effect model

2010-02-15 Thread vaibhav dua
Hi, I'm trying to fit nonlinear mixed effects model using nlme function but getting an error message. Here is what I have: fitted_model = nlme(scores~spline(b1,b2,b3,kt,time), fixed = list(b1~1, b2~1, b3~1, kt~1), random = b1+b2+b3~1, groups= ~id, da

[R] Univariate Normal Mixtures - HELP

2010-02-18 Thread vaibhav dua
Hi, I'm generating a mixture of 2 univariate normal distributions using norMix and rnorMix and would like to put a constraint on Mean (Equal). here is my code snippet: library(nor1mix) X <- norMix(mu=c(50, 60 ), sig2=c(5,10), w=c(0.5, 0.5)) mixData <- rnorMix(1000, X) plot(mixData, type="l"

[R] Generate Univariate Normal Mixtures???

2010-02-19 Thread vaibhav dua
Hi, I'm trying to generate a mixture from 2 univariate normal distributions and would like to put a constrained on the means (Equal) of the two distributions. Here is my shot at the code: library(nor1mix) X <- norMix(mu=c(50, 60 ), sig2=c(5,10), w=c(0.5, 0.5)) mixData <- rnorMix(1000, X) p

[R] remove quotes from the paste output

2010-12-11 Thread vaibhav dua
Hi, I'm generating the name of the variable with paste function and then using that variable name further to get the specific position value from the data.frame, here is the snippet from my code: modelResults <- extractModelParameters("C:/PilotStudy/Mplus_Input/Test", recursive=TRUE) #extract

[R] extractModelParameters HELP!!!

2011-03-05 Thread vaibhav dua
Hello Gurus, I'm using one of the function, i.e. extractModelParameters, from the MplusAutomation package to read parameters fromMplus out files into R. I have 360 cells and each cell has several out files (max. 100). I'm extracting the model parameters from each cell and doing subsequent calcu