Re: [R] getting confidence intervals after lm fitting

2004-01-23 Thread Prof Brian Ripley
On Thu, 22 Jan 2004, Yun-Fang Juan wrote: After fitting with lm, I want to get the 95 % confidence intervals of the expected predicted value. Of the predicted expected value, I assume? I know in S-plus I can use pointwise() to get the confidence intervals . But in R, I couldn't find such a

Re: [R] BVLS

2004-01-23 Thread Prof Brian Ripley
On Fri, 23 Jan 2004, Briggs, Meredith M wrote: Is there an R package that solves linear least squares with upper and lower bounds on the variables. Something like the Parker and Stark algorithm written in Fortran. Probably not, but this is an easy application of optim(method=L-BFGS) or of

Re: [R] Syntay-Highlighting for KDE-Kate

2004-01-23 Thread Arne Henningsen
Hi, I am the new maintainer of the R syntax highlighting plugin for kate. (The initial author, Egon, stopped maintaining it, because he does not use kate anymore.) @ Andy: Thanks for providing the link to this file. (I was very busy this week and did not look at the R mailing list until

Re: [R] Turnbull estimate

2004-01-23 Thread Martyn Plummer
On Fri, 2004-01-23 at 01:33, Denise wrote: Hi everyone, I've been looking for a function that calculates the Turnbull estimate for left, right and interval censored data. None of the data that I am using has exact failure times. The only function I seem to find can handle data that has no

RE: [R] Fitting compartmental model with nls and lsoda?

2004-01-23 Thread Jesus Frias
Now, I am including Richar Upton's 2 cm model implementation and Christoffer Tornoe's nls solution (I recommend Christoffer's nlmeODE package for these problems also if multi-response data is available) The code follows: Multi-response data?, Is there a way of dealing with those

[R] Question

2004-01-23 Thread Wassim Kamoum
Hello I' d like to implement a space-time model as it is formulated Technometrics vol 22 1980 A three stage iterative procedure for space-time modelling P,E Pfeiffer, and S.J Deutsch For model identification the sample space-time autocorrelation and partial autocorrelation function are

[R] (no subject)

2004-01-23 Thread Tamara Porter
I downloaded R onto my computer using Windows XP The program itself works fine, but my problem occurs when I type help.start() and enter the search engine. Whenever I type anything into it, nothing happens. I've tried reinstalling it a few times, but the same thing keeps

[R] predict.lda problem with posterior probabilities

2004-01-23 Thread Mike White
With predict.lda the posterior probabilities only relate to the existing Class definitions. This is fine for Class definitions like gender but it is a problem when new data does not necessarily belong to an existing Class. Is there a classification method that gives posterior probabilities for

Re: [R] predict.lda problem with posterior probabilities

2004-01-23 Thread Marc R. Feldesman
At 04:23 AM 1/23/2004, Mike White wrote: With predict.lda the posterior probabilities only relate to the existing Class definitions. This is fine for Class definitions like gender but it is a problem when new data does not necessarily belong to an existing Class. Is there a classification method

[R] Levels number of a factor object

2004-01-23 Thread Laurent Houdusse
Hi all! How to retrieve the levels number of a factor object? See this code: groups-gl(4,10) I want to retrieve the number of levels (4) of my object groups I tried groups.levels but this don't work Thanks Laurent Houdusse Analyste Programmeur

[R] question

2004-01-23 Thread Wassim Kamoum
Hello I was looking for the source code for the new arima procedure the STARMA model (Space-Time ARMA) If somebody have the source code of this model please send me them I'm student and in my research for my master I'm appling the STARMA model for modelling the pollutant particules. It's

Re: [R] how to take derivatives of a step function

2004-01-23 Thread Eugene Salinas (R)
Thanks. Here is some more info that may help: I am estimating a transformation model H(Y)=Xb+e and have obtained (up to a location parameter) an estimated of H(T), call it H_est(T). (Its a new estimator I am trying out that works with unspecified functional forms etc). Now if the model is a

Re: [R] Levels number of a factor object

2004-01-23 Thread Marc Schwartz
On Fri, 2004-01-23 at 08:24, Laurent Houdusse wrote: Hi all! How to retrieve the levels number of a factor object? See this code: groups-gl(4,10) I want to retrieve the number of levels (4) of my object groups I tried groups.levels but this don't work Thanks

Re: [R] Levels number of a factor object

2004-01-23 Thread Thomas Petzoldt
Laurent Houdusse schrieb: Hi all! How to retrieve the levels number of a factor object? See this code: groups-gl(4,10) I want to retrieve the number of levels (4) of my object groups I tried groups.levels but this don't work Simply try: length(levels(groups)) Thomas P.

[R] Response Surface

2004-01-23 Thread Li Xiaolei
Hi, Is there any existing way in R of doing response surface analyses and plotting the response surface, like what minitab can do? Thanks Xiaolei __ [EMAIL PROTECTED] mailing list https://www.stat.math.ethz.ch/mailman/listinfo/r-help PLEASE do read

Re: [R] question

2004-01-23 Thread Peter Dalgaard
Wassim Kamoum [EMAIL PROTECTED] writes: Hello I was looking for the source code for the new arima procedure the STARMA model (Space-Time ARMA) If somebody have the source code of this model please send me them I'm student and in my research for my master I'm appling the STARMA

[R] STARMA model (Space-Time ARMA)

2004-01-23 Thread Wassim Kamoum
Hello I was looking for the source code for the new arima procedure the STARMA model (Space-Time ARMA) If somebody have the source code of this model please send me them I'm student and in my research for my master I'm appling the STARMA model for modelling the pollutant particules. It's

Re: [R] how to take derivatives of a step function

2004-01-23 Thread Peter Dalgaard
Eugene Salinas (R) [EMAIL PROTECTED] writes: Thanks. Here is some more info that may help: I am estimating a transformation model H(Y)=Xb+e and have obtained (up to a location parameter) an estimated of H(T), call it H_est(T). (Its a new estimator I am trying out that works with unspecified

Re: [R] Axes Ticks

2004-01-23 Thread ivo welch
thanks to everyone for their help. easy now. best regards, /iaw __ [EMAIL PROTECTED] mailing list https://www.stat.math.ethz.ch/mailman/listinfo/r-help PLEASE do read the posting guide! http://www.R-project.org/posting-guide.html

Re: [R] question - STARMA model (Space-Time ARMA)

2004-01-23 Thread Roger Bivand
On 23 Jan 2004, Peter Dalgaard wrote: Wassim Kamoum [EMAIL PROTECTED] writes: Hello I was looking for the source code for the new arima procedure the STARMA model (Space-Time ARMA) If somebody have the source code of this model please send me them I'm student and in my

[R] Problem with hasArg() using R-files

2004-01-23 Thread Siegfried.Macho
Dear R-helpers, I got the following problem with the function hasArg() I would like to use the function to check whether the argument crit is present as an argument of the function (here is the simplified version of the function): SDT.Optim - function(crit = NULL, Hess = F) {

Betr.: [R] Levels number of a factor object

2004-01-23 Thread tobias . verbeke
[EMAIL PROTECTED] wrote on 23/01/2004 15:24:21: Hi all! How to retrieve the levels number of a factor object? See this code: groups-gl(4,10) I want to retrieve the number of levels (4) of my object groups I tried groups.levels but this don't work Is

RE: [R] Problem with hasArg() using R-files

2004-01-23 Thread Liaw, Andy
Please do give reproducible example. The code you gave, which you claimed `works correctly' doesn't: SDT.Optim - function(crit = NULL, Hess = F) + { + q - length(par); x - data + if(hasArg(crit)) + cat(\n Crit

Re: [R] Levels number of a factor object

2004-01-23 Thread Spencer Graves
groups-gl(4,10) levels(groups) [1] 1 2 3 4 length(levels(groups)) [1] 4 Is this what you want? spencer graves Laurent Houdusse wrote: Hi all! How to retrieve the levels number of a factor object? See this code: groups-gl(4,10) I want to retrieve the number of levels

[R] graphical book suggestion

2004-01-23 Thread ivo welch
as I am playing around with parameters to graphics, it seems to me that there would be a market for a book written that is the equivalent of the perl cookbook, but for graphics creation in R/S. just demonstrating how each parameter works, how one can make complex graphics, some graphics

Re: [R] Levels number of a factor object

2004-01-23 Thread Peter Dalgaard
Laurent Houdusse [EMAIL PROTECTED] writes: Hi all! How to retrieve the levels number of a factor object? See this code: groups-gl(4,10) I want to retrieve the number of levels (4) of my object groups I tried groups.levels but this don't work Thanks Do you mean

Re: [R] (no subject): really about search in HTML help under Windows

2004-01-23 Thread Prof Brian Ripley
See the rw-FAQ, which explains this (as you have already been told privately). The problem is in your browser configuration, not in R. You need Java and Javascript installed and enabled. On Fri, 23 Jan 2004, Tamara Porter wrote: I downloaded R onto my computer using Windows XP The

[R] correlation between two time serials

2004-01-23 Thread Jieping
Hi, there, I'd like to get a statistics to reflect the correlation of two time serials observations. for example, one time serials over time 1 to 20, and another one from same time points from 1 to 20. I need one scalar satistics as a correlation of these two observations. Thank you and have a

[R] Fast ARMA

2004-01-23 Thread Andrew Robinson
Greetings R-friends, I'm looking for some code to fit an ARMA, up to (3,3), preferably with ML. I'm doing a simulation experiment, so I need to fit it to the order of millions of times. I'm using Professor Ripley's arima() command right now, but I hope to find soemthing faster in c. Does anyone

[R] lags in regressions

2004-01-23 Thread Tobias Mühlhofer
Hi! I am trying to get R to run regressions for me of a variable on lagged differences of itself. Specifically: x-x(-1) = a + b1(x(-1)-rx(-2))+b2(x(-4)-rx(-5))+e I need to do this a lot of times, altering the value of r. What I've been trying to do was to use the lag() command to create

Re: [R] lags in regressions

2004-01-23 Thread Prof Brian Ripley
You need to do tsunion (or something like it) on the time series. lm() does not know about time series and you are merely supplying vectors as far as it is concerned. On Fri, 23 Jan 2004, Tobias Mühlhofer wrote: Hi! I am trying to get R to run regressions for me of a variable on lagged

[R] Re: matrix __power__ (was exponential)

2004-01-23 Thread VICENTED.CANTO.EXT
Sorry I didn't answer before. Martin, thanks very much for your notes and for tidying the function up!= ! I've learnt a lot with you. About matPower, I didn't realise that where it says pot - pot %*% p ot, consumption of memory was higer than it was on the original algorithm,=

[R] Re: nonlinear regression and Excel solver

2004-01-23 Thread Kristian Omland
Hi all, Here is a summary of substantive replies to my posts (on both s-news and r-help) regarding use of Excel's Solver for a nonlinear regression problem. A number of people replied that Solver performs well as an optimizer based on experience, particularly if given reasonable starting

RE: [R] Fitting compartmental model with nls and lsoda?

2004-01-23 Thread DivineSAAM
Dear Jesus, Yes there is a way and it is via Christoffer Torn\{o}e's package nlmeODE. I checked Chapter 4 of Bates and Watts. As usual, a little work involved, but doable and quite powerful. Thanks, olinares __ [EMAIL PROTECTED] mailing list

[R] ylab argument to plot.stepfun

2004-01-23 Thread Aleksey Naumov
Dear R Users, Does anyone know how to specify Y-label in plot.stepfun()? I do: library(stepfun) x = ecdf(rnorm(30)) plot(x, ylab='CDF') and I get: Error in plot.stepfun(ylab = Fn(x), ..., verticals = verticals) : formal argument ylab matched by multiple actual arguments Thanks, Aleksey

[R] regression to get non-negative parameters

2004-01-23 Thread Yun-Fang Juan
Hi, I am working on a model sort of like a credit assignment problem. The goal is to get non-negative weights for each input variable and discount the predictors (down to 0) that lead to no credits. I am trying out the lm fitting and find out they give a lot of negative weights. I am totally

Re: [R] (no subject)

2004-01-23 Thread Duncan Murdoch
On Fri, 23 Jan 2004 12:08:31 +, you wrote: I downloaded R onto my computer using Windows XP The program itself works fine, but my problem occurs when I type help.start() and enter the search engine. Whenever I type anything into it, nothing happens. I've tried reinstalling