search for: jonesw

Hi, Does anyone know if R has the functionality to calculate a simple moving average. I cant seem to find it in the help menu. thanks, Wayne Dr Wayne R. Jones Statistician / Research Analyst KSS Group plc St James''s Buildings 79 Oxford Street Manchester M1 6SS Tel: +44(0)161 609 4084 Mob: +44(0)7810 523 713 KSS Ltd A division of Knowledge Support Systems Group plc Seventh

Dear all, Is there a test in R for the randomness of a sequence of observations (e.g. to test the random number generator)? Specifically I am looking for autocorrelations which are not necessarily linear in nature, which the acf function does not seem to be flexible enough to detect as it tests for linear autocorrelation. Thanks in advance, Paul.

Hi there fellow R-users, I'm stuck on this seemingly trivial problem. All I want to coerce a character string into a command. For example: x<-rnorm(20) y<-rnorm(20) str<-"lm(y~x)" I want to evaluate the "str" command. I have tried eval(as.expression(str)) But it doesn't seem to work. I am aware of the call command, but for reasons I won't go

Dear all, Could anyone please tell me how to sort a matrix or a data frame against a column/row/component? Many thanks. -MY

Hi there, Does anyone know if R has a function for testing whether a time series is stationary?? Thanks in advance, Wayne Dr Wayne R. Jones Statistician / Research Analyst KSS Group plc St James's Buildings 79 Oxford Street Manchester M1 6SS Tel: +44(0)161 609 4084 Mob: +44(0)7810 523 713 KSS Ltd A division of Knowledge Support Systems Group plc Seventh Floor St James's

Hi there, Does anyone know how to explicitly refer to the p-value of thet test that the chosen number of factors is significant in a factor analysis. It's not in the list of values for the factanal command output yet it is printed out with the results. Thanks in advance. Wayne Dr Wayne R. Jones Statistician / Research Analyst KSS Group plc St James's Buildings 79 Oxford Street

Hi there fellow R-users, Can anyone tell me how to build a prediction interval for a gaussian log-link model for the reponse variable?? I can find the standard error of the predictions but I cant seem to find the prediction interval. Is there a way I can calculate the prediction interval from the standard errors?? Here's the example: logX<-rnorm(100)

Hi there fellow R-users, I am trying to follow an example of modelling a serial correlation structure in the textbook "Mixed Effects Model in S and Splus". However, I am getting some very odd results. Here is what I am trying to run: library(nlme) data(Ovary) fm1<-lme(follicles~sin(2*pi*Time)+cos(2*pi*Time),data=Ovary,random=pdDiag(~s in(2*pi*Time))) ### The example is fine up

Hello there fellow R-users, I am using the RODBC functionality to query a database. I am trying to read in a columns of strings which have a character field lengths greater than 255. The data.frame that I recieve back from the RODBC query only contains the first 255 characters (the rest having been truncated). Any help on how to solve this problem would be greatly appreciated. Reagrds Wayne

www.kssg.com a pricing consulatancy based in Manchester (UK) are seeking an experienced statistical analyst. Candidates must be eligible to work in the EU. If you are interested please send an updated CV to jonesw@kssg.com . Company Overview KSS is a leading provider of pricing and revenue management systems for the Retail and Petroleum sectors. Our business applications help our clients get more out of their existing pricing strategies, as well as improve the efficiency and productivity of their pricing...

Hi there fellow R-users, Can anyone tell me if there exits an R package that deals with serial correlation in the residuals of an lm model. Perhaps, using the Cochrane Orcutt or Praise Wilson methods? Thanks, Wayne Dr Wayne R. Jones Senior Statistician / Research Analyst KSS Limited St James's Buildings 79 Oxford Street Manchester M1 6SS Tel: +44(0)161 609 4084 Mob: +44(0)7810 523 713

Hi there R-Helpers, Does anyone know if it is possible to sort a dataframe? I.e. Sort alphabetically column 1 ( which has some reocurring elements) then sort alphabetically column2 but keeping the order of column 1 constant; much the same way that the sort function works in Excel. Regards, Wayne Dr Wayne R. Jones Statistician / Research Analyst KSS Group plc St James's Buildings 79

Hello there fellow R-users, I have received some data which comes in the following format: example1<-"200301" The first 4 digits correspond to the year and the remaining 2 digits correspond to the week of the year. I have tried to convert this to a date by using strptime as follows: strptime(example1,format="%Y%U") where U (looking up strptime) is the week of the

Does anyone know how to coerce a numeric to a string?? THanks Wayne KSS Ltd A division of Knowledge Support Systems Group plc Seventh Floor St James''s Buildings 79 Oxford Street Manchester M1 6SS England Company Registration Number 2800886 (Limited) 3449594 (plc) Tel: +44 (0) 161 228 0040 Fax: +44 (0) 161 236 6305 mailto:kssg@kssg.com http://www.kssg.com The information in

Hello All, Does anyone out there know a way to decompose time series objects with missing values. A simple "na.omit" will not work since it does not preserve the time differences between succesive observations. Thanks in advance, Wayne Dr Wayne R. Jones Statistician / Research Analyst KSS Group plc St James''s Buildings 79 Oxford Street Manchester M1 6SS Tel: +44(0)161

HI there, Does anyone how I can build my own library of R functions? Regards, Wayne Dr Wayne R. Jones Statistician / Research Analyst KSS Group plc St James's Buildings 79 Oxford Street Manchester M1 6SS Tel: +44(0)161 609 4084 Mob: +44(0)7810 523 713 KSS Ltd A division of Knowledge Support Systems Group plc Seventh Floor St James's Buildings 79 Oxford Street Manchester M1

Hi R-Users, I have very recently switched from using R on windows 2000 to windows XP Professional. With XP I am having problems using the "odbcConnect("")" command. I can use the command perfectly well until I close the database channel with the command "odbcClose(channel)". When I try and re-establish connection following the close command no dialog box appears

> Hi there fellow R-Users, > > I am trying to use the "package.skeleton" to create my own package with > R.1.7.1 on Windows XP Professional. > I have followed the package.skeleton example and have downloaded the > necessary files found at http://www.stats.ox.ac.uk/pub/Rtools/tools.zip. > > and perl5, available via

Hello there fellow R users, Im trying to fit a gls model to data which has serial correlation in the errors e(t)=p*e(t-1). However I dont seem to be having much luck in erradicating the autocorrelation in the residuals. I have created the following example. library(nlme) x<-rnorm(100) y<-3+2*x y<-y+arima.sim(100,model=list(ar=(0.6)))+rnorm(100,0,0.2) #Create a data set with first

Hi there fellow R users, Has anyone got an R example of applying an Ideal demand system, possibly using the library systemfit?? Thanks Wayne Dr Wayne R. Jones Senior Statistician / Research Analyst KSS Limited St James's Buildings 79 Oxford Street Manchester M1 6SS Tel: +44(0)161 609 4084 Mob: +44(0)7810 523 713 KSS Ltd Seventh Floor St James's Buildings 79 Oxford Street