similar to: Memory explosion, plotting nmle grouped data object

Dear R users, I'm applying a correlation structure in a mixed model (nmle function) to control for spatial correlation between land parcels that are adjacent to each other. I generated X,Y coordinates in ArcGIS for each land parcel and used them in the correlation form like this: test.exp<-corExp(1, form = ~ X + Y) test.exp<- Initialize(test.exp,dataset) However, the correlation

Following the Indomethicin example in Pinheiro & Bates, chapter 6, page 277 etc, coming to the following comand: fm2Indom.nlme <- update( fm1Indom.nlme, random = pdDiag(A1 + lrc1 + A2 ~ 1) ) debugging nlme gives the following output: Browse[1]> n debug: modelResid <- ~eval(model, data.frame(data, getParsNlme(plist, fmap, rmapRel, bmap, groups, beta, bvec, b, level,

Hello. I have 16 subjects with 1-4 obs per subject. I am using the package "nlme" to fit a simple random effects (variance components model) with 3 parameters: overall mean (fixed effect), between subject variance (random) and within subject variance (random). I need a 3x3 variance-covariance matrix that includes all 3 parameters in order to compute the variance of a specific

Hi, All: Does 'nlme' require the nonlinear function to be differentiable? I'm fitting structural change models to many related time series, and I'd like to improve the estimates of the change point times through some type of pooling. Unfortunately, I've so far been unable to get 'nlme' to work for me. The following toy example is the closest I've come

A few days ago I posted a question to this discussion group concerning to origin of an error message < subscript out of bounds > while using the nonlinear mixed model (nlme) function in R with a self-starting function. Thanks for those who responded. This posting is to explain what (I think) it causing the error. Pinheiro & Bates (2000, pages 342-347) describe how to construct a

Hi, I am not sure this is a bug but I can repeat it, The functions and data are below. I know this is nasty data, and it is very questionable whether a 4pl model is appropriate, but it is data fed to an automated tool and I would have hoped for an error. Does this repeat for anyone else? My details: > version _ platform i686-pc-linux-gnu

Hi, So, we currently have gc.result.int, gc.result.float. gc.result.ptr, gc.relocate, and gc.statepoint. gc.statepoint's signature is fine with a iPTRAny as the first argument. gc.result is in trouble, because none of the signatures admit even a simple array of integers, and there's no aAny. And certainly no vectors. So we can get a gc.result.vector to add to this mess, and admit [1] to

Philip Reames wrote: >> 1. Introduce aAny. > > Having a generic any type seems fine. I assume you'd create something like > an llvm_any_type in Intrinsics.td? That's precisely what ifavpAny is about: integer, float, array, vector, pointer Any. aAny is meant for array-Any, and I wonder why so few people care about arrays. I'll go ahead with Any and llvm_any_type.

Hi, Am I hitting some limit in match? Consider the following example: > tst<-seq(100,125,by=.2)%in%seq(0,800,by=.1) > sum(tst) [1] 76 > seq(100,125,by=.2) [1] 100.0 100.2 100.4 100.6 100.8 101.0 101.2 101.4 101.6 101.8 102.0 102.2 [13] 102.4 102.6 102.8 103.0 103.2 103.4 103.6 103.8 104.0 104.2 104.4 104.6 [25] 104.8 105.0 105.2 105.4 105.6 105.8 106.0 106.2 106.4 106.6 106.8

Hi, Am I hitting some limit in match? Consider the following example: > tst<-seq(100,125,by=.2)%in%seq(0,800,by=.1) > sum(tst) [1] 76 > seq(100,125,by=.2) [1] 100.0 100.2 100.4 100.6 100.8 101.0 101.2 101.4 101.6 101.8 102.0 102.2 [13] 102.4 102.6 102.8 103.0 103.2 103.4 103.6 103.8 104.0 104.2 104.4 104.6 [25] 104.8 105.0 105.2 105.4 105.6 105.8 106.0 106.2 106.4 106.6 106.8

Hi, I have a many-to-many relationship between two tables, questions and answers. My MySQL schema for the link table creates a primary key from the two fields to prevent duplicate records. When I add the same answer to the same question more than once, rails attempts to create a duplicate record and explodes with a MySQL error: MysqlError: Duplicate entry ''3-3'' for key 1:

Hi, I would also be in favor of a stronger stance on licenses. In industry, where we can really get in big trouble for violating a license, we often maintain internal repositories, or need to be careful about filtering what is used from CRAN. I think that is should really be a requirement the package authors commit to stating what the restrictions are on their packages. Nicholas On 10

Following up to some extent on Friday's discussion regarding the 'validation' of R, could I ask the list group's opinion on possible advantages of R over Splus from a pharma/devices perspective? I wish to exclude the obvious price difference, which doesn’t seem to carry as much weight as I would have thought. Besides, I have noticed many former Splus users gravitating towards R,

I use asterisk 1.2.5 and h323 that comes with addons 1.2.1. Incoming call comes on h323 trunk. Person A (local SIP phone) receives call and tries to make attendant transfer to person B (local SIP phone). They speak. Then A hangs up. Call form h323 trunk doesn't get to person B. This is what I get on CLI. -- Attempting native bridge of OOH323/xxx.xxx.xxx.xxx-5381 and SIP/307-5663 --

Hi all, I have a very simple question that I can't seem to find the answer to. How do I extract rows that meet a certain criteria from a data frame and group them into a new data frame? For example, if I want to make a new data frame that only includes rows of data for which the p values (given by one of the columns in the data frame) are less than a certain value, how do I do this? It seems

Dear list Having just started to use the Windows version, I am very impressed with it's package handling as well as the gui. So I tried to see what was due for update and packages such as Hmisc, Matrix and others came up. But when I had updated them - which took a few goes as something hung between here and Bristol - I noticed that the default packages such as nmle, MASS had disappeared. I

If I have an IAX client (Firefly softphone in this example), and the client is not registered at the moment because they are not connected to the network and someone dial that extension, they get the user's "I'm on the phone at the moment" message vs. the "I'm unavailable" message. Is this by design? Here's the extension in question's dialplan:

I am running lxc via libvirt1.1.1, the cpu usage is heavy of process libvirt_lxc. PID USER PR NI VIRT RES SHR S %CPU %MEM TIME+ COMMAND 2384 root 20 0 10480 2048 1200 R 99.6 0.1 1291:39 libvirt_lxc

I'm trying to set up a mixed model to solve using lme. It will have 3 fixed effects, two random effects and two interaction terms. I've been reading Pinheiro's and Bates's book on the nmle library, but find the part about display functions to be unclear. When creating a groupedData object from a data.frame, you need to enter a function of the form: response ~primary|grouping

Hi guys! I have a problem regarding a binary logistic hierarchical model I am trying to use. The model contains various covariates that depend on the location the response was measured at but do not depend on time (year). I also have a spatial covariate that depends both on location and time. I have been trying to use the lme4 pack but the package only allows me to model variance covariance