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.
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
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
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
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
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
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
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
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
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
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
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
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
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
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
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
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
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
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,
?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.
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!
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
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
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
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)
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
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
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)
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
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
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
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
-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
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
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
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
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
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
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
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
__
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
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
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 -
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
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
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
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,
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
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
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
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
51 matches
Mail list logo