similar to: NLS fit for exponential distribution

Dear R-help, Can anybody tell me which R package has Lo-Mendell Rubin LR test and Bootstrap LR test to compare the model fit between k class and k+1 class model for Latent class analysis? Thanks in advance,  warn regards,Ms.Karunambigai M PhD Scholar Dept. of Biostatistics NIMHANS Bangalore India [[alternative HTML version deleted]]

Dear R listers, I know I'm breaking the rules by asking a "homework" related question-- I hope you'll forgive me. I am a social psychology graduate student, and the only one in my department who uses R. I successfully completed my multiple regression and structural equation modeling courses using R (John Fox's car and sem packages were a big help, as was his book).

I'm hoping that this is a relatively easy question for someone familiar with the lme4 package. I'm accustomed to using HLM software and writing a simple 2 level [null] equation like this: L1 - Yij = b0 + e L2 - b0 = B00 + u0 The following command in R provides results that are identical to the HLM program. results <- lmer( Y ~ 1 |id , PanelData4) I can't

Hi there, I am trying to fit the following model with a sum of exponentials - y ~ Ae^(-md) + B e^(-nd) + c the model has 5 parameters A, b, m, n, c I am using nls to fit the data and I am using DEoptim package to pick the most optimal start values - fm4 <- function(x) x[1] + x[2]*exp(x[3] * -dist) + x[4]*exp(x[5] * -dist) fm5 <- function(x) sum((wcorr-fm4(x))^2) fm6 <- DEoptim(fm5,

Hi there, I would like to solve a simple equation in R a^2 - a = 8.313 There is no real solution to this problem but I would like to get an approximate numerical solution. Can someone suggest how I can set this up? Thanks in advance, Diviya [[alternative HTML version deleted]]

Dear all, I have a question regarding setIs and method dispatch in S4 classes: Consider the following example: ##################################################### ## A02 "contains" A01 by setIs: setClass("A01", representation(a="numeric",b="numeric")) setClass("A02", representation(a="numeric",d="numeric"))

On 01/08/2013, at 24:37, Gene Cumm <gene.cumm at gmail.com> wrote: > On Wed, Jul 31, 2013 at 7:36 AM, Daniel O'Connor <doconnor at gsoft.com.au> wrote: >> Hi, >> I am trying to create a hybrid ISO image for FreeBSD 9.2 using syslinux 6.01 and isohybrid.pl I found at http://www.overclockix.com/sources/isohybrid.pl (I know there is a C version in syslinux but I

Hello, I get some file names by list.files(). These names are in alphabetical order. I want to change it to logical numeric order. Example: > fileNames <- c("A10", "A1", "A2", "B1", "B2", "B10") > sort(fileNames) [1] "A1" "A10" "A2" "B1" "B10" "B2" I want to have:

--- Matt Zimmerman <mdz@debian.org> wrote: > On Fri, Jul 20, 2001 at 05:15:21PM -0700, Josh Coalson wrote: > > > I rethought it and it seemed like a bad idea to post > > a big file, so you can get it here: > > > > > http://cvs.sourceforge.net/cgi-bin/viewcvs.cgi/flac/junk/flac-1.0-src-candidate.tar.bz2 > > With the attached patch, a complete

Hello there, Can anyone point me to the code for logLik of an nls object? I found the code for logLik of an lm but could not find exactly what function is used for calculating the logLik of nls function? I am using the nls to fit the following model to data - Model 1: y ~ Ae^(-mx) + Be^(-nx) +c and want to understand what is the likelihood function used by nls. Presumably it is using -

Hi there, I have a complex math equation which does not have a closed form solution. It is - y <- (p*exp(-a*d)*(1-exp((d-p)*(a-x[1]))))/((p-d)*(1-exp(-p*(a-x[1])))) For this equation, I have all the values except for x[1]. So I need to solve this problem numerically. Can anyone suggest an optimization package that I can use to estimate the value for x[1]? Thanks in advance, Diviya

Hello there, I have a question regarding bar plots. I am trying to plot the data from the following matrix as a barplot - # input data mdat <- matrix(c(0.1,0.9,0.9,0.1,0.5,0.5,0.45,1-0.45,0.6,0.4,0.8,0.2), nrow = 6, ncol=2, byrow=TRUE, + dimnames = list(c("Mon", "Mon", "Tues", "Tues", "Thurs", "Friday"), +

Hello there, I have a matrix with some data and I want to split this matrix based on the values in one column. Is there a quick way of doing this? I have looked at cut but I am not sure how to exactly use it? for example: I would like to split the matrix "a" based on the spending such that the data is binned groups [0..99],[100..199]...and so on. a <- data.frame(patient=1:7,

Hello there, I am new to using regression in R. I wanted to solve a simple regression problem where I have 2 equations and 2 unknowns. So lets say - y1 = alpha1*A + beta1*B y2 = alpha2*A + beta2*B y1 <- runif(100000, 0,1) y2 <- runif(100000,0,1) alpha1 <- 0.6 alpha2 <- 0.75 beta1 <- 1-alpha1 beta2 <- 1-apha2 I now want this equation to estimate the values of A and B. Both A

Hi there, I was wondering if there is any R package that one can use for plotting that has more legend symbols - the standard pch has 18 symbols but I need ~30 for my application- and just using different colors is not an option. Thank you in advance, Diviya [[alternative HTML version deleted]]

Hello there, I recently wrote some code to perform pairwise correlations between all samples in a large dataset. So we are talking about performing pairwise correlations between 400K vectors. Since R has a very rich library of functions, it was very easy to code this in R. However, R was probably not the best choice as it is super slow for this large job. So my plan is to recode it in C. I was

CentOS 8 Samba Version 4.12.5 compiled from source configured as DC and joined to domain SAMDOM.EXAMPLE.COM. smb.conf -- # Global parameters [global] netbios name = B10 realm = SAMDOM.EXAMPLE.COM server role = active directory domain controller workgroup = SAMDOM interfaces = lo,enp15s0 bind interfaces only = Yes dns forwarder = 8.8.8.8

Hello fellow R people, I can not figure out a pretty way to use strplit with vectors Imagine that I got the following data from someone with ID's representing several factors ID data A1-B1-t1 0 A1-B1-t2 1 A1-B2-t1 5 A1-B2-t2 10 A1-B10-t1 0 A1-B10-t2 1 A1-B20-t1 5 A1-B20-t2 10 ... I would like to turn this dataframe to station substation time data A1 B1 t1 0 A1 B1 t2 1 A1 B2 t1 5

Hi there, I would like to solve the following equation in R to estimate 'a'. I have the amp, d, x and y. amp*y^2 = 2*a*(1-a)*(-a*d+(1-a)*x)^2 test data: amp = 0.2370 y= 0.0233 d= 0.002 x= 0.091 Can anyone suggest how I can set this up? Thanks, Diviya [[alternative HTML version deleted]]

Hello there, I am using NLS for fitting a complex model to some data to estimate a couple of the missing parameters. The model is - y ~ (C+((log(1-r))*exp(-A*d)*(-1+r+exp(d*(A-B)))/(r*(-A*d+d*B+log(1-r))))) where A, B and C are unknown. In order to test the model, I generate data by setting values for all parameters and add some noise (C). A <- 20 B <- 500 r <- 0.6881 d <-