similar to: question about nlme...

Help R list servers, I have 500 external data sets for a simulation study that I would like to bring into R to analyze. They have the names data1.dat, data2.dat, ..., dataN.dat Is there a way to automatize the bringing in of these data sets in R using the read.table function within a looping cycle? For example... for (j in 1:N){ data_"j" = read.table("data"j".dat,

Respected All, I hope you are enjoying good health, I am tring to write a program in R but could not be very sucessful. My program draws random sample form bivariate normal distribution and then compute a variable PIJ. For certian samples some entries of variable PIJ is apearing as negative, which result in negative variance estimator. I want to introduce a loop in my program that verify the each

Douglas Bates wrote: > > The GLMM function in the lme4 package allows you to specify crossed > random effects within the random argument without the need for the > pdBlocked and pdIdent constructions. Simply ensure that your grouping > factors are defined in such a way that each distinct group has a > different level in the grouping factor (this is usually not a problem

Hi, I need help regarding estimating a linear model where restrictions are imposed on the coefficients. An example is as follows: Y_{t+2}=a1Y_{t+1} + a2 Y_t + b x_t + e_t restriction a1+ a2 =1 Is there a function or a package that can estimate the coefficient of a model like this? I want to estimate the coefficients rather than test them. Thank you for your help Ahmad Abu Hammour --------------

We don't currently support SVR4-style fd passing which involves I_SENDFD/I_RECVFD ioctls. I'm not sure if that might help some p platforms with privsep. Some platforms that use STREAMS ptys use a set-uid root helper to handle the privileged portions of pty creation. grantpt(3) has a bit more info. on Solaris 8 it's: /usr/lib/pt_chmod on HP-UX it's: /usr/lbin/chgpt Perhaps this

The ELF file header isn't always covered by a segment but still affects loading. I think everything else that effects loading/dynamic linking is always covered by a PT_LOAD segment. As evidence this is basically what --strip-sections in llvm-strip and eu-strip do and they produce perfectly runnable binaries. Having a hash of the actual memory map is interesting IMO. Build IDs can't really

Hi Toni, I'm sorry, I haven't had much time to work on this today. When I run sshd (from the patched snapshot) in a debugger, with a breakpoint early in setup_sia(), this is what I find after connecting with a client: (1) There are two sshd processes. One is running as root, and the other as the user I logged with using the client. The root process is the one in the debugger,

Hi, all, I'd like to find one or more changepoints in a dataset eg for data like 0.9 1 0.8 0.8 0.85 0.3 0.9 0.9 0.9 1.1 0.3 is unlikely a changepoint but for data like 0.9 1 0.8 0.8 0.85 0.3 0.2 0.2 0.1 0.2 0.3 0.3 is a change point I am wondering is there any function can do changepoint analysis? Ying Fan

Hello, I would like to evaluate the relationship between flows and phytoplankton abundance (or Chlorophyll a concentrations) using a changepoint analysis.? Specifically, I have two study questions: Study Question 1: Are there certain flow thresholds that result in spikes in phytoplankton abundance? Study Question 2: Are the duration of certain flows important for phytoplankton abundance (e.g.,

Dear All, I'm trying to do the estimation in a changepoint regression problem via R, but never found any suitable function which might help me to do this. Could someone give me a hand?on this matter? Thank you.

Hello Does anyone know how to get the p values for changepoints in the "changepoint" package? There is guidance in the CRAN vignette https://cran.r-project.org/web/packages/changepoint/changepoint.pdf but I can't make it work to get the p value out... Thanks Nick Wray [[alternative HTML version deleted]]

Hello. I have made many pngs in the past, but started having trouble with them yesterday. I just downloaded R 2.13.0 this morning, and was able to make some pngs earlier this morning. Now I am trying to remake them, but get the same error message I got yesterday: > u='in' > w=10.5 > h=8 > png('winbugs changepoint graphs ind hibs gt 400

Hello The other say I posted a question about altering confidence limits in the seqMK() function in the pheno package. I didn't get any bites but in any case I've been trying to use this function. It seems to me to give very odd results - sometimes identifying 8 or 9 changepoints in a relatively short vector of maybe 40 elements, and at other times not giving any changepoints at all,

I am trying to understand my population abundance data and am looking into analyses of change point to try and determine, at approximately what point do populations begin to change (either decline or increasing). Can anyone offer suggestions on ways to go about this? I have looked into bcp and strucchange packages but am not completely convinced that these are appropriate for my data. Here is

Hello All, I've written basic support for creating AIX native packages (ie installp/SMIT installable .bff's) that lives in contrib/aix. It is based on Ben Lindstrom's Solaris buildpkg.sh. It requires no special tools and all code was either written by me or derived from the openssh distribution itself. It is made available under the same terms as buildpkg.sh (ie public domain).

Dear friends. I would like to link an object file made in C. Dyn.load rather wants to link a whole library and I only wanted to link this small object file. I followed instructions given for S - by Clive Loader in his software for changepoint handling and it is his object file I made from http://cm.bell-labs.com/cm/ms/departments/sia/clive/smcp.shar Any dirty tricks or easy oversights ? Yours

Hello does anyone know how to set the confidence level within the seqMK() function in the pheno package. It seems to be set automatically at 0.05 and there seems to be neither an input function to set a different level, as there are with some changepoint functions, nor an output to give one the p-value of the changepoint identification which you get with functions ie func$p_value etc ... Thanks

Dear R-Help list, I have a nonlinear least squares problem, which involves a changepoint; at the beginning, the outcome y is constant, and after a delay, t0, y follows a biexponential decay. I log-transform the data, to stabilize the error variance. At time t < t0, my model is log(y_i)=log(exp(a0)+exp(b0)) at time t >= t0, the model is log(y_i)=log(exp(a0-a1*(t_i - t0))+exp(b0=b1*(t_i -

The 6th annual Modern Modeling Methods Conference will be held on May 24-25, 2016 at the University of Connecticut. This year's keynote speakers are Andrew Gelman and Bengt Muthen. May 23rd, Bengt Muthen is offering a full-day pre-conference workshop on Advances in Latent Variable Modeling in Mplus. May 26, Jeffrey R. Harring is offering a full-day post-conference workshop on Mixture Modeling

Hi, I am having trouble installing a previous version of a package, and multiple version of the package in R 2.12 or 2.15. I have download the tar.gz packages from the CRAN archive, and created an folder holding all the files. install.packages("changepoint_0.3.tar.gz", lib="~/Rtools/dir/of/pkgs/", repos=NULL, dependencies=FALSE,