Re: [R] (Statistics question) - Nonlinear regression and simultaneous equation

question, please provide more detail of your specific application, preferably including commented, minimal, self-contained, reproducible code, as requested at the end of every email forwarded by r-help. Hope this helps. Spencer Graves [EMAIL PROTECTED] wrote: Hi,I have

Re: [R] Empirical copula in R

=fgac). Hope this helps. Spencer Graves GWeiss wrote: Hi, I would like to implement the empirical copula in R, does anyone know if it is included in a package? I know it is not in the Copula package. This one only includes a gof-test based on the empirical copula process

Re: [R] how to solve a min problem

: I'm confused. Hope this helps. Spencer Graves domenico pestalozzi wrote: I know it's possible to solve max e min problems by using these functions: nlm, optimize, optim but I don't know how to use them (...if possible...) to solve this problem. I have a personal function

Re: [R] unequal variance assumption for lme (mixed effect model)

capabilities in the 'nlme' package is the best tool I know for that purpose. If someone thinks there are better tools available for estimating heterscedasticity, I hope s/he will enlighten us both. Hope this helps. Spencer Graves shirley zhang wrote: Thanks for Spencer

Re: [R] unequal variance assumption for lme (mixed effect model)

. Spencer Graves Thanks, Shirley On 6/27/07, Simon Blomberg [EMAIL PROTECTED] wrote: The default settings for lme do assume equal variances within groups. You can change that by using the various varClasses. see ?varClasses. A simple example would be to allow unequal variances across

Re: [R] Dominant eigenvector displayed as third (Marco Visser)

answer to your question. Hope this helps. Spencer Graves Marco Visser wrote: Dear R users Experts, This is just a curiousity, I was wondering why the dominant eigenvetor and eigenvalue of the following matrix is given as the third. I guess this could complicate automatic

[R] R.matlab questions

, Spencer Graves Using Matlab 7.3.0 (R2006b) under Windows XP Pro. sessionInfo() R version 2.5.0 (2007-04-23) i386-pc-mingw32 locale: LC_COLLATE=English_United States.1252;LC_CTYPE=English_United States.1252;LC_MONETARY=English_United States.1252;LC_NUMERIC=C;LC_TIME=English_United States

Re: [R] stepAIC on lm() where response is a matrix..

models and you want a common one, take the models you get, and interpolate manually a list of alternative plausible models between the two best univariate models. Then fit those manually and select the one with the smallest AIC. Hope this helps. Spencer Graves vinod gullu wrote

Re: [R] Gaussian elimination - singular matrix

as Moshe suggested, but I suspect that 'ginv', for example, might return the minimum length solution automatically. Hope this helps. Spencer Graves Moshe Olshansky wrote: All the nontrivial solutions to AX = 0 are the eigenvectors of A corresponding to eigenvalue 0 (try eigen function). The non

Re: [R] nlme correlated random effects

your question, but I hope it helps. Best Wishes, Spencer Graves Daniel O'Shea wrote: I am examining the following nlme model. asymporig-function(x,th1,th2)th1*(1-exp(-exp(th2)*x)) mod1-nlme(fa20~(ah*habdiv+ad*log(d)+ads*ds+ads2*ds2+at*trout)+asymporig(da.p,th1,th2), fixed=ah

Re: [R] PCA for Binary data

functions for 'Multiple Factor Analysis for Mixed [quantitative and qualitative] Data' Hope this helps. Spencer Graves Josh Gilbert wrote: I don't understand, what's wrong with using prcomp in this situation? On Sunday 10 June 2007 12:50 pm, Ranga Chandra Gudivada wrote: Hi

Re: [R] How do I obtain standard error of each estimated coefficients in polr

the code. If this does NOT answer your question, please provide commented, minimal, self-contained, reproducible code, as requested in the posting guide www.R-project.org/posting-guide.html. Hope this helps. Spencer Graves [EMAIL PROTECTED] wrote: Hi, I obtained all

Re: [R] Fwd: Using odesolve to produce non-negative solutions

in line Martin Henry H. Stevens wrote: Hi Jeremy, First, setting hmax to a small number could prevent a large step, if you think that is a problem. Second, however, I don't see how you can get a negative population size when using the log trick. SG: Can lsoda estimate complex or imaginary

Re: [R] Nonlinear Regression

'. That gave me the residuals. # I then tried 'sum(residuals(Pur.wt)^2)', which returned 14.59690. Hope this helps. Spencer Graves p.s. Did this answer your question? Your example did not seem to me to be self contained, which makes it more difficult for me to know if I'm

Re: [R] Fwd: Using odesolve to produce non-negative solutions

in the past. Often, it can even make the log likelihood or sum of squares surface more elliptical, which means that the standard normal approximation for the sampling distribution of parameter estimates will likely be more accurate. Hope this helps. Spencer Graves p.s. Your

Re: [R] factor analysis

. Hope this helps. Spencer Graves p.s. I might have answered this a day or two earlier, but the lack of a simple, self-contained example meant that I would have to work harder to understand your question and craft an answer. bunny , lautloscrew.com wrote: Hi there, i´ve trouble

Re: [R] R-squared in mixed-effects models

RSiteSearch(R^2 in lme) produced 34 hits for me just now, including 5 from an exchange on this list dated 2007.05.14. They should answer your question better than I can. Hope this helps. Spencer Graves Richard Gunton wrote: Hello, I'm fitting general linear models

Re: [R] Can I treat subject as fixed effect in linear model

The short answer is that you could fit that fixed-effect model using 'lm', for example. That would make sense if you wanted to make inference only about that particular group of 20 subjects AND you thought it was inappropriate to consider that their contribution to a model would follow

[R] 'varFunc' class with additive variance? (was: can I get same results using lme and gls?)

(not multiplicative) variance structure. Might there be other tools for modeling an additive variance structure? Alternatively, does it sound sensible to create 'varSumSq' functions similar to 'varComb'? Thanks, Spencer Graves Douglas Bates wrote: On 5/20/07, [EMAIL

Re: [R] [OT] Is data copyrightable?

at least in the US. Spencer Graves Brian's reply seems more consistent with what I've heard than Peter's. The briefest summary I know of copyright law is that expression but not ideas can be copyrighted. Copyright law exists to promote useful arts

Re: [R] Nonlinear constrains with optim

a little outside your constraints, shrink a little the place at which the penalty starts. Hope this helps. Spencer Graves They might be less confused with smaller penalty rates. However if the penalty rate is too small, then you can get a solution that breaks one or more penalties

Re: [R] Optim

for another method. Also consider 'nlminb'. hope this helps. spencer graves Ravi Varadhan wrote: Let us first assume that you have enumerated all the local maxima, which is by no means a trivial thing to assure. How different are the likelihood values? If they are significantly

Re: [R] getting informative error messages

in the 'base' package. I've heard that it has more extensive capabilities, but I've never used it.} I suspect you may already know 'debug', but for those who don't, I think it's worth noting its utility for this kind of thing. Hope this helps. Spencer Graves Tony Plate wrote

Re: [R] getting informative error messages

. 3. Thanks again for your many contributions to the R project and to my education more generally. Best Wishes, Spencer Graves Prof Brian Ripley wrote: It is not clear to me what you want here. Errors are tagged by a 'call', and f(1:3) is the innermost 'call' (special primitives

Re: [R] factanal AIC?

this be something closer to the 0.4755 for fit1$criteria? Thanks, Spencer Graves JENS: For your purposes, I suggest you try to compute (-n/2)*log(det(tcrossprod(fit$loadings)+diag(fit$uniquenesses)). If you do this, you might get more sensible results with your examples

Re: [R] Matlab import

. Spencer Graves Henrik Bengtsson wrote: Hi, as already mentioned, do not save MAT files in ASCII format but save to binary formats, i.e. do *not* use -ascii. Moreover, from ?readMat, you find that: From Matlab v7, _compressed_ MAT version 5 files are used by default [3

Re: [R] Matlab import

Have you considered exporting from Matlab using the '-ascii' option, then reading it with 'read.table', after perhaps checking the first few rows using, e.g., 'readLines', and confirming that all records have the same length using 'count.fields'? Hope this helps. Spencer

Re: [R] R in Industry

whose fee structure is becoming prohibitive. In many situations, just analyzing the legal issues could cost you more than paying someone to modify R code to support your changing needs. Spencer Graves Charilaos Skiadas wrote: On Feb 8, 2007, at 12:48 PM, Ben Fairbank wrote: If my

Re: [R] Rating competitors

Have you considered Bradley-Terry models? RSiteSearch(bradley, functions) just returned 31 hits for me. Hope this helps. Spencer Graves Jeff Newmiller wrote: I am looking for hints on how to estimate ratings for competitors in an ongoing pairwise competition using R

Re: [R] Find S4 Generic?

Dear Prof. Ripley: Thanks very much. That works. I got stuck on the help page for dumpMethods and failed to check See Also. Best Wishes, Spencer Graves Prof Brian Ripley wrote: Do you want E (type 'E') or its methods (getMethods(E) works for me)? On Sun, 17 Dec 2006

[R] Find S4 Generic?

How can I get the R code for E in the distrEx package? The function call 'dumpMethods(E, E.R)' created E.R in the working directory. Unfortunately, it apparently contains 0 bytes. Thanks, Spencer Graves __ R-help

Re: [R] DOE teaching suggestions?

Hi, Erin: Also, have you seen the BHH2 package, companion to Box, Hunter and Hunter (2005) Statistics for Experimenters, 2nd ed. (Wiley)? Hope this helps. Spencer Graves Are you planning to have them design and conduct an actual

Re: [R] DOE teaching suggestions?

Hi, Erin: Also, have you seen the BsMD package (Bayes Screening and Model Discrimination), also discussed in Box Hunger and Hunter (2005). Spencer Graves ## Also, have you seen the BHH2 package, companion to Box, Hunter and Hunter (2005) Statistics

Re: [R] Help with response CIs for lme

it. I know this doesn't solve your problem, but I hope it helps. Spencer Graves Michael Kubovy wrote: Hi, Can someone please offer a procedure for going from CIs produced by intervals.lme() to fixed-effects response CIs. Here's a simple example: library(mlmRev) library

Re: [R] Quadratic Optimization

). SpG Maximize a'x subject to x'Ax=c. Not I believe the usual meaning (nor what Googling 'quadratic constraints' came up with for me). Thanks for the clarification. Spencer Graves __ R-help@stat.math.ethz.ch mailing list https

Re: [R] lme function

in the book [e.g., x^2 must be I(x^2)]. Before I found the script files, I couldn't understand why I got substantially different results from the book when just typing the commands into R. Hope this helps. Spencer Graves Mark Wilson wrote: Hello. As advised by Mick Crawley in his

Re: [R] Quadratic Optimization

the problem. Of course, to do this, you have to make sure your objective function returns valid numbers outside the constrained region. How's this? Spencer Graves amit soni wrote: Hi, I need to solve an optimization problem in R having linear objective function and quadratic

Re: [R] help

as commercial software. Hope this helps. Spencer Graves Aimin Yan wrote: consider p as random effect with 5 levels, what is difference between these two models? p5.random.p - lmer(Y ~p+(1|p),data=p5,family=binomial,control=list(usePQL=FALSE,msV=1)) p5.random.p1 - lmer(Y ~1+(1|p

Re: [R] problem with the plm package (2)

until I found this. Hope this helps. Spencer Graves Giangiacomo Bravo wrote: Thanks a lot for the quick answerd I received and that helped me to solve my first (simple) problem with the plm package. Unfortunatley, here's another one: I do not know why, but I'm unable to estimate

Re: [R] My own correlation structure with nlme

in your installation. If I got the same error message, there is a good chance that I figure out how to get around that and provide more helpful suggestions in less time than I've been thinking about this. Hope this helps. Spencer Graves Mohand wrote: Dear all, I am trying

Re: [R] Nonlinear statistical modeling -- a comparison of R and AD Model Builder

. The CRAN summary of that package led me to an Adobe Acrobat version of a PPT slide presentation that seemed to consider just this type of problem (e.g., http://mcmcpack.wustl.edu/files/MartinQuinnMCMCpackslides.pdf). Have you considered that? Hope this helps. Spencer Graves

Re: [R] how to forecast the GARCH volatility?

are probably described many places, but one I just checked was sec. 3.1 in Venables and Ripley (2000) S Programming (Springer). Hope this helps. Spencer Graves __ R-help@stat.math.ethz.ch mailing list https://stat.ethz.ch/mailman

Re: [R] SNA packages/network package

) with the 'sna' package in the second position on the search path, I got 201 different objects. I didn't see an obvious entry point for the package. Hope this helps. Spencer Graves Bagatti Davide wrote: Hello everyone, I am an italian student who is working with packages SNA (social

Re: [R] stl and the relative scale data-seasonal

. Hope this helps. Spencer Graves p.s. By 'debug' here, I mean the 'debug' in the 'base' package, not the 'debug' package, which I've never used but is doubtless more powerful (though may also require more effort and thought from the user). For more information on debug{base}, I

Re: [R] nlme

and Its Applications (Wiley). Spencer Graves Douglas Bates wrote: On 11/19/06, Fluss [EMAIL PROTECTED] wrote: Hello! I am trying to fit a mixed non-linear model using nlme. How can I constrain the fixed parameter space (add bounds) as in nls. By rewriting the nlme function

Re: [R] Need help in nlmeODE

your needs (as suggested in the posting guide www.R-project.org/posting-guide.html). Hope this helps. Spencer Graves Ron Fisher wrote: I am trying to use nlmeODE. But I don't know what is ObsEq. How you can get it from your ODEs system? Could someone help me out? Best, Ron

Re: [R] deriv when one term is indexed

case. Sorry I couldn't be more encouraging. Spencer Graves (*) Reference: library(fortunes) Ftns - read.fortunes() sel - regexpr(This is R, Ftns$quote) which(sel0) [1] 109 fortunes(109) Ken Knoblauch wrote: I don't know why I get a different result than you do, though I'll play

Re: [R] Repeated measures by lme and aov give different results

. Spencer Graves Vicki Allison wrote: I am analyzing data from an experiment with two factors: Carbon (+/-) and O3 (+/-), with 4 replicates of each treatment, and 4 harvests over a year. The treatments are assigned in a block design to individual Rings. I have approaches

[R] nmle for time of structural change?

step functions or linear or quadratics in time with any number of breakpoints. I'm currently estimating the required parameters using the 'strucchange' package. That work, but I'm having trouble with this enhancement. Thanks, Spencer Graves

Re: [R] Confidence interval for relative risk

*log(likelihood ratio) being approximately chi-square? Just my 2e-9 Euros (or 2e-7 Yen or Yuan). Spencer Graves Michael Dewey wrote: At 12:35 12/11/2006, Peter Dalgaard wrote: Michael Dewey [EMAIL PROTECTED] writes: At 15:54 10/11/2006, Viechtbauer Wolfgang (STAT

Re: [R] POSIXlt converted to POSIXct in as.data.frame()

)? Hope this helps. Spencer Graves Roger Bivand wrote: In trying to use as.Date(), I've come across the conversion of POSIXlt to POSIXct when a POSIXlt variable is included in a data frame: my_POSIX - strptime(c(11-09-2006, 11-10-2006, 11-11-2006, 11-12-2006, 11-13-2006), %m

Re: [R] Problems with metaMDS from vegan

Have you sent this question to the vegan maintainer, identified, e.g., with help(package=vegan)? I've added that name as a cc to this email. Hope this helps. Spencer Graves Peter Roosen wrote: Hello all, I recently used the Vegan library quite extensively

Re: [R] Value at Risk historical simulation

this helps. Spencer Graves Benjamin Dickgiesser wrote: Hi Has someone got a package/script at hand to do a historical simulation to calculate the Value at Risk? If your not sure what Historical Simulation is: In simple terms, Historical Simulation (HS) is just taking sample

Re: [R] help with nlme function

will likely find the problem. Hope this helps. Spencer Graves Bill Shipley wrote: Hello. I am trying to fit a nonlinear mixed model involving 3 parameters. I have successfully made a self-starting function. getInitial() correctly outputs the initial estimates. I can also use the nlsList

Re: [R] Problems with metaMDS from vegan

. Spencer Graves Peter Roosen wrote: Hello Spencer, Have you sent this question to the vegan maintainer, identified, e.g., with help(package=vegan)? I wrote a mail directly to Jari prior to posting in this group. My posting here was no sooner than waiting (in vain

Re: [R] Variable limit in nlm?

much simpler (and self contained)? If that does not lead you to an answer to your question, it should at least provide something that could help someone else help you. I know my comments don't solve your problem, but hopefully they will get you closer. Spencer Graves Robert

Re: [R] Comparison between GARCH and ARMA

, as suggested in the posting guide www.R-project.org/posting-guide.html. Hope this helps. Spencer Graves [EMAIL PROTECTED] wrote: Hi i was just going by this thread, i thought of igniting my mind and got something wierd so i thought of making it wired. i think whether you take

Re: [R] Making a case for using R in Academia

of this latter point, I now look for an R companion for technical books and articles I read: I will learn more, faster if R code exists than if it doesn't. Bottom line: You need to prepare your students for the future, not the past. Hope this helps. Spencer Graves Doran

Re: [R] Comparing models in multiple regression and hierarchical linear regression

could consult or even possibly collaborate on this? Hope this helps. Spencer Graves Jenifer Larson-Hall wrote: I don’t know if this question properly belongs on this list, but I’ll ask it here because I’ve been using R to run linear regression models, and it is only in using R

Re: [R] multivariate splines

): It produced 45 hits, some of which might interest you. Hope this helps. Spencer Graves p.s. I received permission from Paul Dierckx to build an R package using the the software associated with his book, Curve and Surface Fitting with Splines (Oxford University Press, 1993

Re: [R] Which genetic optimization package allows customized crossover and mutation operation

Are you familiar with www.bioconductor.org? They have a listserve with people who may know more about the question you asked. Spencer Graves sun wrote: Hi, I am looking for genetic optimization package that allow to define my own chromosome solution and crossover/mutation

Re: [R] Measuring the effects of history on event probabilities

observations on 2-4 individuals. Then tell us what you'd like to do in terms of this toy example, possibly including things you've tried and why they didn't seem to give you what you wanted. Hope this helps. Spencer Graves Jon Minton wrote: This is probably very simple but my brain

Re: [R] help went away from top emacs window :-(

over the world. Best Wishes, Spencer Graves Richard M. Heiberger wrote: I agree with both Duncan and Martin. With Duncan, that this is not a Windows FAQ. With Martin, that it is an FAQ. I recommend that it go on the R FAQ page, prior to the current 7.24. Here is a draft

Re: [R] help went away from top emacs window :-(

I concur: Thanks very much for looking into this. Spencer Graves Duncan Murdoch wrote: Thanks for looking into this. On 11/3/2006 9:53 AM, Richard M. Heiberger wrote: I located what might be the problem. There is a different misbehavior with Gnu emacs. I don't have Xemacs. I have

Re: [R] nlme Error: Subscript out of bounds

a chain of earlier recommendations on debug, I tried 'RSiteSearch(graves debug), then 'RSiteSearch(graves debug lme)'. The fifth hit from the latter is http://finzi.psych.upenn.edu/R/Rhelp02a/archive/77298.html;. Hope this helps. Spencer Graves Lisa Avery wrote: Hello, I am new to non

Re: [R] Variance Component/ICC Confidence Intervals via Bootstrap or Jackknife

probably try the second, though I might try both and compare. If you try them both, I'd be interested in the comparison. If the answers were substantially different, I'd worry. Hope this helps. Spencer Graves Rick Bilonick wrote: I'm using the lme function in nmle to estimate

Re: [R] Mixed conditional logit model

.) Hope this help. Spencer Graves Ryuichi Tamura wrote: Hi, Sven Müller sven.mueller at tu-dresden.de writes: i wonder whether it is possible to estimate a mixed (random parameters) logit model in R. I wrote a collection of R functions for estimating discrete choice

Re: [R] Help with random effects and smoothing splines in GAMM

in the posting guide www.R-project.org/posting-guide.html. Your example was not self contained. Hope this helps. Spencer Graves Liang, Hua wrote: Try to fit a longitudinal dataset using generalized mixed effects models via the R function gamm() as follows: library(mgcv) gamm0.fit

Re: [R] Get the names of the columns in a tserie

it will likely take me a lot longer to even understand what you are asking, and I will likely not be very confident that I even understand your question even if I attempt an answer. Hope this helps. Spencer Graves lvdtime wrote: Please, I need this information, it's important

Re: [R] Variance Component/ICC Confidence Intervals via Bootstrap or Jackknife

see in line Rick Bilonick wrote: On Sun, 2006-10-29 at 11:06 -0800, Spencer Graves wrote: I can think of two ways to get confidence intervals on intraclass correlations (ICCs) and more accurate intervals for variance snip The way I've done the bootstrapping

Re: [R] plot.POSIXct plot.POSIXlt

;. Hope this helps. Spencer Graves Patrick Giraudoux wrote: Hi, Just to signal that when I want to plot POSIXct variable on x using format within plot(), I get what I want on the plot but with a number of warnings: plot(y~x,format=%y-%m) Warning messages: 1: format

Re: [R] Tukey-Kramer test

-project.org/posting-guide.html. Hope this helps. Spencer Graves A.R. Criswell wrote: Hello All, I found the TukeyHSD() function. Is there a Tukey-Kramer test for unbalanced data? Andrew __ R-help@stat.math.ethz.ch mailing list https

Re: [R] linear mixed effects models with breakpoints

that to the 'lmeChgPt' function. If you'd like further help from this listserve, please provide commented, minimal, self-contained, reproducible code, as suggested in the posting guide www.R-project.org/posting-guide.html. Hope this helps. Spencer Graves Ewan Wakefield wrote: Hi

Re: [R] MARS help?

. Spencer Graves Liaw, Andy wrote: Spencer, MARS fits splines, not disconnected lines. Perhaps the strucchange package has facility to fit your data better. Cheers, Andy *From:* [EMAIL PROTECTED] on behalf

[R] MARS help?

traces in terms of a fixed number of knots and associated slopes for each linear spline segment max(0, t1-t.knot). Thanks, Spencer Graves __ R-help@stat.math.ethz.ch mailing list https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read

Re: [R] Interval pls

. Spencer Graves Dirk De Becker wrote: Hey R-helpers, Does anybody know of an implementation of interval PLS in R? Thx, Dirk __ R-help@stat.math.ethz.ch mailing list https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read the posting guide http

Re: [R] extracting nested variances from lme4 model

and ch05.R in ~library\nlme\scripts, it might take me a while to figure out how to do it (and even answer the question of whether it can). Thanks again. Spencer Graves Douglas Bates wrote: On 10/15/06, Douglas Bates [EMAIL PROTECTED] wrote: On 10/14/06, Spencer Graves [EMAIL

Re: [R] split-plot analysis with lme()

a method for estimating an interaction term with 5 degrees of freedom, I hope I've at least answered your basic question. Best Wishes, Spencer Graves i.m.s.white wrote: Dear R-help, Why can't lme cope with an incomplete whole plot when analysing a split-plot experiment

Re: [R] Gradient problem in nlm

the quantity). Hope this helps. Spencer Graves singyee ling wrote: Hello everyone! I am having some trouble supplying the gradient function to nlm in R for windows version 2.2.1. What follows are the R-code I use: fredcs39-function(a1,b1,b2,x){return(a1+exp(b1+b2*x))} loglikcs39

Re: [R] nlme_varcov matrix

. If you'd like further help from this listserve, please provide commented, minimal, self-contained, reproducible code as suggested in the posting guide www.R-project.org/posting-guide.html. Para servirle. Spencer Graves gabriela escati peñaloza wrote: Hi all

Re: [R] Optim: Function definition

Have you tried making contour plots of your objective function, e.g., using expand.grid and 'contour' or 'contourplot', as described in Venables and Ripley (2002) Modern Applied Statistics with S (Springer)? Hope this helps. Spencer Graves Serguei Kaniovski wrote: Hi all

Re: [R] extracting nested variances from lme4 model

like 'trth' in your example. If we are lucky, someone else might educate us both. I'm sorry I couldn't answer your question, but I hope these comments might help in some way. Spencer Graves Frank Samuelson wrote: I have a model: mod1-lmer( x ~ (1|rtr)+ trth/(1|cs) , data

Re: [R] Nonlinear mixed effects model

ways to do things like this, and this book is probably the best-known product of that work. Hope this helps. Spencer Graves Cam wrote: In nonlinear mixed effects models, SAS doesn't allow for free manipulation of the covariance matrix (you can only specify a type, and our type

Re: [R] mixed models: correlation between fixed and random effects??

I suspect you've observed an outcome that occurs once out of every 2^4 = 16 times. Have you tried Monte Carlo, e.g., using 'simulate.lme' in library(nlme) or 'mcmcsamp' in library(lme4)? Hope this helps. Spencer Graves Bruno L. Giordano wrote: Hello, I built 4 mixed

Re: [R] Once again: aov vs. lme/lmer

' in the R Wiki. Hope this helps. Spencer Graves Rafael Laboissiere wrote: First of all, I apologize for asking a question that has appeared recurrently in this mailing list. However, I have googled for it, have looked at the mailing list archives, and also looked at Pinheiro Bates book

Re: [R] Helmert contrasts for repeated measures and split-plot expts

both on how sensible this might be. Hope this helps. Spencer Graves Many thanks in advance Roy --- Roy Sanderson Institute for Research on Environment and Sustainability Devonshire Building

Re: [R] Plackett-Dale Model in R

. Spencer Graves Pryseley Assam wrote: Dear R users, Can someone inform me about a library/function in R that fits a Plackett-Dale model ? Thanks in advance Pryseley - [[alternative HTML version deleted

Re: [R] warning message in nlm

Statistics with S (Springer), discussing 'expand.grid' and 'contour'. Hope this helps. Spencer Graves singyee ling wrote: Dear R-users, I am trying to find the MLEs for a loglikelihood function (loglikcs39) and tried using both optim and nlm. fredcs39-function(b1,b2,x){return(exp(b1

Re: [R] treatment effect at specific time point within mixed effects model

1.4340322 10 -0.4474502 0.6641 factor(Time)6:drugP -2.1396105 1.4340322 10 -1.4920240 0.1666 Does this answer your question? Hope this helps. Spencer Graves Afshartous, David wrote: All, The code below is for a pseudo dataset of repeated measures on patients where

Re: [R] variance-covariance structure of random effects in lme

. Hope this helps. Spencer Graves Mi, Deming wrote: Dear R users, I have a question about the patterned variance-covariance structure for the random effects in linear mixed effect model. I am reading section 4.2.2 of Mixed-Effects Models in S and S-Plus by Jose Pinheiro and Douglas Bates

Re: [R] variance functions in glmmPQL or glm?

-guide.html. Hope this helps, even if it doesn't answer your question. Spencer Graves Gretchen wrote: Hello R users- I am new to R, and tried searching the archives and literature for an answer to this - please be patient if I missed something obvious. I am fitting a logistic

Re: [R] Need help to estimate the Coef matrices in mAr

data set, e.g., 10 simulated observations on 2 variables, with comments describing what you tried and what you don't understand about it. Hope this helps. Spencer Graves Arun Kumar Saha wrote: Dear Spencer, Thank you very much for your attention on my problem. According

Re: [R] plotting all subgroups with augPred

need, you can plot them any way you want. Hope this helps. Spencer Graves Afshartous, David wrote: All, I have a question RE plotting the prediction lines of a random effects model via augPred. I'll illustrate via the Pixel dataset within the nlme package: library(nlme

Re: [R] EBAM ERROR

suggest you send your email From something more revealing of who you are. Some of the most knowledgeable and frequent contributors to this listserve rarely if ever answer questions from anonymous email addresses like [EMAIL PROTECTED]. Hope this helps. Spencer Graves Learning R

Re: [R] Formula aruguments with NLS and model.frame()

this, please submit another post including commented, minimal, self-contained, reproducible code, as suggested in the posting guide 'www.R-project.org/posting-guide.html'. Hope this helps. Spencer Graves and provide commented, minimal, self-contained, reproducible code. Joe W. Byers

Re: [R] time series simulation

this listserve, please provide commented, minimal, self-contained, reproducible code, as suggested in the posting guide www.R-project.org/posting-guide.html. Hope this helps. Spencer Graves march wrote: Hi everybody I'm trying to simulate a stochastic process in R. I would like consider n

Re: [R] Simulation of a Dirichlet Process.

MCMC)', try something. If it is not what you want, please convert it to a commented, minimal, self-contained, reproducible code and send that to this listserve, as suggested in the posting guide www.R-project.org/posting-guide.html. Hope this helps. Spencer Graves [EMAIL PROTECTED

Re: [R] ISO8601 week-of-year to date

describing something you've tried and did not find adequate (as suggested in the posting guide 'www.R-project.org/posting-guide.html'). Hope this helps. Spencer Graves Ott-Siim Toomet wrote: Hi, are there any way to convert ISO8601 weeks to gregorian dates? Something like coverttodate

Re: [R] Power analysis for repeated measures ANCOVA

I haven't seen a reply to this post, and you may have moved beyond this question by now, but I will nevertheless offer one brief question: Have you looked at 'simulate.lme'? This might simplify coding your power analysis. Hope this helps. Spencer Graves Steve Su wrote

Re: [R] nls

and relatively simple example. Without it, I could have done little more that suggest 'traceback' (which provided zero information when I tried it) and 'debug'. Hope this helps. Spencer Graves fit.gaus-function(coordinates,values,guess.c0,guess.c1,guess.a) { long-rep(coordinates[,1

  1   2   3   4   5   6   7   8   9   10   >