Displaying 17 results from an estimated 17 matches for "soyer".

Did you mean:
boyer

2013 Oct 16

7

Is there something wrong with R version 3.0.2 (2013-09-25) -- "Frisbee Sailing"?

Hi,
pnorm(-1.53,0,1) under version 3.0.2 gives 0.05155075. I am pretty sure it
should be 0.063. Is there something wrong with this version of R?
I am using:
R version 3.0.2 (2013-09-25) -- "Frisbee Sailing"
Copyright (C) 2013 The R Foundation for Statistical Computing
Platform: i686-pc-linux-gnu (32-bit)
--
Tom
[[alternative HTML version deleted]]

2008 May 09

1

Which gls models to use?

Hi,
I need to correct for ar(1) behavior of my residuals of my model. I noticed
that there are multiple gls models in R. I am wondering if anyone
has experience in choosing between gls models. For example, how
should one decide whether to use lm.gls in MASS, or gls in nlme for
correcting ar(1)? Does anyone have a preference? Any advice is appreciated!
Thanks,
--
Tom
[[alternative HTML

2007 Apr 27

1

acf and pacf plot

Hi,
I noticed that whenever I ran acf or pacf, the plot generated by R always
includes two horizontal blue doted lines. Furthermore, these two lines are
not documented in the acf documentation. I don't know what they are for, but
it seems that they are important. Could someone tell me what they are and
how are they calculated?
Thanks,
--
Tom
[[alternative HTML version deleted]]

2008 Feb 02

1

ARCH LM test for univariant time series

Hi,
Does anyone know if R has a Lagrange multiplier (LM) test for ARCH
effects for univariant time series?
Thanks!
--
Tom
[[alternative HTML version deleted]]

2007 Dec 27

1

questions about plot.zoo

Hi,
I have been having very good results using plot.zoo to chart time series
data. But I have three questions about plot.zoo and I am wondering if anyone
knows the answers.
(1) when I tried to use semi-log scale, via log="y", R issued a warning,
although it looked like plot.zoo plotted in semi-log scale anyway:
Warning message:
In plot.xy(xy.coords(x, y), type = type, ...) :

2006 Nov 01

4

Problem with data type recognition and conversion

Hi,
I have a CSV file with two columns; the first column is date, second column
is numbers. I used read.csv() to load the file into the variable temp.
Somehow, R could not recognize my numbers as double. Instead, it thinks
these numbers are integer even though they all have decimal points (isn't
that strange?). The problem I ran into is that if I tried to convert the
numbers to double using

2008 Jan 06

1

aggregate.ts help

Hi,
I have a ts object with a frequency of 4, i.e., quarterly data, and I would
like to calculate the mean for each quarter. So for example:
> ts.data=ts(1:20,start=c(1984,2),frequency=4)
> ts.data
Qtr1 Qtr2 Qtr3 Qtr4
1984 1 2 3
1985 4 5 6 7
1986 8 9 10 11
1987 12 13 14 15
1988 16 17 18 19
1989 20
If I do this manually, the mean

2008 Jan 26

1

Is there a safe mode?

Hi,
My R just froze. I can't get it to do anything. It gives "Error: band value"
message to everything I type. Does anyone know if R has a safe mode that I
could check for errors and perform diagnostics? I am using R 2.6.1 on
Windows XP.
> ls()
Error: bad value
> search()
Error: bad value
> ?ls
Error: bad value
> q()
Error: bad value
--
Tom
[[alternative HTML version

2008 Jan 24

2

plot help

Hi,
Suppose I already have two plots on the same screen, and I want to draw
lines in each of them. Is that possible in R? It seems that once you have
two plots on the screen, you can only draw lines in the the last plot, never
the 1st. Here is what I mean:
#some data
y1=rnorm(1:3)
y2=rnorm(1:3)
#draw two plots on the same screen
par(mfrow=c(2,1),oma = c(6, 0, 5, 0))
par(mar=c(0, 5.1, 0, 5.1))

2008 Jan 19

1

show overstruck in plot legend

Hi,
I have a plot with type="o", or overstruck. Now I am trying to add the
legend, but I couldn't figure out how to show the overstruck type in the
legend. It seems that the legend only allows one to set lty. Does anyone
know how to show overstruck in the legend?
Thanks!
--
Tom
[[alternative HTML version deleted]]

2007 Apr 28

1

pacf

Hi,
I wanted to understand exactly how acf and pacf works, so I tried to
calculate ac and pac manually. For ac, I used the standard acf formula:
acf(k) = sum(X(t)-Xbar)(X(t-k)-Xbar))/sum(X(t)-Xbar)^2. But for pac, I could
not figure out how to calculate it by hand. I understand that in both R and
EVIEWS, it is done using the Durbin-Levinson algorithm by the computer.
However, I don't

2008 Apr 08

2

diagonally fill a rectangle with color gradient

Hi,
Is it possible to diagonally fill a rectangle with a color gradient? I
noticed that the gradient.rect of plotrix could fill a rect either up and
down or from side to side. I am looking for something similar but fills
diagonally instead, e.g., from the upper left corner to the bottom right.
Does anyone know how to do it in R?
Thanks,
--
Tom
[[alternative HTML version deleted]]

2007 Nov 30

1

what to do if residuals produced by lm() have long tails?

Hi,
I am using lm() for regression analysis of my data set. My regression
results look pretty good, i.e., the coefficient is significant and the p
value is much less than 0.05. But when I checked the residuals, both using
qqnorm() and hist(), the distribution does not look normal. It looks like
the residuals have long tails. I assume that lm() uses OLS, and since one of
the assumptions of OLS is

2008 May 09

1

RMSE and lm

Hi,
Does anyone know if the RMSE is one of the values provided by the lm model,
or do we have to calculate it by hand from the residuals?
Thanks,
--
Tom
[[alternative HTML version deleted]]

2008 Jul 27

1

help with durbin.watson

Hi,
I have two time series, y and x. Diff(y) and Diff(x) both show no
autocorrelation. But durbin.watson(lm(Diff(y)~lag(Diff(x),k=-4)) gives a DW
value of zero. How come the residule is autocorrelated while Diff(y) and
Diff(x) are not? Does anyone know if in my case a DW of zero indicates
serial correlation, or is it telling me that the DW statistics is not the
appropriate statistics to use here?

2006 Nov 19

1

Dovecot (unofficial) patches

...beta
and 1.0 allows remote attackers to list files and directories under the
mbox parent directory and obtain mailbox names via ".." sequences in the
(1) LIST or (2) DELETE IMAP command.
I'm not a spec file wizard, so I change as less as possible. Unfortunately
the maintainer (Jerome Soyer) has never responded to my emails.
My server runs Mandriva Official 2007.0 and Postfix 2.3.3 and dovecot rc14
(will move to rc15 soon).
TNX
EgbertJan (NL)
-------------- next part --------------
A non-text attachment was scrubbed...
Name: smime.p7s
Type: application/x-pkcs7-signature
Size: 38...

2008 Jan 14

0

question about order selection of ar

Hi,
I have an AR(1) series, so I thought that the order of the series should be
1. A simple lm fit with one period lag predicts the series pretty well. But
when I tried ar, I got different orders: ar.mle selected order 6,
ar.burgselected order 14, and
ar.yw selected order 6. So I am wondering maybe I miss understood what
"order" means in R? Does anyone know why would ar select order 6