Re: [R] bVar slot of lmer objects and standard errors

). Ulrich: Would you mind checking this either with a published example or a Monte Carlo and reporting the results to us? Viel Glück, spencer graves Ulrich Keller wrote: Hello, I'm sorry to resurrect this thread that I started almost two months ago

Re: [R] panel-VAR

Have you looked at the nlme package? If you provide more detail about your application (as suggested in the posting guide! www.R-project.org/posting-guide.html), you might receive more useful replies more quickly. hope this helps, spencer graves [EMAIL PROTECTED

Re: [R] OT Futility Analysis

and then comparing them with appropriate simulations -- summarized as confidence intervals, as you suggest. And I hope that someone else will enlighten us both if there are better options available. Best Wishes, spencer graves p.s. For any attorneys who may read

Re: [R] lme, nlsList, nlsList.selfStart

Have you tried making mydata2$Dose the same length as everything else? spencer graves Patrick Giraudoux wrote: Well, right, Dose was indeed in the global environment and not in the data.frame. Changing it with mydata2$Dose-100 # the real dose at the beginning

Re: [R] effect sizes in lme/ multi-level models

, and VarCorr. hope this helps. spencer graves Leo Gürtler wrote: Dear alltogether, I am searching for a way to determine effect size in multi-level models by using lme(). Coming from Psychology, for ordinary OLS there are measures (for meta-analysis, etc.) like CohensD

Re: [R] glmmPQL and random effects

. Espero que estos commentos le ayude. Spencer Graves [EMAIL PROTECTED] wrote: Hello R users, I am trying to run a model with a binary response variable (nesting success: 0 failure, 1 success) and 8 fixed terms. Nesting success was examined in 72 cases in 34 territories (TER

Re: [R] Lmer with weights

this helps. spencer graves ## I agree: The lmer weights argument seems not to have any effect. To check this, I modified the first example in the lmer documentation as follows: Sleep - sleepstudy Sleep$wts - 1:180 (fm1 - lmer(Reaction ~ Days

Re: [R] Baysian model averaging method

or more functions and try to walk through the code line by line until I thought I understood it well enough to modify it to do what I wanted. With luck, someone may reply to this email with hints or even code that might make that task easier. hope this helps, spencer graves p.s

Re: [R] Plotting contour filled.contour in one graph

, or are they still looking for volunteers to create the necessary help files, etc.? Best Wishes, spencer graves Gabor Grothendieck wrote: Could you walk us through, in detail, what that graph is showing? On 2/12/06, Michael Prager [EMAIL PROTECTED] wrote: Besides

Re: [R] Plotting contour filled.contour in one graph

and y must be in the upper right corner. Thanks for this. spencer graves Michael Prager wrote: GG Yes, gladly. It is an idealized example of the following data situation: There are two control or independent variables. They are represented here as x and y

Re: [R] Help with the ts function

can work through line by line, testing modifications, etc., as you go. This is my preferred way to learn about new capabilities in R. If you are not familiar with vignettes, see http://finzi.psych.upenn.edu/R/Rhelp02a/archive/67006.html;. hope this helps. spencer graves p.s

Re: [R] decomposed.ts class and method

plot.dendrogram* snip...? Thus, 'getAnywhere' should give you the plot method in this case. The class definition is outlined under 'Value' in the decompose help page. For an example, try str(m). hope this helps. spencer graves Erin Hodgess wrote: Dear R

Re: [R] marginal distribution wrt time of time series ?

increases, the number of observations to obtain reasonable convergence will increase. Standard goodness of fit test will NOT work with dependent series, but that's another issue. Perhaps someone else will provide further details. hope this helps. spencer graves

Re: [R] lme4: Error in getResponseFormula(form) : Form must be a two sided formula

tried this in nlme? If you absolutely need some feature of lme4, at least you could do this kind of preliminary work in nlme, then switch to lme4 (after quitting and restarting R to avoid potential conflicts between nlme and lme4). hope this helps. spencer graves

Re: [R] post-hoc comparisons following glmm

. hope this helps, spencer graves Michaël Coeurdassier wrote: Dear R community, I performed a generalized linear mixed model using glmmPQL (MASS library) to analyse my data i.e : y is the response with a poisson distribution, t and Trait are the independent variables which

Re: [R] Fixing AR coefficients in VAR model

, spencer graves Daniel Medina wrote: Dear Colleague, I would like to set a few AR coefficients (not order) to zero in the multivariate AR function (mAr.est; mAr library); however, the manual for this function does not provide this information. I would appreciate any suggestions along this line

Re: [R] Fixing AR coefficients in VAR model

. Best Wishes, spencer graves Paul Gilbert wrote: You can do this with dse. See require(dse1) ?fixConstants Paul Gilbert Spencer Graves wrote: I know of no existing functions in R that support fitting a multivariate autoregression while fixing some

Re: [R] Significance of degrees of freedom in nlme

if you make it easier for people to help you. hope this helps, spencer graves gabriela escati peñaloza wrote: Dear Dr. Bates, Thank you very much for your response. I had consulted the algorithm described in Pinheiro and Bates. However, what I don't understand (among other

Re: [R] how to use mle?

) returned 130 hits for me just now. hope this helps. spencer graves ronggui wrote: Y [,1] [,2] [,3] [1,]010 [2,]010 [3,]001 [4,]100 [5,]001 [6,]001 [7,]100 [8,]10

Re: [R] glm-logistic on discrete-time methods with individual and aggregated data

regression). Sehr interessant, nicht? hope this helps. spencer graves Camarda, Carlo Giovanni wrote: Dear R-Users, without going into details I tried to prepare a simple example to show you where I would need help. In particular I prepare two examples-template

Re: [R] Passing additional paramaters to nlsList(nlme) fit function

)'; the asterisk (*) says that this function is non-visible, which means that just typing its name won't get it. Then I might use debug to figure out what it's doing and what I want to change. hope this helps. spencer graves Dieter Menne wrote: Hello, nls-users

Re: [R] weights argument in the lmer function in lme4

can use this info as he sees fit. hope this helps. spencer graves sessionInfo() R version 2.2.1, 2005-12-20, i386-pc-mingw32 attached base packages: [1] methods stats graphics grDevices utils datasets [7] base other attached packages: lme4 lattice

Re: [R] Density estimation with monotonic constaints

useful information more quickly. spencer graves Debayan Datta wrote: Hi All, I have a sample x={x1,x2,..,xn} fom a distribution with density f. I wish to estimate the density. I know a priori that the density is monotonically decreasing. Is there a way to do this in R? Thanks

Re: [R] output hessian matrix in constrOptim

-guide.html. hope this helps, spencer graves liu wrote: Hi, Is there any way to get the hessian matrix from the constrOptim function without supplying gradient function? Thanks. - Bring words and photos together

Re: [R] non linear 3SLS with constraints

likelihood, or if you don't think transformations are appropriate, please submit another question to this listserve explaining as clearly and succinctly as you can why you think so. hope this helps. spencer graves Van Campenhout Bjorn wrote: hi, i am new here and wanted

Re: [R] fExtreme packages

, but they are not unfortunately well organized. hope this helps. spencer graves [EMAIL PROTECTED] wrote: Hello, I am a new user of R. I am trying to use the packages fBasics and fExtremes when i am running the examples I get few error. Could someone tell me what

Re: [R] Filtering the time series

management that it comes most highly recommended, and that the purchase of a copy will be an exceptionally good investment, in my judgment. Hope this helps. spencer graves anil kumar rohilla wrote: Hi List, I have a time series of 122 values, actualy it is a time series

Re: [R] How do I normalise a power spectral density analysis?

that posting guide will, I believe, tend to increase the speed and utility of response.) hope this helps. spencer graves Tom C Cameron wrote: Hi everyone Can anyone tell me how I normalise a power spectral density (PSD) plot of a periodical time-series. At present I get

Re: [R] Complex Matrix Exponentials.

question behind this, please let us know. hope this helps. spencer graves Andrew wrote: Hello, I was curious if there was a complex valued matrix exponential function available for R? I have some Laplace transforms of occupation times for a hidden Markov model

Re: [R] Question about Aggregate

. hope this helps, spencer graves Matthieu Cornec wrote: hello, Suppose you a monthly series you want to aggregate at a quaterly frequency with the start and the end of your series in the middle of the quarter. For example 2001M2 2001M3 2001M4 2001M5 2001M6 2001M7 12

Re: [R] fitting generalized linear models using glmmPQL

the posting guide! www.R-project.org/posting-guide.html, especialy the bit about providing a simple, self contained example. I sometimes solve problems like this in the course of trying to prepare a simple example. hope this helps. spencer graves [EMAIL PROTECTED] wrote

Re: [R] Linearize a Function

). Similarly, to compute log.returns from simple.returns, compute log1p(simple.returns) = log(1+simple.returns) [making the obvious conversions between percentages and proportions]. Or am I missing something? hope this helps, spencer graves Gottfried Gruber wrote: hi, i

Re: [R] draft of Comment on UCLA tech report

from this one comment, what you wrote sounds sensible to me. Best Wishes, Spencer Graves Patrick Burns wrote: You may recall that there was a discussion of a technical report from the statistical consulting group at UCLA. I have a draft of a comment on that report, which

Re: [R] weighted likelihood for lme

recommeded it to more people than any other single work over the past five years. hope this helps. spencer graves Marco Geraci wrote: Dear R users, I'm trying to fit a simple random intercept model with a fixed intercept. Suppose I want to assign a weight w_i to the i

Re: [R] nlme in R v.2.2.1 and S-Plus v. 7.0

. hope this helps, spencer graves Paolo Ghisletta wrote: Dear R-Users, I am comparing the nlme package in S-Plus (v. 7.0) and R (v. 2.2.1, nlme package version 3.1-68.1; the lattice, Matrix, and lme4 have also just been updated today, Jan. 23, 2006) on a PC (2.40 GHz Pentium 4

[R] D(dnorm...)?

, this should be something like ((mean-x)/sd^2)*dnorm(...). Thanks for your help. Spencer Graves __ R-help@stat.math.ethz.ch mailing list https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read the posting guide! http://www.R-project.org

Re: [R] D(dnorm...)?

With dnorm, I got 0!. Whe didn't I get an error message? Thanks for the reply. Spencer Graves Berton Gunter wrote: dnorm() is an internal function, so I don't see how D (or deriv) can do anything with it symbolically. Am I missing something? -- Bert

Re: [R] D(dnorm...)?

Hi, Prof. Ripley: Thanks for the explanation. If I had read your book more carefully, I would not have needed this email exchange. Thanks again, Spencer Graves Prof Brian Ripley wrote: On Wed, 25 Jan 2006, Spencer Graves wrote: Hi, Bert: I think I

Re: [R] regression with nestedness

Have you considered lme in library(nlme)? The companion book Pinheiro and Bates (2000) Mixed-Effects Models in S and S-Plus (Springer) is my favorite reference for this kind of thing. From what I understand of your question, you should be able to find excellent answers in this

Re: [R] D(dnorm...)?

like Excel, but ... . Spencer Graves [EMAIL PROTECTED] wrote: Yes Bert, this time you are missing something (unusually) ... As Brian Ripley pointed out 'dnorm' is in the derivative table, *but* only as a function of one variable. So if you want to find the derivative

Re: [R] Within-Subjects ANOVA comparisons of individual means

might have a not-quite-fully debugged version someplace that could help you. Best Wishes, Spencer Graves Steffen Katzner wrote: I am having problems with comparing individual means in a within-subjects ANOVA. From my understanding, TukeyHSD is not appropriate

Re: [R] R-help Digest, Vol 35, Issue 24

standard for academic excellence in comportment and the use of language? Sincerely, Spencer Graves [EMAIL PROTECTED] wrote: Dear Prof Ripley, First of all, unless you are an english professor, then I do not think you have any business policing language. I'm still very much

Re: [R] mutlivariate normal and t distributions

it may not be clear which is best. However, if you and Torsten don't have time energy for that, I would still encourage you to submit your package to CRAN. Best Wishes, Spencer Graves Adelchi Azzalini wrote: Dear R-help list members, I have created a package 'mnormt

Re: [R] In which application areas is R used?

an entry. What do you think? spencer graves John Maindonald wrote: In this context extensive might be use of R in at least maybe 2% or 5% of the published analyses in the area, enough to make waves and stir awareness. The immediate subtext is the demand of a book

Re: [R] lme model specification

wishes, spencer graves # Bill Simpson wrote: Thanks very much Spencer for your helpful reply. On Fri, 2006-01-20 at 07:38 -0800, Spencer Graves wrote: Does each subject get only one LED per session or all 4 LEDs? I simplified a bit. Each subject gets pairs

Re: [R] Influence measure + lme ?

to write an 'influence' method for an 'lme' object. However, that is far beyond what you asked and what I have time now to disucss. Hope this helps. spencer graves Chia, Yen Lin wrote: Hi all, Does lme has function to compute the cook's distance or influence

Re: [R] mcmcsamp() in lmer

useful. I'd be much obliged to anyone who could help understand how I could diagnose this issue. Best Wishes, Spencer Graves Andrew Gelman wrote: I am working with lmer() in the latest release of Matrix, doing various things including writing a function

Re: [R] mcmcsamp() in lmer

Hi, Doug: Thanks for all your hard work and creativity with mixed-effects models -- and for your quick reply to my question about how to use 'showMethods' to view the source. Best Wishes, spencer graves Douglas Bates wrote: On 1/21/06, Spencer Graves [EMAIL

Re: [R] numericDeriv() giving a vector when multiple variables input

. - A%*%G) (AG.%*%1:4) This seemed to work for me and give me the correct answer. If this is not enough, RSiteSearch(numericDeriv) gave me 115 hits; many of these are not relevant to your question, but I would expect that some might be. hope this helps. spencer graves

Re: [R] lme model specification

) spencer graves Bill Simpson wrote: I have been asked to analyse the results of (what is to me) a very complicated experiment. The dependent measure is the estimated distance, which is measured as a function of the actual distance. There are also several other IVs. The plot of log

Re: [R] Help with mixed effects models

a replicable example, as suggested in the posting guide (www.R-project.org/posting-guide.html), I might have been able to help more -- or someone else might have answered your question sooner. hope this helps, spencer graves Pryseley Assam wrote: Dear R-users I have

Re: [R] glmmPQL: Na/NaN/Inf in foreign function call

, this person must also have the time and interest to reply. Don't gamble on a long shot. Play the odds. PLEASE do read the posting guide! www.R-project.org/posting-guide.html. spencer graves David Reitter wrote: I'm using glmmPQL, and I still have a few problems

Re: [R] Homogenous groups building - Randomisation

this and would still like help from this listserve, please think very carefully about what it would take for someone to produce a useful response to your question. Je voudrais vous aidez, mais ce demand plus que j'ai. Hope this helps, spencer graves [EMAIL PROTECTED] wrote

Re: [R] singular convergence(7)?

darts blindly, hoping one will find the bulls-eye of a target you only hope is there. hope this helps. spencer graves Chia, Yen Lin wrote: Hi all, I just wonder what singular convergence means. Thanks. Yen Lin Error in lme.formula(Data ~ 1, random = ~1

Re: [R] fitted values from lmer (lme4 0.98)

Hi, Doug: I think it would help me to have an option like you just suggested -- with your comment below included in the help file for that option. Thanks for all your hard work and creativity in this. Best Wishes, spencer graves Douglas Bates wrote

Re: [R] find mean of a list of timeseries

()-start.time) [1] 0.00 0.00 0.07 NA NA This seems to run slightly faster on this miniscule data. However the answer is transposed, and I don't want to take the time to understand and fix that problem. hope this helps. spencer graves tom wright wrote: Can

Re: [R] first derivative of a time series

to solve, for which you think the first derivative of an original time series might be useful. spencer graves p.s. The fact that you included an almost reprodicible example was critical in permitting me to reply. Without that, the most constructive comment I might have been able

Re: [R] lmer with nested/nonnested groupings?

. spencer graves Andrew Gelman wrote: I'm trying to figure out how to use lmer to fit models with factors that have some nesting and some non-nested groupings. For example, in this paper: http://www.stat.columbia.edu/~gelman/research/published/parkgelmanbafumi.pdf we have a logistic

Re: [R] trouble with extraction/interpretation of variance structure para meters from a model built using gnls and varConstPower

by line, until I found the following: (val - as.vector(model3$modelStruct$varStruct)) Variance function structure of class varConstPower representing const power 0.6551298 0.8913665 Answer the questions? spencer graves Rand, Hugh wrote: I have been using gnls

Re: [R] glmmPQL / system is computationally singular

can achieve perfect separation, software of this type might generate error messages similar to what you describe. hope this helps. spencer graves David Reitter wrote: Hi, I'm having trouble with glmmPQL from the MASS package. I'm trying to fit a model with a binary

Re: [R] Homogenic groups generation - Randomisation

of getting help from someone who understands the concepts but not your particular jargon. J'espere que ceci vous aidera. spencer graves [EMAIL PROTECTED] wrote: Dear R-users, We expect to create N homogenic groups of n features from an experimentation including N*n mesures

Re: [R] glmmPQL and variance structure

weights.lme to Venables and Ripley's glmmPQL and uses that in place of 'quote(varFixed(~invwt))' when provided; see below. hope this helps. spencer graves glmmPQL. - function (fixed, random, family, data, correlation, weights, weights.lme, control, niter = 10, verbose = TRUE

Re: [R] panel data unit root tests

] } UnitRootSims AR1 p.value [1,] NA 8.277712e-11 [2,] NA 1.174823e-08 # Why are the p-values so small? Shouldn't they be insignificant? Thanks, Spencer Graves If replies to this post will no longer be useful for you

Re: [R] glmmPQL and variance structure

hope this error doesn't seriously inconvenience you or anyone else who might read these comments. spencer graves Thanks for providing a partially reproducible example. I believe the error message you cite came from lme. I say

Re: [R] CRAN versions of lme4/Matrix don't appear to work with R 2.1.1

This can almost eliminate the possibilities for typographical errors in specifying the version name. The posting guide (www.R-project.org/posting-guide.html) also asks posters to mention mention the platform (Windows2000, Linux, MacOS X). Best Wishes, spencer graves

Re: [R] panel data unit root tests

, your code includes print(levinlin(ws, year, id, lags = 3)). What is the levinlin function? RSiteSearch(levinlin) produced zero hits. hope this helps. spencer graves jukka ruohonen wrote: When finally got some time to do some coding, I started and stopped right after

Re: [R] CRAN versions of lme4/Matrix don't appear to work with R 2.1.1

.) Doug Bates has worked very hard to create this, and I for one would not want to ask him to take the time to make it backward compatible for those who for whatever reason haven't yet upgraded to R 2.2.1, especially since for most of us it is fairly easy to upgrade. Spencer Graves White

Re: [R] CRAN versions of lme4/Matrix don't appear to work with R 2.1.1

It worked fine for me using http://cran.fhcrc.org/ Fred Hutchinson Cancer Research Center, Seattle, WA spencer graves Peter Dalgaard wrote: White, Charles E WRAIR-Wash DC [EMAIL PROTECTED] writes: I blew away my existing R directory structure, reinstalled R 2.2.1 from

Re: [R] jointprior in deal package

. hope this helps. spencer graves PLEASE do read the posting guide! http://www.R-project.org/posting-guide.html Tim Smits wrote: Dear all, I recently started using the deal package for learning Bayesian networks. When using the jointprior function on a particular dataset

Re: [R] CRAN versions of lme4/Matrix don't appear to work with R 2.1.1

I don't know, but I had a similar problem a few weeks ago with the one or both of the California sites (I don't remember which now). The problems disappeared after I switched to http://cran.fhcrc.org/; (Fred Hutchinson Cancer Research Center, Seattle, WA). spencer graves

Re: [R] information

esta guia tipicamente receiben contestaciones mas rapido y mas utiles. spencer graves angel toledo wrote: Hi. My name is Angel, I am Mexican, and I write by the following thing: I am in search of commands or options in R that can be used in regional economics. Specially I am

Re: [R] glmmPQL error message (was 'data order affects glmmPQL')

how I might change the code to get it past this error and on to something close to convergence. Hope this helps. spencer graves Jack Tanner wrote: From: Spencer Graves The correlation between the predictions from your two model fits is 0.95. This suggests to me

Re: [R] Selecting significant peaks in periodograms

evidence suggests that posts more consistent with that guide tend to get more useful replies quicker. spencer graves Pete Cap wrote: Greetings all, I am using Fourier analysis to search for periodicities in IP network traffic by generating periodograms and then visually

Re: [R] data order affects glmmPQL

. spencer graves p.s. To got this correlation, I modified the latter portion of your script as follows: fit.foo - fit.model(foo, poisson) o.bar - order(id, score, test, coder) bar - foo[o.bar,] # Reorder data frame fit.bar - fit.model(bar, poisson) predict.foo - predict(fit.foo) predict.bar

Re: [R] lmer error message

, Spencer Graves # R version 2.2.0, 2005-10-06, i386-pc-mingw32 attached base packages: [1] methods stats graphics grDevices utils datasets [7] base other attached packages: lme4 latticeMatrix 0.98-1 0.12-11 0.99-4 Douglas Bates wrote: Also

Re: [R] lmer error message

every case that I've tried suggestion 4, I've been able to figure out how to overcome that particular difficulty. In addition, I've often learned useful things about R that I didn't know befor. hope this helps. spencer graves Abderrahim Oulhaj wrote: Dear All, I have

Re: [R] Use Of makeARIMA

. spencer graves p.s. If you'd like more help from this group, please submit another question. Before you do, however, I suggest you first read the posting guide! www.R-project.org/posting-guide.html. Anecdotal evidence suggests that posts more consistent with that guide are more

Re: [R] Glimmix and glm

this group, PLEASE do read the posting guide! www.R-project.org/posting-guide.html. Anecdotal evidence suggests that posts that conform more closely to the suggestions there tend to get quicker more useful replies. Best Wishes, Spencer Graves [EMAIL PROTECTED] wrote: Hello

Re: [R] unexpected false convergence

the code fo lme.formula to force it to call optim, but this generated a different error. I am therefore copying Professors Bates Ripley in case one of them might want to look at this. hope this helps. spencer graves Jack Tanner wrote: I've come into some code that produces

Re: [R] glmmPQL and variance structure

commits a fatal error, you will be evicted from glmmPQL and will have to start over. This is the quickest way I know to understand and debug R code. hope this helps. spencer graves Patrick Giraudoux wrote: Dear listers, glmmPQL (package MASS) is given to work by repeated

Re: [R] Different ARCH results in R and Eviews using garch from tseries

.) Or write a function to compute the negative of the log(likelihood) and use optim to minimize it, with hessian = TRUE to get the observed information, whose inverse is the variance for the Wald approximation. hope this helps. spencer graves Constantine Tsardounis wrote

Re: [R] Parameter Constraints in nls.lm()

or returned NAs. That was a few years ago, and the functions may be better today, but that's what I've done. hope this helps. spencer graves Kilian Plank wrote: Hello, Is there a possibility of setting parameter constraints in nls.lm() ? The documentation does not say

Re: [R] bVar slot of lmer objects and standard errors

. If not, maybe he can find another way to check the answer, e.g., via Monte Carlo. Best Wishes, Spencer Graves Douglas Bates wrote: On 12/29/05, Doran, Harold [EMAIL PROTECTED] wrote: Uli: The graphic in the paper, sometimes called a catepillar plot, must be created with some

Re: [R] bVar slot of lmer objects and standard errors

) fm1.s - coef(fm1OrthF.)$Subject fm1.s.var - [EMAIL PROTECTED] fm1.s0.s - sqrt(fm1.s.var[1,1,]) fm1.s0.a - sqrt(fm1.s.var[2,2,]) fm1.s[,1]+outer(fm1.s0.s, c(-2,0,2)) fm1.s[,2]+outer(fm1.s0.a, c(-2,0,2)) hope this helps. Viel Glueck. spencer graves Ulrich Keller

Re: [R] Testing a linear hypothesis after maximum likelihood

that follow more closely the suggestions in that guide tend to get more useful replies quicker. hope this helps. spencer graves Peter Muhlberger wrote: I'd like to be able to test linear hypotheses after setting up and running a model using optim or perhaps nlm. One hypothesis I

Re: [R] Testing a linear hypothesis after maximum likelihood

. hope this helps. spencer graves On 12/29/05 7:04 AM, Spencer Graves [EMAIL PROTECTED] wrote: Why can't you use a likelihood ratio? I would write two slightly different functions, the second of which would use the linear constraint to eliminate one

Re: [R] Testing a linear hypothesis after maximum likelihood

, for testing, when I have the time, I use likelihood ratio procedures. spencer graves Peter Muhlberger wrote: On 12/29/05 1:35 PM, Spencer Graves [EMAIL PROTECTED] wrote: I think the question was appropriate for this list. If you want to do a Wald test, you might consider asking optim

Re: [R] parameterization of factor in R

hope it will help you solve, etc. Spencer Graves Prof Brian Ripley wrote: ?contr.sum And Chapter 6 of MASS would have explained all this to you, so if you have further questions please consult it. On Tue, 27 Dec 2005, Oana Mocila wrote: I encountered this problem

Re: [R] how to put constraint in R?

- Manuals). If you would like more help from this group, I believe you will increase your chances of getting quickly the information you seek in you PLEASE do read the posting guide! www.R-project.org/posting-guide.html. hope this helps. spencer graves Oana Mocila wrote

Re: [R] glmmADMB: Generalized Linear Mixed Models using AD Model Builder

with the singular hessian. Without that, I have to program optim or something else separately to get the answers and the hessian in order to do my own diagnosis -- if I know enough to do that. Just my 0.02 Euros. spencer graves Roel de Jong wrote: Of course it is generally

Re: [R] nlme problems

Are you familiar with Pinheiro and Bates (2000) Mixed-Effects Models in S and S-Plus (Springer)? I suspect that book, especially the latter half, might contain the information you seek. spencer graves p.s. PLEASE do read the posting guide! www.R-project.org/posting

Re: [R] given a mid-month date, get the month-end date

. Instead, I use the following: Stangle(zoo.vignette$file) This writes the R commands to file zoo.R in the working directory. From there, I then use menu File - Open - “zoo.R”. hope this helps. Spencer Graves p.s. In case it's not completely obvious from the above, I

Re: [R] autocorrelation test

submit another question. Before you do, however, I suggest you read the posting guide! www.R-project.org/posting-guide.html. Anecdotal evidence suggests that posts more consistent with this guide tend to receive quicker, more useful replies. hope this helps. spencer graves

Re: [R] partially linear models

. Best Wishes, spencer graves Elizabeth Lawson wrote: Hey, I am estiamting a partially linear model y=X\beta+f(\theta) where the f(\theta) is estiamted using wavelets. Has anyone heard of methods to test if the betas are significant or to address model fit

Re: [R] ANCOVA Post-hoc test

What search terms did you use? Have you considered the multcomp and multtest packages? spencer graves p.s. If you'd like more help from this list, I suggest you read the posting guide! www.R-project.org/posting-guide.html. Anecdotal evidence suggests that posts more closely

Re: [R] suggestions for nls error: false convergence

) Nonlinear Regression Analysis and Its Applications (Wiley, esp. ch. 6 on Graphical Summaries) hope this helps. spencer graves Rajarshi Guha wrote: Hi, I'm trying to fit some data using a logistic function defined as y ~ a * (1+m*exp(-x/tau)) / (1+n*exp(-x/tau

Re: [R] suggestions for nls error: false convergence

comments on this example. Spencer Graves Christian Ritz wrote: Hi. An alternative is to use the package 'drc' on CRAN to fit your data! x - 1:100 y - c(0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0, 0,0,0,0,0,1,1,1,2,2,2,2,2,3,4,4,4,5, 5,5,5,6,6,6,6,6,8,8,9,9,10,13,14,16,19,21

Re: [R] residual and null deviance of an lme object with correlation structure

this helps, spencer graves Corey Bradshaw wrote: Hello, I am attempting to calculate the residual and null deviance of an lme object that includes a corAR1 correlation structure. I tried deviance(lme.object) and it only returned NULL. Can anyone help? Thank you. Corey Bradshaw

Re: [R] Kalman Filter Forecast using 'SSPIR'

guide! www.R-project.org/posting-guide.html. Anecdotal evidence suggests that posts more consistent with that guide tend to get quicker, more useful replies. hope this helps. spencer graves -- Sumanta Basak wrote: Dear R Users, I am new

Re: [R] help with multivariate analysis

consistent with that guide tend to get more useful replies quicker. hope this helps. spencer graves Matteo Vidali wrote: dear R users, I need some help for multivariate analysis. I have 2 anaesthetic treatment groups (20 patients/group) where I register heart frequency

<    1   2   3   4   5   6   7   8   9   10   >