Re: [R] Internal NA removal out of Time Series with na.omit.ts()

2004-03-08 Thread Adrian Trapletti
Hi R specialists, The na.omit.ts() method fails when the time series contains internal NA's. How can these automatically be removed? try na.remove from tseries. This is, e.g., useful when removing weekends (NA prices) from financial data, i.e., switching from physical time to business time.

Re: [R] Internal NA removal out of Time Series with na.omit.ts()

2004-03-08 Thread Adrian Trapletti
Prof Brian Ripley wrote: On Mon, 8 Mar 2004, Adrian Trapletti wrote: The na.omit.ts() method fails when the time series contains internal NA's. How can these automatically be removed? try na.remove from tseries. This is, e.g., useful when removing weekends (NA prices) from financial

[R] lm - significance disappears

2004-03-08 Thread Stuart Leask
Hi there. I am trying to model the effect of Age at onset of a condition (AGEONSET) upon IQ (ie. do they have more problems if they get ill younger?). I also want to see if there is an interaction with the TYPE of test (reading maths). I am struggling to understand why the effect of AGEONSET, so

[R] getting the std errors in the lm function

2004-03-08 Thread Fulvio Copex
Hello, I have a simple question for you: making: mylm-lm(y~x) summary(mylm) I get the following results: ** Coefficients: Estimate Std. Error t value Pr(|t|) (Intercept) 16.540870.19952 82.91 2e-16 *** x[1:19] -2.32337

Re: [R] getting the std errors in the lm function

2004-03-08 Thread Prof Brian Ripley
On Mon, 8 Mar 2004, Fulvio Copex wrote: Hello, I have a simple question for you: making: mylm-lm(y~x) summary(mylm) I get the following results: ** Coefficients: Estimate Std. Error t value Pr(|t|) (Intercept) 16.54087

[R] dsn

2004-03-08 Thread Margarida Júlia Rodrigues Igreja
Hi, I have a data base in oracle and need to link R using ODBC. When i use the command odbcConnect a dsn is needed. I would like to create a dsn.Do you know how can i do it? Thank you. Margarida,Portugal __ [EMAIL PROTECTED] mailing list

Re: [R] lm - significance disappears

2004-03-08 Thread Prof Brian Ripley
On Mon, 8 Mar 2004, Stuart Leask wrote: Hi there. I am trying to model the effect of Age at onset of a condition (AGEONSET) upon IQ (ie. do they have more problems if they get ill younger?). I also want to see if there is an interaction with the TYPE of test (reading maths). I am

RE: [R] dsn

2004-03-08 Thread Wayne Jones
Hi If you are using windows then do this: Basically you need to add your oracle database as an ODBC data source. Control Panel - Administrative tools - Data Sources (ODBC) Select the System DSN tab and then Add. Select which type of Data source i.e. Oracle database. Then follow the

Re: [R] getting the std errors in the lm function

2004-03-08 Thread Peter Dalgaard
Fulvio Copex [EMAIL PROTECTED] writes: Hello, I have a simple question for you: making: mylm-lm(y~x) summary(mylm) I get the following results: ** Coefficients: Estimate Std. Error t value Pr(|t|) (Intercept) 16.54087

[R] Graphical Test

2004-03-08 Thread JVP L
Dear R mailing list, Greetings!!! I'm doing a graphical test of suitability for the aircondit data. I actually just followed the exercise in Davison and Hinkley's book of Graphical Test for model checking. The program for testing if it comes from an exponential distribution seems to be

RE: [R] boot package

2004-03-08 Thread Hanke, Alex
Hi, The function anosim() in vegan package or sample() in base may be of help to you. Alex -Original Message- From: Rogério Rosa da Silva [mailto:[EMAIL PROTECTED] Sent: March 4, 2004 4:57 AM To: rhelp Subject: [R] boot package Dear all As part of an ongoing study on the

[R] a question on CSV file

2004-03-08 Thread Cindy Juncker
I have a csv file, heartrate excercise medicine Time 12 00 4:30 am 23 5 0 5:00 am 34 6 0 5:30 am

Re: [R] applying data generating function

2004-03-08 Thread Spencer Graves
With gc() right before each call to proc.time, as Brian Ripley and Gabor Grothendieck suggested, the times were substantially more stable. For the for loop, extending the vector with each of 1e5 iterations, I got 181.25, 181.27, 182.72, 182.44, and 182.56. The averages of the last 3 of

RE: [R] Statistical Quality Control

2004-03-08 Thread apjaworski
A new package named qcc just appeared on CRAN. It has standard QC charts, OC curves and process capability stuff. It even has a Minitab-like process capability sixpack. It looks like it is pretty nicely done. Andy __ Andy Jaworski 518-1-01 Process

[R] years from as.POSIXlt

2004-03-08 Thread Christian Schulz
Hi, how it's possible to extract the year and the number of days from Julian date. i'm little confused about the last two functions and ?years . EDATE comes from sqlQuery with as.is=T EDATE - as.POSIXlt(datvears$ENROLLDAY) Many thanks, Christian EDATE[1:5] [1] 2000-06-30 11:25:01

Re: [R] Internal NA removal out of Time Series with na.omit.ts()

2004-03-08 Thread Jan Verbesselt
Thanks for the tips and advice! I found out how the NA's came into the time series. By a threshold mechanism, extreme values (outliers) are removed from the time series (environmental data, remote sensing data) and set as NA. A solution could be to detect an outlier and replace it by the

[R] question

2004-03-08 Thread cfinet
I do not manage to make a Fisher´s exact test with the next matrix : [,1] [,2] [,3] [,4] [,5] [,6] [,7] [,8] [,9] [,10] [,11] [,12] [1,]130129002 5 8 6 [2,]033000050 3 0 0 [3,]0000

[R] abline

2004-03-08 Thread Grace Conlon
if I want to specify y-coordinates for the heights of horizontal lines to go across a plot. x - c(1,2,3,6,4,8,4,7) y - c(3,2,7,4,5,4,5,6) h - c(3,5,7) plot(x,y) abline(y=h) However I got error message: Warning message: parameter y couldn't be set in high-level plot() function (I tried

Re: [R] unusual name in .Rd file: documenting an infix function

2004-03-08 Thread David Brahm
Christian Hoffmann [EMAIL PROTECTED] wrote: I am having difficulties documenting an infix function: %% - function(x,y) { paste(x,y,sep=) } We have the exact same function, and the following documentation works fine. Note the documentation filename is g.am.Rd, but that's just a random name,

Re: [R] abline

2004-03-08 Thread Spencer Graves
?abline in R 1.8.1 includes the following: Usage: abline(a, b, untf = FALSE, ...) abline(h=, untf = FALSE, ...) abline(v=, untf = FALSE, ...) abline(coef=, untf = FALSE, ...) abline(reg=, untf = FALSE, ...) It looks to me like you want abline(h=h). hope this helps.

Re: [R] abline

2004-03-08 Thread XIAO LIU
x - c(1,2,3,6,4,8,4,7) y - c(3,2,7,4,5,4,5,6) h - c(3,5,7) plot(x,y) abline(y=h) abline(h=h) will work. __ [EMAIL PROTECTED] mailing list https://www.stat.math.ethz.ch/mailman/listinfo/r-help PLEASE do read the posting guide!

Re: [R] graphic device MetaPost

2004-03-08 Thread Na Li
On 7 Mar 2004, Jinsong Zhao said: use Chinese character in plots with a minor modification to the MetaPost file. I never tried this, but maybe it is possible to use the psfrag LaTeX package to insert Chinese characters into an eps file (with the help of proper ps fonts and the CJK packages of

[R] Creating CART: A Personal Retrospective, Opening Session of CART Data Mining 2004

2004-03-08 Thread Lisa Solomon
OPENING SESSION OF CART DATA MINING 2004, Monday March 22nd (no charge) Welcoming Reception: 5:30 - 7:30 PM Session Begins promptly at 7:30 PM Personal recollections of Leo Breiman, Jerome Friedman, Richard Olshen, Charles Stone. Join us to hear Breiman, Friedman, Olshen and Stone discuss their

RE: [R] abline

2004-03-08 Thread Gabor Grothendieck
abline does not have a y parameter so it is letting you know that you have specified an invalid parameter. See ?abline Issuing the command args(abline) is another way to find out the arguments to a command. At any rate, based on your description I think you want segments. See ?segments

[R] Decision Trees

2004-03-08 Thread Benjamin . STABLER
I am familiar with the rpart and tree packages for classification and regression trees. However, quite a bit of the research in the transportation community relating to decision trees uses the C4.5 family of algorithms by Quinlan. Are there any plans to make a C4.5 (or a derivative of it)

Re: [R] normal scores test

2004-03-08 Thread Torsten Hothorn
On Sat, 6 Mar 2004, Padmanabhan, Sudharsha wrote: Hello, I need help in performing a Van_der_Waerden normal scores test in R. I have two arrays of scores(final on therapy scores from drug and placebo) and want to use the normal scores procdeure to test for significance. (observations are

RE: [R] a question on CSV file

2004-03-08 Thread Gabor Grothendieck
Note that a CSV file is one in which data elements are separated by commas. There are no commas in what you display so I assume that you are displaying the data after you have read it in, not the CSV file. If you read a file into a data frame x you can refer to row i column j of the data frame

[R] test

2004-03-08 Thread apjaworski
I am trying to test whether I can post to the list. I responded to a couple of messages recently and I have not seen my response posted. __ Andy Jaworski 518-1-01 Process Laboratory 3M Corporate Research Laboratory - E-mail: [EMAIL PROTECTED] Tel: (651)

Re: [R] memory problem

2004-03-08 Thread Peter Dalgaard
Joshi, Nina (NIH/NCI) [EMAIL PROTECTED] writes: I am trying to upload into R 143 Affymetrix chips onto using R on the NIH Nimbus server. I can load 10 chips without a problem, however, when I try to load 143 I receive a error message: cannot create a vector of 523263 KB. I have expanded the

[R] Different missing links on Windows in 'check' vs. 'install'

2004-03-08 Thread Jeff Gentry
Hello ... Using R-1.9.0 alpha, I'm having some problem getting a few packages to pass check under Windows - specifically with the 'missing link(s)' section of the package install phase. I started trying to track down how the missing link was showing up as in some cases I could not see why the

Re: [R] Internal NA removal out of Time Series with na.omit.ts()

2004-03-08 Thread kjetil
On 8 Mar 2004 at 15:29, Jan Verbesselt wrote: Thanks for the tips and advice! I found out how the NA's came into the time series. By a threshold mechanism, extreme values (outliers) are removed from the time series (environmental data, remote sensing data) and set as NA. A solution could

RE: [R] a question on CSV file

2004-03-08 Thread Gabor Grothendieck
In rereading what I had written there is an error. The first statement selecting Med as 1 should be as follows: x[x$Med == 1,] --- Date: Mon, 8 Mar 2004 20:53:48 -0500 (EST) From: Gabor Grothendieck [EMAIL PROTECTED] To: [EMAIL PROTECTED], [EMAIL PROTECTED] Subject: RE: [R] a question

RE: [R] applying data generating function

2004-03-08 Thread Liaw, Andy
-Original Message- From: [EMAIL PROTECTED] [mailto:[EMAIL PROTECTED] On Behalf Of Spencer Graves Sent: Monday, March 08, 2004 11:21 AM To: Prof Brian Ripley Cc: Fred J.; r-help; Peter Dalgaard Subject: Re: [R] applying data generating function With gc() right before

RE: [R] abline

2004-03-08 Thread Liaw, Andy
So exactly what do you want? h - c(3,5,7) abline(h=h) draws three horizontal lines at y=3, 5, and 7, as the help page for abline says it should. [Note: that's abline(h=h), not abline(y=h).] Andy From: Grace Conlon if I want to specify y-coordinates for the heights of horizontal lines

[R] About reading data into R

2004-03-08 Thread Rui Song
I have a problem about reading data into R. There is a \n between each pair of data, like: -155.65 -155.77 -155.40 -155.46 -155.52 -155.34 ... Could anyone tell me how to read in such data? Thanks! Rui

[R] Re: R-help Digest, Vol 13, Issue 8

2004-03-08 Thread Yong Wang
On Mon, 8 Mar 2004 [EMAIL PROTECTED] wrote: Send R-help mailing list submissions to [EMAIL PROTECTED] To subscribe or unsubscribe via the World Wide Web, visit https://www.stat.math.ethz.ch/mailman/listinfo/r-help or, via email, send a message with subject or body 'help' to

[R] how to replace NaN in a vector

2004-03-08 Thread Yong Wang
Hi, all a vector such as (1,2,4,-1,NaN,2,4,NaN) if try to replace all the NaN with a numerical value, what's the easiest way? thanks a lot best yong __ [EMAIL PROTECTED] mailing list https://www.stat.math.ethz.ch/mailman/listinfo/r-help PLEASE do

[R] Content violation

2004-03-08 Thread Symantec Alerts
The email message: From: [EMAIL PROTECTED] To: [EMAIL PROTECTED] includes the file your_document.pif which is prohibited by Bechtel Information Security policy. The file was blocked and did not reach the intended recipient. If

[R] vector extraction

2004-03-08 Thread Jens Hainmueller
Hello, I could need some help on this one: From the data.frame Test.dataset2 below (TSCS data for 151 countries.to.map for year 1973-95; each country.to.map is described by a unique code), I would like to extract a vector color that for each country.to.map takes on the value of dv (a categorical

Re: [R] how to replace NaN in a vector

2004-03-08 Thread Duncan Murdoch
On Mon, 8 Mar 2004 18:18:28 -0600 (CST), you wrote: Hi, all a vector such as (1,2,4,-1,NaN,2,4,NaN) if try to replace all the NaN with a numerical value, what's the easiest way? x[is.nan(x)] - 42 See ?NaN for related functions. Duncan Murdoch __

Re: [R] how to replace NaN in a vector

2004-03-08 Thread David Scott
On Mon, 8 Mar 2004, Yong Wang wrote: Hi, all a vector such as (1,2,4,-1,NaN,2,4,NaN) if try to replace all the NaN with a numerical value, what's the easiest way? thanks a lot Is this what you wanted? x-c(1,2,4,-1,NaN,2,4,NaN) y-ifelse(is.nan(x),-999,x) y [1]124

[R] SSfpl question

2004-03-08 Thread Yun-Fang Juan
Hi, I am trying to do a nls fitting using the SSfpl self-restart function but I got an error like the following: nls.2 - nls(map.rr ~ SSfpl(mail.pv, phi1, phi2, phi3, phi4), data=pv.rr.200) Error in nls(y ~ cbind(1, 1/(1 + exp((xmid - x)/exp(lscal, data = xy, : step factor

[R] More factor names in x axis, how?

2004-03-08 Thread Richard A. O'Keefe
I have some data which records, amongst other things, age (recoded by me to be in months), and drug group (a factor). The drug group is a 2 digit number, but there is no numeric relationship, and only 50 of the possible groups occur So asd - read.table(asd.dat, header=TRUE) asd$group -

Re: [R] Different missing links on Windows in 'check' vs. 'install'

2004-03-08 Thread Duncan Murdoch
On Mon, 8 Mar 2004 18:20:29 -0500 (EST), you wrote: In trying to figure out how these things are determiend, I noticed that in Windows there was a difference in the reported missing links when one does 'Rcmd check' vs 'Rcmd install' (or 'Rcmd install --build'). In the example above, using either

Re: [R] Monte Carlo p-value (was question)

2004-03-08 Thread Spencer Graves
What is the standard convention for a Monte Carlo p-value when the observed outcome is more extreme than all simulations? The example provided by Cédric Fine produced a Monte Carlo p-value, according to Kjetil Halvorsen, of 2.2e-16 (based on 2000 replicates). This seems inappropriate to

[R] Am failing on making lagged residual after regression

2004-03-08 Thread Ajay Shah
Folks, I'm most confused in trying to do something that (I thought) out to be mainstream and straightforward R. :-) Could you please help? I am doing an ordinary linear regression. My goal is: After a regression, to make residuals, and make a new variable which is the lagged residuals (lagged by

RE: [R] vector extraction

2004-03-08 Thread Gabor Grothendieck
My understanding is that you want to extract the dv values for year 1980: color - with( Test.dataset2, dv[ year == 1980 ] ) --- Date: Mon, 8 Mar 2004 19:36:45 -0800 From: Jens Hainmueller [EMAIL PROTECTED] To: R-Help [EMAIL PROTECTED] Subject: [R] vector extraction Hello,

RE: [R] More factor names in x axis, how?

2004-03-08 Thread Gabor Grothendieck
You can omit the x axis from plot and create it yourself using axis: plot(age ~ group, data=asd, xaxt=n) axis(1, 1:nlevels(group), levels(group)) --- Date: Tue, 9 Mar 2004 17:09:23 +1300 (NZDT) From: Richard A. O'Keefe [EMAIL PROTECTED] To: [EMAIL PROTECTED] Subject: [R] More

RE: [R] About reading data into R

2004-03-08 Thread Gabor Grothendieck
scan has no problem with blank lines. read.table has an argument that controls how it handles blank lines and the default setting is to ignore them. --- Date: Mon, 8 Mar 2004 20:27:48 -0600 (CST) From: Rui Song [EMAIL PROTECTED] To: [EMAIL PROTECTED] Subject: [R] About reading data

RE: [R] time-series

2004-03-08 Thread antonio rodriguez
Hi Eduardo, Probably a good idea is to do a multivariate analysis through spectral analysis. With function spectrum you can look for coherence and phase between several time series. In your case with different lengths, you can use the option: na.action=na.omit

Re: [R] Decision Trees

2004-03-08 Thread Prof Brian Ripley
Please check the licence conditions for C4.5: it is basically commercial. You are not allowed to use the published source code, and I am not allowed to show you the source code I bought, which I am not allowed to `use for commercial purposes or gain'. WEKA includes a re-implementation of the