search for: soyer

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