similar to: nlme starting values are not the correct length

Windows 2000 . 5.00.2195 with Service Pack 1. R 1.5.1 Output from my split-split plot aov "alerted" me that I have done something wrong. I designed an experiment with all combinations of all levels of each treatment, but lost a little data (3 out of 192 plots). With the following data, I run the following model: > collim[c(1:6,187:192),c(1,3:6,9)] plot Litter Fert

Hi, I've data from an incomplete fatorial design. One level of a factor doesn't has the levels of the other. When I use lm(), the summary() return NA for that non estimable parameters. Ok, I understant it. But I use contrast::contrast(), gmodels::estimable(), multcomp::glht() and all these fail when model has NA estimates. This is becouse vcov() and coef() has different dimensions. Is

Dear R sages, I used the function rbind to combine a matrix (M) and a vector (Fert) to get a new matrix (A). This was fine. The issue is however, that the new matrix A has as its row names the name of the vector Fert, even though I set teh new vector A to have dimnames=NULL. See short code below fyi. *Fert<-c(0,1,5) * *M <- matrix(0, 2, 3) diag(M) <- c(0.3,0.5) * *A<-

Windows 2000 (v.5.00.2195), R 1.5.1 I have an lme object, fm0, which I examine with anova() and intervals(). The anova output indicates one of the interaction terms is significant, but the intervals output shows that the single parameter for that term includes 0.0 in its 95% CI. I believe that the anova is a conditional (sequential) test; is intervals marginal or approximate? Which should I

Hi R users, I sent a message yesterday about NA in model estimates ( http://r.789695.n4.nabble.com/How-set-lm-to-don-t-return-NA-in-summary-td3722587.html). If I use aov() instead of lm() I get no NA in model estimates and I use gmodels::estimable() without problems. Ok! Now I'm performing a lot of contrasts and I need correcting for multiplicity. So, I can use multcomp::glht() for this.

Hi R-Helpers: I?m trying to use the function GLMM from lme4 package, (R-1.8.1, Windows 98),and I get the following error: > pd5 = GLMM(nplant~sitio+ + fert+ + remo+ + sitio:fert+ + remo:sitio+ + remo:fert+ + remo:fert:sitio + data=datos, + family=binomial, +

Using R v. 1.7.0 on Windows 2000 I would like to plot the fitted values of a model as a function of a continuous covariate, augmented with data (e.g., augPred) grouping by combinations of fixed effects. I have not been able to use augPred effectively, and am wondering if it does not handle unbalanced data (3 out of 192 missing). I include below the model and an xyplot that almost does the

Hello, I am a Masters student and I am working on my thesis modelling smallholder farms using a program in R. I have modified the original code and I am having some issues with lists that I cannot figure out. Originally, I had list file defining lists such as: Param, Crop1, Crop1, Soil, etc. (ex. Param <- list() ). Their subsets were listed as Crop1$CContent for example and there was quite

Hi I have a file saved in R, named agrexp.Rdata, shown below > agrdata fert yield 1 25 84 2 50 80 3 75 90 4 100 154 5 125 148 If I double clicked on this file, the data is displayed without problem. However if I tried to import using: > agrdata<-read.table("agrexp.Rdata") or >

Dear all I'm not sure if I did the right analysis for my specific split splot design. We are studying biomass increase with different CO2 concentrations with four different functional plant groups (e.g. grasses, herbs, broad-leafed trees and conifers). Of each functional plant group we have four species. The design is orthogonal. The design is: Blocks: 2 (climate chambers, called

Hi Rneters, I'm running a GLM model with a full factorial design in blocks and binomial error distribution. I would like to have the F values for this model but I got a message that "using F test with a binomial family is inappropriate in: anova.glm(model, test = "F")". Should I not report F statistics on this kind of analysis? I would appreciate any comment on this.

Hi, I'm trying to fit a logit model to a set of data that was collected under a split plot scheme. The structure of the data is Whole plot factors: Watering Frequency (2 levels: Hi/Lo) and Fertilizer type (3 factors A/B/C) Subplot factor: slope type (2 factors up/down) Response: Proportion of infected leaves(Infected leaves/Total leaves) of the plant (2 plants recorded from each plot)

Thanks a lot Dona and jg for your inputs. I'll try to find some way to do this from Dialplan or AMI and let you guys know soon. Please share if you have some more ideas. Regards, Rajib Date: Tue, 11 Jun 2013 18:34:46 +0200 From: jg <webaccounts at jgoettgens.de> Subject: Re: [asterisk-users] announcement to be played for attended transfer call To: Asterisk Users Mailing List -

Unsubscribe Elvin G. Nodalo -----Original Message----- From: asterisk-users-request at lists.digium.com Sent: 7/10/2013 1:00 AM To: asterisk-users at lists.digium.com Subject: asterisk-users Digest, Vol 108, Issue 14 Send asterisk-users mailing list submissions to asterisk-users at lists.digium.com To subscribe or unsubscribe via the World Wide Web, visit

Unsubscribe Elvin G. Nodalo -----Original Message----- From: asterisk-users-request at lists.digium.com Sent: 7/10/2013 1:00 AM To: asterisk-users at lists.digium.com Subject: asterisk-users Digest, Vol 108, Issue 14 Send asterisk-users mailing list submissions to asterisk-users at lists.digium.com To subscribe or unsubscribe via the World Wide Web, visit

Hi, It is a shame I feel obliged to write this email. I sent a single email with a possible bug in it, asking if it was working in a newer release. I got 3 spam emails from mindrot.org. Two asking me to confirm, and then when I was forced to confirm I got another auto-responder spam. If I had wanted the 3rd spam I would have set "return receipt". I got this (below) auto-responder

"sip show channels" shows some info about active sip channels, the current codec included. What does it say? jg" jg, sip show channels reports the Format as being ulaw for 17 active calls. Holds - no Peer User/ANR Call ID Format Hold Last Message Expiry Peer xxxxxxxxxx kbrown xxxxxxxx (ulaw) No Rx:

hi, i've always been using mount -t smbfs and recently "upgraded" to mount.cifs because i have a very aggressive application which always killed the smb connection (getting timeouts, ...). my setup: * some passwordless public/guest shares * one admin share over a x-over link to be able to use the admin share i had to switch my server from "security = share" to

> On Thursday 05 Feb 2015, jg wrote: >> Calling from ServerB to ServerA works, but not vice versa. The only odd >> thing that appears to me is the different perceived port on ServerA. >> >> Does someone have an idea at what to look in detail? > > Look in /etc/asterisk/iax.conf in the first instance. > Basically I used the example from the Asterisk book

2015-02-09 14:36 GMT-06:00 jg <webaccounts173 at jgoettgens.de>: > Hi! > > Sometimes IAX peers are not reachable and with "iax2 set debug on" I get > something like this > > Tx-Frame Retry[ No] -- OSeqno: 000 ISeqno: 001 Type: IAX Subclass: PONG > Timestamp: 00014ms SCall: 00001 DCall: 01200 79.233.155.174:49153 > Rx-Frame Retry[ No] -- OSeqno: 001