similar to: Installing packages.

Has anyone written any code for the Hodrick-Prescott filter? I have a some uncompiled FORTRAN code from Ed Prescott but I'd like to save myself some programming time if possible. Thanks for your help. Nick Davis Crown Financial Policy Asset and Liability Management Branch The New Zealand Treasury Direct: +64-4-471-5924 Fax: +64-4-499-0143 Email:

I've had a play with this and, due to my own short-comings, remain none the wiser. In particular, I'm not sure what value of 'spar' is consistent with the magic lambda=1/1600 for quarterly data. I initially interpreted spar as lambda and tried setting spar=1/1600. This results in almost no smoothing while spar=1600 causes an error. The smooth.spline function seems to want

I am using the ts package to fit ARIMA models, and the tsdiag() function to plot diagnostics. In doing so I'm generating an awful lot of diagnostic plots of different models and different data sets all within the same R session. So my question is, is there an option in tsdiag() similar to <main="Title"> that I can use? This would be quite helpful when I print out the plots,

Hi folks, we are experiencing a problem when AutoCAD 2000+ drawing files are saved on a Samba (2.0.8) file share. Regardless of what the file is 'saved as', the name is always converted to all capitals. e.g. Drawing1 becomes DRAWING1 Just wondering if anyone else has experienced this problem. Seems strange that it is only AutoCAD files, but I think it is Samba related, not AutoCAD,

We are experiencing long delays when loading MX files stored on our samba fileserver. We are running version 2.2.7a of Samba on Solaris 8 servers. Using MX version 2.5 and higher on windows 2000 and XP workstations When the files are stored locally they load in a few seconds, whilst taking more than a minute to run from the Samba server. Our current workaround is to store the files locally,

The Stern Environmental Review, a British Government Report into tackeling global climate change has been published online. Leader page with presentation and speaking notes :- http://www.hm-treasury.gov.uk/independent_reviews/stern_review_economics_climate_change/sternreview_index.cfm The Stern Review Report :-

HI All, I am trying to analyse the results of a heap of samples, to determine how many samples are required to meet a defined confidence interval. I checked out the R manual for bootstrapping and re-sampling methods, but did not find anything. Are there any contributed procedures that deal with this? thanks, Matt Redding Soil Scientist Intensive Livestock Environmental Management-Research

Dear All, I like to know how to simulate survival data using median (or mean) survival time. Any help will be greatly appreciated. Best wishes, Kere Kerenaftali Klein PhD| Biostatistician | Queensland Clinical Trials & Biostatistics Centre The University of Queensland | School of Population Health | Building 33, Level 1| Princess Alexandra Hospital |Ipswich Road | Woolloongabba QLD 4102 |

I have experienced a couple of (minor) problems with Samba Thought I would report these and ask if they are known problems. 1. Problem with "!" in username map file 2. Samba tools e.g. smbstatus, testparm not accounting for "config file" parameter 1. There appears to be a problem with a minor part of username map file interpretation. When I put an exclamation mark

Could someone who is familiar with the book "Samba-3 HOWTO and reference guide and who is also familiar with printer setup please give me some advice... I am summarising the procedure for setting up a printer for automatic driver download which is described in chapter 17. There are a couple of processes mentioned that look very similar and I am wondering if both steps are required for

-----BEGIN PGP SIGNED MESSAGE----- Hash: SHA1 I have the following data: > y [1] 0.000 0.004 0.008 0.016 0.024 0.032 0.044 0.064 0.072 0.088 0.108 0.140 [13] 0.156 0.180 0.208 0.236 0.264 0.296 0.320 0.360 0.408 0.444 0.472 0.524 [25] 0.576 > x [1] 100 200 300 400 500 600 700 800 900 1000 1100 1200 1300 1400 1500 [16] 1600 1700 1800 1900 2000 2100 2200 2300 2400 2500 I'd

Hi All R guru's, A couple of simple questions which I cannot find answers to in "an introduction to R" or in the reference manuals. I am running R on Win95, and have plotted up quite a nice box and whisker plot: char<-read.table("n:/reddinm/grifpr~1/survey/char.txt") attach(char) boxplot(Al,Fe,OC,pH, range=0,xaxis=FALSE, ylab="Proportion of Unirigated

Hi we're running 3.0.1 on Solaris 9 ( with NIS/flat files as the NS ) as a member server of the AD domain ( via kinit and then net join ). there's a couple of things we've noticed and I'm not sure if they're just the way it works or configuration problems: (1) we assign the gid an uid mappings with idmap in smb.conf and I thought that winbindd would not assign uid/gids if

We have a problem with file security settings changing when a file is edited by someone who is not the owner. When the edited file is saved, ownership is transferred to the person who edited the file. Also, access control settings for the file are reset to a default setting. This is a real problem, since in some cases it means that the original owner of the file no longer has permission to edit

I have 2 zoo objects - 1) Interest rate spread between 10-YR-US-Treasury and 2-YR-US-Treasury (object name = sprd) 2) S&P 500 index (object name = spy) > str(spy) ?zoo? series from 1976-06-01 to 2011-03-31 Data: num [1:8791] 99.8 100.2 100.1 99.2 98.6 ... Index: Class 'Date' num [1:8791] 2343 2344 2345 2346 2349 ... > str(sprd) ?zoo? series from 1976-06-01 to 2011-03-31

Hi, I'm having trouble getting an answer to the following problem on -questions - I hope someone here has done something similar and can help. I'd like to compile support for FreeBSD OPIE into sshd. Presently I have to use PAM to achieve one-time password support. On a 4.x system I have in /etc/ssh/sshd_config ChallengeResponseAuthentication yes and in /etc/pam.conf sshd auth

Dear R users, I'm having trouble using aggregate() and would greatly appreciate your advice. I am using R 1.5.1 on Windows 2000. I want to call my function in the following way extract.data(x=dat[, "Age", "Year"]) where extract.data() uses aggregate() to count the number of cases for each combination of "Age" and "Year". I've defined

-----BEGIN PGP SIGNED MESSAGE----- Hash: SHA1 I'm generating a symmetric correlation matrix using a data matrix as input: mat <- cor(data.mat) My question is: Is there a more memory efficient way to store this data? For instance, since: all(mat == t(mat)) every value is duplicated, and I should be able to almost half the memory usage for large matrices. Any thoughts/comments? Cheers,

Dear list, I have a data frame of survey respondents, a little like this: set.seed(20081215) n <- 100 dat <- data.frame(id=1:100, addr1=sample(LETTERS, n, replace=TRUE), addr2=sample(LETTERS, n, replace=TRUE), addr3=sample(LETTERS, n, replace=TRUE)) head(dat) id addr1 addr2 addr3 1 1 R H Q 2 2 H C K 3 3

Hi All, I've been trawling through the documentation and listserv archives on this topic -- but as yet have not found a solution. I'm sure this is pretty simple with R, but I cannot work out how without resorting to ugly nested loops. As far as I can tell, grep, match, and %in% are not the correct tools. Question: given these vectors -- patrn <- c(1,2,3,4) exmpl <-