Re: [R] how do I make a correlation matrix positive definite?

not reasonable to assume multivariate normal, then I don't know if it's reasonable to estimate a correlation matrix.) I have not tried it, but I believe that missMDA should solve this problem assuming all the uniquenesses are equal. Hope this helps, Spencer Graves On 10/21/2010 5:03 PM

Re: [R] visualize TukeyHSD results

, Torsten Hothorn and Peter Westfall (2010), Multiple Comparisons Using R, CRC Press, Boca Raton. Hope this helps. Spencer Graves On 10/21/2010 7:50 PM, Timothy Spier wrote: I am a new R user but a long time SAS user. I searched for a response to this question but no luck, so forgive me

Re: [R] cube root of a negative number

install.packages('sos')# if you don't have it already library(sos) rs - ???roots # 216 matches summary(rs) # in 106 packages rs # opens a web browser with all 216 matched in a table # listing the package with the most matches first. # This included roots{signal}, which referenced polyroot{base},

Re: [R] Constrained Regression

__ R-help@r-project.org mailing list https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read the posting guide http://www.R-project.org/posting-guide.html and provide commented, minimal, self-contained, reproducible code. -- Spencer Graves, PE, PhD President

Re: [R] Constrained Regression

... On Sun, Oct 31, 2010 at 1:44 PM, David Winsemiusdwinsem...@comcast.netwrote: On Oct 31, 2010, at 12:54 PM, Spencer Graves wrote: Have you tried the 'sos' package? I have, and I am taking this opportunity to load it with my .Rprofile to make it more accessible. It works very well. Very clean display

Re: [R] RMySQL on Windows 2008 64 Bit -Help!

://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read the posting guide http://www.R-project.org/posting-guide.html and provide commented, minimal, self-contained, reproducible code. -- Spencer Graves, PE, PhD President and Chief Operating Officer Structure Inspection and Monitoring, Inc. 751

Re: [R] where are my pspline knots?

Hello: I hope that someone more knowledgeable will confirm or correct what I'm about to say. I don't have time now to check this by studying the pspline code. From reading the pspline{survival} help page, I believe it uses standard B-splines, penalising the integrated second

Re: [R] kalman filter in sspir

To learn why sspir does not have a filter function, you need to ask the package maintainer, Claus Dethlefsen c...@rn.dk. My belief is that he, Soren Lundbye-Christensen and Anette Luther Christensen found other outlets for their time since they completed the package and the companion

Re: [R] a philosophy R question

...@gmail.com __ R-help@r-project.org mailing list https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read the posting guide http://www.R-project.org/posting-guide.html and provide commented, minimal, self-contained, reproducible code. -- Spencer

Re: [R] Can't invert matrix

mailing list https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read the posting guide http://www.R-project.org/posting-guide.html and provide commented, minimal, self-contained, reproducible code. -- Spencer Graves, PE, PhD President and Chief Operating Officer Structure Inspection

Re: [R] Find in R and R books

Other people like R Site Search (http://search.r-project.org/nmz.html), which is available via the standard R function RSiteSearch. For me, the fastest literature search on virtually anything statistical is the findFn function in the sos package. (Disclaimer: I'm the lead

Re: [R] Find in R and R books

Hi, Mike, et al.: in line On 11/22/2010 5:43 PM, Mike Marchywka wrote: Date: Mon, 22 Nov 2010 17:28:57 -0800 From: spencer.gra...@structuremonitoring.com To: g...@well.ox.ac.uk CC: r-h...@stat.math.ethz.ch Subject: Re: [R] Find in R and R books

Re: [R] how to calculate derivative

. __ R-help@r-project.org mailing list https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read the posting guide http://www.R-project.org/posting-guide.html and provide commented, minimal, self-contained, reproducible code. -- Spencer Graves, PE, PhD President

Re: [R] help

://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read the posting guide http://www.R-project.org/posting-guide.html and provide commented, minimal, self-contained, reproducible code. -- Spencer Graves, PE, PhD President and Chief Operating Officer Structure Inspection and Monitoring, Inc. 751

Re: [R] looking for the RMySQL package for R 2.12.0 under XP

Hello: That package needs to be compiled from source on the computer you are using. Instructions for doing so can be found at http://biostat.mc.vanderbilt.edu/wiki/Main/RMySQL;. I'm not sure of all the issues, but I believe it needs to be configured to the specific version of

Re: [R] Apropos the day...

ditto. Spencer Graves On 11/25/2010 1:20 PM, Dennis Murphy wrote: Since today is American Thanksgiving, I want to thank: (a) R-core for all of their efforts to produce what is, IMHO, the best statistical software around, not simply for the convenience of doing more, better

Re: [R] Where is gdata?

/posting-guide.html and provide commented, minimal, self-contained, reproducible code. -- Spencer Graves, PE, PhD President and Chief Operating Officer Structure Inspection and Monitoring, Inc. 751 Emerson Ct. San José, CA 95126 ph: 408-655-4567

Re: [R] Where is gdata?

and install the perl gdata: libaries needed to support Excel XLS and XLSX formats. Read ?installXLSXsupport -- Spencer Graves, PE, PhD President and Chief Operating Officer Structure Inspection and Monitoring, Inc. 751 Emerson Ct. San José, CA 95126 ph: 408-655-4567

Re: [R] how to divide each column in a matrix by its colSums?

prop.table divides every element by the matrix total, not its colSums: m - matrix(1:4, 2) m [,1] [,2] [1,]13 [2,]24 prop.table(m) [,1] [,2] [1,] 0.1 0.3 [2,] 0.2 0.4 m/rowSums(m) divides every row by its rowSum: m/rowSums(m) [,1] [,2]

Re: [R] how to divide each column in a matrix by its colSums?

Thanks: I didn't read the prop.table help page with sufficient care. On 11/28/2010 1:48 PM, casperyc wrote: I am using prop.table(m,1) and prop.table(m,2) my aim. which I think is the most 'easy' way. Thanks. Casper -- Spencer Graves, PE, PhD President and Chief Operating Officer

Re: [R] Where is gdata?

. Then everything seemed to work appropriately. Thanks to all. Spencer On 11/28/2010 7:49 AM, Gabor Grothendieck wrote: On Sun, Nov 28, 2010 at 10:40 AM, Spencer Graves spencer.gra...@structuremonitoring.com wrote: I tried install.packages('gregmisc'), which Stephen said he had

Re: [R] Where is gdata?

to all. Spencer On 11/28/2010 7:49 AM, Gabor Grothendieck wrote: On Sun, Nov 28, 2010 at 10:40 AM, Spencer Graves spencer.gra...@structuremonitoring.comwrote: I tried install.packages('gregmisc'), which Stephen said he had. Then help(pac=gdata) seemed to work normally

Re: [R] Where is gdata?

commented, minimal, self-contained, reproducible code. -- Spencer Graves, PE, PhD President and Chief Operating Officer Structure Inspection and Monitoring, Inc. 751 Emerson Ct. San José, CA 95126 ph: 408-655-4567 __ R-help@r-project.org mailing list

Re: [R] Where is gdata?

://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read the posting guide http://www.R-project.org/posting-guide.html and provide commented, minimal, self-contained, reproducible code. -- Spencer Graves, PE, PhD President and Chief Operating Officer Structure Inspection and Monitoring, Inc. 751

Re: [R] [R-pkgs] New package list for analyzing list survey experiments

identify help pages whose name contains a particular term. The R Journal from last December contains an article describing this: http://journal.r-project.org/archive/2009-2/RJournal_2009-2_Graves~et~al.pdf;. Hope this helps. Spencer Graves On 7/12/2010 7:08 AM

Re: [R] [R-pkgs] New package list for analyzing list surveyexperiments

://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read the posting guide http://www.R-project.org/posting-guide.html and provide commented, minimal, self-contained, reproducible code. -- Spencer Graves, PE, PhD President and Chief Operating Officer Structure Inspection and Monitoring, Inc. 751

Re: [R] Odp: R forum

faster than posting a question to r-help. Hope this helps. Spencer Graves On 7/22/2010 12:39 AM, Petr PIKAL wrote: Hi r-help-boun...@r-project.org napsal dne 22.07.2010 05:58:54: Hi, Are there any other forums on R where I can get quick responses to my questions? I've lot

Re: [R] Hurst Exponent Estimation

functions are designed to produce an Excel file with a summary page that can help you identify the package that seems to be most actively maintained among the relevant packages, as explained in a vignette. Hope this helps. Spencer Graves On 7/22/2010 12:25 AM, Hans W Borchers wrote

Re: [R] tsdiag

you ask. This may not help, but it might. Best Wishes, Spencer On 7/22/2010 5:00 AM, nuncio m wrote: HI list, I want to know whether tsdiag uses k-(p+q) as the lag in ljung box test. How is it possible to save those values nuncio -- Spencer Graves, PE, PhD President

Re: [R] optimization subject to constraints

in the sos vignette. A shameless plug from the lead author of the sos package. Spencer Graves On 8/9/2010 10:01 AM, Ravi Varadhan wrote: constrOptim can only handle linear inequality constraints. It cannot handle equality (linear or nonlinear) as well as nonlinear inequality constraints

Re: [R] Kalman filter

other than the a linear combination of components of the state vector plus normal noise, then I suggest you consider the sspir package, which is number 6 on this list in terms of the number of help pages matching the search terms. Hope this helps. Spencer Graves On 8/13/2010 2:06

Re: [R] R reports

in a variety of ways for whatever you want. It has union capabilities to support combining the results of different searches. Spencer Graves, lead author of sos On 8/19/2010 6:49 AM, Frank Harrell wrote: What do low level proc print and proc report have on Sweave or http

Re: [R] Aspect Ratio

The documentation is not clear. It would help if it had an example like the following: plot(1:2, 1:2/10) plot(1:2, 1:2/10, asp=1) Does looking at these two plots answer the question? Spencer Graves On 8/19/2010 2:36 PM, David Winsemius wrote: On Aug 19, 2010, at 5

Re: [R] Aspect Ratio

of different examples until sufficient enlightenment is achieved. Best Wishes, Spencer Graves On 8/19/2010 3:38 PM, r.ookie wrote: Thanks for your example as well. Ted's example was exactly what I needed. On Aug 19, 2010, at 3:18 PM, Gavin Simpson wrote: On Thu, 2010-08-19 at 14:28 -0700

Re: [R] First and second derivatives of raw observations

). The fda package includes scripts to work all but one of the 76 examples in the book. Hope this helps. Spencer Graves On 8/30/2010 4:13 PM, Duncan Murdoch wrote: On 30/08/2010 6:40 PM, mtor...@math.carleton.ca wrote: Dear all, I was asked to send the following question: We have

Re: [R] First and second derivatives of raw observations

package includes scripts to work all but one of the 76 examples in the book. Hope this helps. Spencer Graves On 8/30/2010 4:13 PM, Duncan Murdoch wrote: On 30/08/2010 6:40 PM, mtor...@math.carleton.ca wrote: Dear all, I was asked to send the following question: We have some (raw

Re: [R] dirichlet models

to that produced in the browser. Hope this helps. Spencer Graves (lead author of sos) On 9/5/2010 8:07 AM, David Winsemius wrote: On Sep 5, 2010, at 10:53 AM, Donald Braman wrote: Does anyone know of a package (or workaround) for fitting a dirichlet distribution by maximum likelihood

Re: [R] queue implementation?

an extended summary by package in one sheet with the search results in another, documented in a vignette. Hope this helps. Spencer Graves lead author of sos On 9/7/2010 12:56 AM, Barry Rowlingson wrote: On Tue, Sep 7, 2010 at 8:11 AM, raje...@cse.iitm.ac.in raje...@cse.iitm.ac.in

[R] Wiki entries on package development process and software repository

comments, suggested improvements, etc., via email. Thanks, Spencer Graves __ R-help@r-project.org mailing list https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read the posting guide http://www.R-project.org/posting-guide.html and provide

Re: [R] Dampening the spline interpolation for contours

levels. The output of predict can then be fed to contour to produce what you want, I think. I'm sorry I don't have time now to fill in the details. Hope this helps. Spencer Graves p.s. Disclaimer: I'm the lead author of the sos package and third author on the Functional Data

Re: [R] The future of R - Ross Ihaka stirs discussions around the web

commented, minimal, self-contained, reproducible code. -- Spencer Graves, PE, PhD President and Chief Operating Officer Structure Inspection and Monitoring, Inc. 751 Emerson Ct. San José, CA 95126 ph: 408-655-4567 [[alternative HTML version deleted

Re: [R] Prime Factorization

? Yes. At least two. The obvious search strategy with your favrite search tool should work well. -- Spencer Graves, PE, PhD President and Chief Operating Officer Structure Inspection and Monitoring, Inc. 751 Emerson Ct. San José, CA 95126 ph: 408-655-4567

Re: [R] How many R packages are not free?

]] __ R-help@r-project.org mailing list https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read the posting guide http://www.R-project.org/posting-guide.html and provide commented, minimal, self-contained, reproducible code. -- Spencer Graves, PE, PhD President and Chief Operating Officer

Re: [R] How many R packages are not free?

Thanks. Spencer On 10/2/2010 10:51 AM, Dirk Eddelbuettel wrote: On Sat, Oct 02, 2010 at 10:38:22AM -0700, Spencer Graves wrote: Is there anything on CRAN that is NOT free? Lots. Also notice that 'free' means more than redistribution. I assumed that CRAN had a policy

Re: [R] Debye function

* page 12 of reference manual of gsl. http://cran.r-project.org/web/packages/gsl/gsl.pdf -- Spencer Graves, PE, PhD President and Chief Operating Officer Structure Inspection and Monitoring, Inc. 751 Emerson Ct. San José, CA 95126 ph: 408-655-4567 __ R

Re: [R] Tukey HSD Test als Post Hoc Test nach einem GLM inkl. Anova

to receive the news when it was removed from CRAN, because the package had an obsolete email address for me. The package will be back on CRAN by the end of this month. In the meantime, I trust you can install.packages from the attached zip file. Hope this helps. Spencer Graves

Re: [R] Package for converting R datasets into SQL Server (create table and insert statements)?

-project.org/posting-guide.html and provide commented, minimal, self-contained, reproducible code. -- Spencer Graves, PE, PhD President and Chief Operating Officer Structure Inspection and Monitoring, Inc. 751 Emerson Ct. San José, CA 95126 ph: 408-655-4567

Re: [R] GPS data!

-help@r-project.org mailing list https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read the posting guide http://www.R-project.org/posting-guide.html and provide commented, minimal, self-contained, reproducible code. -- Spencer Graves, PE, PhD President and Chief Operating Officer Structure

Re: [R] Rmysql - dbWritetable

, b20101010, x2010, overwrite = T) If this does not work, please note that there is a mailing list specifically for questions regarding database issues with R (R-SIG-DB: https://stat.ethz.ch/mailman/listinfo/r-sig-db;). Hope this helps, Spencer Graves On 10/12/2010 5:57 PM, stephen

[R] Interrupt R?

? Beyond this, what would you suggest to update a real-time report when new data arrives in a certain directory? A generalization of the above works, but I'd like something more graceful. Thanks, Spencer Graves sessionInfo() R version 2.11.1 (2010-05-31) i386-pc-mingw32 locale

Re: [R] Interrupt R?

more skill from the programmer or more user training than using escape under Rgui or ctrl-g/c under Emacs. Best Wishes, Spencer On 7/11/2010 12:02 PM, Duncan Murdoch wrote: On 11/07/2010 2:29 PM, Spencer Graves wrote: How can one interrupt the following gracefully: while

Re: [R] ARCH LM test for univariant time series

the function -- or use a function in another package. Tsay (2005) Analysis of Financial Time Series, 2nd ed. (Wiley) includes an example on p. 103 that could be used for a reference. Hope this helps. Spencer Graves tom soyer wrote: Hi, Does anyone know if R has a Lagrange multiplier

Re: [R] ARCH LM test for univariant time series

= denominator degrees of freedom; http://en.wikipedia.org/wiki/F-distribution], while the expected value for a chi-square is the number of degrees of freedom Unfortunately, I don't feel I can afford the time to dig into this further right now. Thanks for your help. Spencer Graves

Re: [R] ARCH LM test for univariant time series

with this. Spencer Graves tom soyer wrote: Spencer, Sorry, I forgot that the default lag in arch is 16. Here is the fix. Can you try it again and see if it gives the correct (or at least similar compared to a true LM test) result? archTest=function(x, lags=12){ #x is a vector require(vars) s=embed(x

Re: [R] ARCH LM test for univariant time series

inflations, Econometrica 50: 987-1007, cited by Tsay (p. 101). However, that's not on my critical path for today. SG On 2/2/08, Spencer Graves [EMAIL PROTECTED] wrote: Dear Tom: Your revised function eliminates the discrepancy in the degrees of freedom but is still very different from

Re: [R] ARCH LM test for univariant time series

work... hmm... On 2/2/08, Spencer Graves [EMAIL PROTECTED] wrote: Dear Tom: Your revised function eliminates the discrepancy in the degrees of freedom but is still very different from the numbers reports on Tsay, p. 102: archTest(log(1+as.numeric

Re: [R] how to connect S-Plus to Matlab?

1. Have you tried a search at http://www.insightful.com/;? 2. Have you asked Insightful tech support and / or [EMAIL PROTECTED]? 3. Do you have S-PLUS 8? Insightful claims that with S-PLUS 8, open source analytics developed in R ... can be imported

Re: [R] [OT] good reference for mixed models and EM algorithm

with Pinheiro and Bates. Spencer Graves Erin Hodgess wrote: Dear R People: Sorry for the off-topic. Could someone recommend a good reference for using the EM algorithm on mixed models, please? I've been looking and there are so many of them. Perhaps someone here can narrow things

Re: [R] Using R in a university course: dealing with proposal comments

can learn much quicker how to use, modify and extend the statistical tools discussed. Spencer Graves [EMAIL PROTECTED] wrote: Comment 1 raises a real issue. R is just a tool. Too often people do confuse the tool with the real skill that the people who use it should have

Re: [R] Is there a simple way to use one-way ANOVA comparing the means of groups?

see also the 'multcomp' and 'multcompView' packages. Henrique Dallazuanna wrote: See ?TukeyHSD On 13/02/2008, Kes Knave [EMAIL PROTECTED] wrote: Dear all, I have finally managed to get a Analysis of Variance Table: Response: LogHand Df Sum Sq MeanSq F

[R] [R-pkgs] FinTS_0.2-7

of many ACF plots in Tsay (2005) Analysis of Financial Time Series, 2nd ed. (Wiley) that follow this secondary standard. Spencer Graves ___ R-packages mailing list [EMAIL PROTECTED] https://stat.ethz.ch/mailman/listinfo/r-packages

Re: [R] History of R

Hi, Kathy, John, et al.: Has there been an answer to the question of why R has been much more successful than Octave? In this regard, can anyone provide a price comparison of student versions for S-Plus and Matlab during R's gestation period, 10-15 years ago? I had the

Re: [R] Kalman Filter

Have you looked at the 'dlm' package? It has a vignette to help you learn to use it. Also, I've heard that a book about that package is scheduled to appear in the next few months. I have looked at the Kalman functions in the 'stats' package but have not found documentation that

Re: [R] About R-Project

I'm not familiar with the Schuster system, but RSiteSearch(Schuster) produced 34 hits and RSiteSearch(river, fun) produced 133 for me just now. If this does not lead you to an answer to your question, PLEASE do read the posting guide http://www.R-project.org/posting-guide.html and

Re: [R] differentiating a numeric vector

. It may include more than you want to consider, but for me, too much is usually better than nothing. hope this helps. Spencer Graves p.s. If you try it and have trouble, please submit another question including commented, minimal, self-contained, reproducible code, as requested

Re: [R] Interesting remarks about R back in 1999

Many people love Mathematica, but it's strength is symbolic mathematics, not data analysis. Googling for data analysis in Mathematica and R led me to an advertisement for a Mathematica add-on called RLink, which is an exciting new tool [to] leverage the statistical analysis power of R

Re: [R] Irregular Time Series Issue

Have you considered using 'corCAR1' with 'lme' in the 'nlme' package? Hope this helps. Spencer Graves p.s. If you are not familiar with this, I highly recommend Pinheiro, J.C., and Bates, D.M. (2000) Mixed-Effects Models in S and S-PLUS (Springer). The ~R\library\nlme

Re: [R] Overloading %*%

again for your reply regarding %*%. Spencer Graves knoblauch wrote: Joe Cainey jcainey at gmail.com writes: Is it possible to supply a new method for the %*% operator? clipped I've tried to do the same thing with %*%: %*%.ad - function(a,b) { # further code

Re: [R] Overloading %*%

Dear Martin: This is wonderful. Thank you very much. It would be a great help if your suggestions could be added to See Also and Examples for methods. Thanks again, Spencer Graves Martin Morgan wrote: Spencer Graves wrote: Thanks, Ken. 1. How can I

[R] function for the average or expected range?

- (F(x))^n - (1-F(x))2] dx. Thanks, Spencer Graves __ R-help@r-project.org mailing list https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read the posting guide http://www.R-project.org/posting-guide.html and provide commented, minimal

[R] function for the average or expected range?; CORECTION

: The exact distribution of the range, p. 339). Unfortunately, either I don't understand this formula, or it's wrong. Using expression (14.1) in the same reference, I get the following: E(R) = n*integral{-Inf to Inf of [(F(x))**(n-1) - (1-F(x))**(n-1)]dF(x). Thanks, Spencer Graves

[R] function for the average or expected range?; CORECTION TO CORRECTION

). Thanks, Spencer Graves __ R-help@r-project.org mailing list https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read the posting guide http://www.R-project.org/posting-guide.html and provide commented, minimal, self-contained, reproducible

Re: [R] function for the average or expected range?; CORECTION

Hi, Greg: Thanks very much for the reply. 1. The 'ptukey' and 'qtukey' function are the distribution of the studentized range, not the range. I tried sum(ptukey(x, 2, df=Inf, lower=FALSE))*.1 and got 1.179 vs. 1.128 in the standard table of d2 for n = 2 observations per

Re: [R] function for the average or expected range?; CORECTION

. *From:* Spencer Graves [mailto:[EMAIL PROTECTED] *Sent:* Fri 3/21/2008 3:39 PM *To:* Greg Snow *Cc:* r-help@r-project.org *Subject:* Re: [R] function for the average or expected range?; CORECTION Hi, Greg: Thanks very much for the reply. 1. The 'ptukey' and 'qtukey

Re: [R] scaling problems in optim

' -- while perhaps illustrating the need to use these parameters. Hope this helps. Spencer Graves library(mvtnorm) (mle2 - optim(rep(0, 2), dmvnorm, method='CG', control=list(fnscale=-10, parscale=c(.3, 3), trace=9), mean=1:2, hessian=TRUE, log=TRUE

Re: [R] Simulate ARX model.

Have you tried the vignette: (dse - vignette(dse-guide)) Stangle(dse$file) This should open dse-guide.pdf in Acrobat and create a file dse-guide.R in the working directory. The latter file should contain the R commands used in the vignette. Do you

Re: [R] dreaded p-val for d^2 of a glm / gam

variable, collect anova(glm.D93t, glm.D93)[2, Deviance] in a vector, and then find out how extreme the deviance you actually got is relative to this permutation distribution. Hope this helps. Spencer Graves p.s. Regarding your 'dread', please see fortune(children) Monica Pisica

Re: [R] dreaded p-val for d^2 of a glm / gam

be misused. Hope this helps. Spencer Graves Monica Pisica wrote: Hi Spencer and David, Thanks for your answers first - yes it is deviance but just before i just spoke and explain that it is the equivalent of r square from the normal regression. I hope i can do

Re: [R] operators %/% bug?

Consider also the following: a - 1.2 b - 1.0 (a-b)%%.1 [1] 0.1 (a-b)%/%.1 [1] 1 .2%%.1 [1] 0 .2%/%.1 [1] 2 This performs as described on ?%/%. In particular, the second paragraph under Value: '%%' indicates 'x mod y' and '%/%' indicates integer division. It is

Re: [R] optimize simultaneously two binomials inequalities using nlm( ) or optim( )

. Spencer Graves [EMAIL PROTECTED] wrote: Dear R users, I´m trying to optimize simultaneously two binomials inequalities (used in acceptance sampling) which are nonlinear solution, so there is no simple direct solution. Please, let me explain shortly the the problem and the question as following

Re: [R] model mix problem. FALSE CONVERGENCE

replies quicker. ¡Buena Suerte! Spencer Graves mergullo wrote: Hi R users I try to use the lme but I can´t! My script is (some words in french, sorry!!): rm(list=ls(all=TRUE)) #Efface tous les objets en mémoire pour éviter des erreurs library(MASS) #Chargement des Librairies

Re: [R] predict for garch

Hi, Erin: I just got 46 hits from RSiteSearch('predict garch'). Might any of these help you? Spencer Erin Hodgess wrote: Dear R People: Is there a predict method for garch methods, please? I tried the usual predict(d1.garch,n.ahead=3) but only got values for the original

Re: [R] how to read multiple lines per case

Dear Michael: There are doubtless better ways to do this, but the I just got the following ugly hack to work for me: rawDat - readLines('linthall.dat') (N - length(rawDat)) N2 - ((N-1)/2) outDat - character(N2+1) outDat[1] - rawDat[1] for(i in 1:N2) outDat[i+1] -

[R] ARMA random effects?

a different set for each Mare in this example? Thanks, Spencer Graves __ R-help@r-project.org mailing list https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read the posting guide http://www.R-project.org/posting-guide.html and provide

Re: [R] False convergence in LME

from my study of it. Hope this helps. Spencer Graves Rebecca Sela wrote: I tried to use LME (on a fairly large dataset, so I am not including it), and I got this error message: Error in lme.formula(formula(paste(c(toString(TargetName), as.factor(nodeInd)), : nlminb problem

Re: [R] Book on R graphics

' graphics from self-study of a book. Hope this helps. Spencer Graves Roberto Laforgia wrote: See the following link: http://www.amazon.com/Graphics-Computer-Science-Data-Analysis/dp/158488486X/ref=wl_itt_dp?ie=UTF8coliid=I3TIYS7LRDPIEDcolid=1ZQCB91VGW5UR Regards, R.L. Marcin Kozak

Re: [R] constrOptim with method = L-BFGS-B

I believe that 'optim' will not accept equality constraints. However, you do not need the generality of 'optim' to minimize a quadratic function with boundary conditions and one equality condition. This type of problem is called quadratic programming, and RSiteSearch(quadratic

Re: [R] Differential Equations

I just got 37 hits from RSiteSearch('differential equation', 'fun') including odesolv{fda}, rk4{deSolve}, rk4{odesolve}, in addition to lsoda{odesolve} and the Rsundials package. hope this helps. Spencer Graves Michael Lawrence wrote: On Tue, Jun 17, 2008 at 6:25 PM, David

Re: [R] Complex Time Series

If this were my problem, I think I'd express it as a Kalman filter and use the 'dlm' package. Fortunately, the package contains a vignette. Hope this helps. Spencer Graves Bernardo Rangel Tura wrote: Hi R masters, In my work I analyse a time serie of number of birth

Re: [R] Comparison between R and MATLAB

as I seen. Spencer Graves stephen sefick wrote: This discussion has already occurred- to my knowledge at least once. I would suggest searching the archived list, and see what you get. If you still have questions after you have a look then fire a couple of specifics, but speaking for myself I

Re: [R] Programming Concepts and Philosophy

Maybe I missed it, but I haven't seen in this tread mention of Venables and Ripley (2000) S Programming (Springer). I found it interesting and useful, though I have not used it as much as MASS -- partly because it is more specialized and it's coverage is not as broad. Spencer

Re: [R] Hessian in box-constraint problem - concern OPTIM function

in the posting guide http://www.R-project.org/posting-guide.html;. Hope this helps. Spencer Graves Cleber Nogueira Borges wrote: Hello all useRs, I am using the OPTIM function with particular interest in the method L-BFGS-B, because it is a box-constraint method. I have interest

Re: [R] running R-code outside of R

believe in magic, but sometimes it requires a little creative thought. Hope this helps. Spencer Graves Rolf Turner wrote: On 25/06/2008, at 6:01 AM, Roger Leenders wrote: R 2.7.0, WinXP Hi list, I have written some code in R that I would like to share with others who

Re: [R] running R-code outside of R

If you want to hide the fact that you are using R -- especially if you charge people for your software that uses R clandestinely -- that's a violation of the license (GPL). I doubt if anyone associated with R would bother with a lawsuit, but a competitor who offers related software

Re: [R] Hessian in box-constraint problem - concern OPTIM function

in S and S-Plus (Springer). Hope this helps. Spencer Graves Cleber Nogueira Borges wrote: Hello all useRs, I am using the OPTIM function with particular interest in the method L-BFGS-B, because it is a box-constraint method. I have interest in the errors estimates too. I make

Re: [R] Test for multiple comparisons: Nonlinear model, autocorrelation?

The question seems too general for me to offer specific suggestions. What problem are you trying to solve that you think 'multiple comparisons' will answer? Can you produce a similar problem that is completely self-contained example that eliminates complexity that may not be needed to

Re: [R] Test for multiple comparisons: Nonlinear model, autocorrelation?

If I had only a very limited time to do this, I might include 'month' as another effect, probably coded as 'sin' and 'cos' on an annual cycle rather than as 12 individual Indicators. This would allow you to explore not only main effects but interactions with plots. Before I did

[R] [R-pkgs] Functional Data Analysis, fda_1.2.4

', with methods for 'plot', 'lines', and 'plotfit'. In addition, there have also been improvements to the 'plotfit' function and the 'CanadianWeather' data. Also, a function 'norder' has been added, which returns the order of a B-spline. Spencer Graves

Re: [R] Operator overloading

The 'polynom' package provides one example of how to do this. However, getAnywhere(+.polynomial) just told me no object named ‘+.polynomial’ was found; I don't know why. The documentation says it uses S3 classes. The 'Matrix' package should provide other examples, using S4 classes.

Re: [R] R bug in the update of nlme?

Did you try 'RSiteSearch(unable to restore saved data in .RData)'? This just produced 61 hits for me, the second of which http://finzi.psych.upenn.edu/R/Rhelp02a/archive/85451.html; looks to me like it might answer your question. Hope this helps. Spencer Susana

  1   2   3   4   5   6   >