similar to: mFilter package

Dear all, firstly I would like to say I am a beginner user for R and also naturally a new member of this mail list. This means I have newer read previous mails before. I need a small help for my analysis. I want to perform a multinominal logistic regression for my data set. But as I told before I am a beginner for R and I couldn't find the code for multinominal logistic

This is not the correct group to post this question. Sorry, but I subscribe nowhere else. What is the BLAS routine to multiply a matrix by a scalar? Thanks. -- Mehmet Balcilar, PhD Assistant Professor Manas University College of Economics & Administrative Sciences Department of Economics Prospect Tinctik 56 Bishkek Kyrgyzstan Tel: +996 (312) 54 19 42 +996 (312) 54 19 43 +996

Dear R users I am currently investigating time series analysis using an irregular time series. Our study is looking at vegetation change in areas of alien vegetation growth after clearing events. The irregular time series is sourced from Landsat ETM+ data, over a six year period I have 38 scenes. For certain periods I have monthly data while for others, images are up to three months apart. So far

I think filter function in ts library has a typo. Line 34 in filter.R reads if(NROW(init) != 1 && NROW(init) != nser) I guess it should be if(NCOL(init) != 1 && NCOL(init) != nser) Please let me know, if I am wrong. Mehmet Balcilar ______________________________________________________ Get Your Private, Free Email at http://www.hotmail.com

Full_Name: Henrik Bengtsson Version: 1.5.1 OS: WinMe Submission from: (NULL) (217.210.0.243) In the Perl script $R_HOME/bin/check there is a bug under the section "Check R code for syntax errors" where the 'Rfiles <- c(...)' is build up. If there are too many files in @Rfiles the source code line generated will be too long and weird things will happen, e.g. strange

hello everyone we are developing a package that has worked up until R3.0 which we just tested. the issue is as above .... when we call a function that works in R 2.15.2 from R 3.0 we get an error Error in .Call("antsImageRead", filename, pixeltype, dimension) : "antsImageRead" not resolved from current namespace (ANTsR) this Error occurs when the .Call is made from

On Wed, 14 Aug 2002, Henrik Bengtsson wrote: > Sorry, but it was indeed the redirection of the standard output in > Cygwin/bash that cause the first problem, not R (I should stop doing > troubleshooting at 1:00 AM). So please forget about the problems reported in > R_CMD_check.out. However, it would still be nice if you still update R CMD > check to do join with "\n".

This message is in MIME format. The first part should be readable text, while the remaining parts are likely unreadable without MIME-aware tools. Send mail to mime@docserver.cac.washington.edu for more info. ------=_NextPart_000_0006_01C2432C.7ACD6BA0 Content-Type: TEXT/PLAIN; CHARSET=iso-8859-1 Content-Transfer-Encoding: 8BIT Content-ID: <Pine.GSO.4.44.0208140920312.15226@auk.stats>

Hi, here are three results of random forest (version 4.0-1). The results seem to be more or less the same which is strange because I changed the classwt. I hoped that for example classwt=c(0.45,0.1,0.45) would result in fewer cases classified as class 2. Did I understand something wrong? Christian x1rf <- randomForest(x=as.data.frame(mfilters[cvtrain,]),

Hi, the following strange error appears when I use qda: > qda1 <- qda(as.data.frame(mfilters[cvtrain,]),as.factor(traingroups)) Error: function is not a closure That's also strange: > qda1 <- qda(mfilters[cvtrain,],as.factor(traingroups)) Error in qda.default(mfilters[cvtrain, ], as.factor(traingroups)) : length of dimnames must match that of dims Some backgroud: >

heya i've been using freebsd's ipfw for quite a while and recently on a new server i've got this issue with ipfw that i can't understand ... something is wrong ... 01000 8042 1947866 allow ip from any to any via fxp0 01010 0 0 allow ip from any to any via lo0 01014 9886 4170269 divert 8668 ip from any to any in via vr0 01015 0 0 check-state 01130 14679 5695969 skipto 1800 ip from

Hello, This is my first post in this e-mail list and I hope it's enough to justify calling for help. In case it's not, sorry. I'm trying to do analysis and graphics using a factor as a criteria to split data and do the analysis/graphics for each subset of data. Right now what I'm trying to do is to fit and plot the following logistic model, according to a third variable named

Hi all, If I use a moving average, it will smooth the choppy time series, but it will lead to lagging... How do I smooth timeseries without the lagging effect? Thanks!

Hi everyone, I am trying to fit a sine function on one year of wind data. I have two questions below. Looking around on the net I managed to get the following: Sine Equation: y = a + b * sin( c + d*x ) b is the amplitude, c is the phase shift, d is something deal with periodicty of data*.* This can be linearised by sin( c+dx ) = cos(c) * sin(dx) + sin(c) * cos(dx). If one calls dx = x1 y

Dear All, How to do Hodrick-Prescott Filter in R? Thanks! [[alternative HTML version deleted]]

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

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:

Could someone, please, write me, how to compute the spar-value for the smooth.spline-routine to get the same HP-filtered time-series with a parameter lambda for a function (see mail "Hodrick-Prescott-Filter example" from ggrothendieck at yifan.net): hpf <- function(y,lambda{ eye <- diag(length(y)) d <- diff(eye,d=2) z <- solve(eye+lambda*crossprod(d),y)} ? Second, is

Can someone please show me how to smooth time series data that I have in the form of a zoo object? I have a monthly economies series and all I really need is to see a less jagged line when I plot it. If I do something like s <- smooth.spline(d.zoo$Y, spar = 0.2) plot(predict(s,index(d.zoo)), xlab = "Year") # not defined for Date objects and if I do something like

Hola a todos Estoy tratando de analizaruna serie temporal con datos diarios de temperaturas (40 años). He creado un objeto zoo (con ayuda de la lista, gracias) sobre el que encuentro la regresión lineal. He probado también a crear un objeto ts a partir del zoo. El problema que encuentro es que nose puede aplicar la función stl para hallar la componente estacional y la tendencia. Rdice que la