similar to: lag() and lm()

Full_Name: Ahmad Abu Hammour Version: rw0651 OS: windows 95 Submission from: (NULL) (63.23.128.44) Although I know that "ts package" is preliminary, I wanted to compare the results from R and SPSS. I ran ARIMA(2,1,2) in both softwares. I got NaN in standard errors of coefficients from R and real figures from SPSS. I changed "delta" in R to match that used by SPSS, I received

Hi, P-values in Durbin-Watson test obtained through the use of functions available in packages "lmtest" and "car" are different. The difference is quite significant. function "dwtest" in "lmtest" is much faster than "burbinwatson" in "car". Actually, you can take a nap while the latter trying to calculated Durbin-Watson test. My question

Hi, I used RODBC to import data from an Excel spreadsheet "*.xls", but some columns were returned as zeros. When I looked at these columns in Excel, I found that thery are results of formula calculations and not entry. My question is: Is there any parameter or command I need to use in order to overcome this problem? Thank you Ahmad Abu Hammour

Hi, Is there a function in any of R-packages that can produce and plot the impulse response function for any model.. Thank you Ahmad Abu Hammour -.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.- r-help mailing list -- Read http://www.ci.tuwien.ac.at/~hornik/R/R-FAQ.html Send "info", "help", or "[un]subscribe" (in the "body",

Hi, I noticed that paste and expression in the following manner does not produce the desired result. vnames=c("pi","y","delta") vx=vn=3 m=16 # mydata any matrix that takes for the following form mydata=array(somedata,c(m,vn,vx)) # producing plots for (j in 1:vx){ for (i in 1:vn) { plot(mydata, type="l", main=expression(paste("Effect of",

Hi, I have the following hypothetical optimization problem: -det(t(x-A%*%x1)%*%(x-A%*%x1)) where A,x,x1 are matrices. A coefficients and x and x1 are variable matrices or vectors. I tried to apply optim and nlm functions but I kept receive the following massage: Error in A%*%x1 : non-conformable arguments. The massage appears even the -det() can be calculated and the dimensions are checked. here

Hi, I have a matrix and time series "xmat". I have no problem executing any matrix functions except t(xmat) which gives the following error message. My question if "xmat" is a matrix why I can not execute this matrix function?? converting the time series and matrix into a data frame solves the problem >dim(xmat) [1] 98 7 > t(xmat) Error in "tsp<-"(*tmp*,

Hi, I need help regarding estimating a linear model where restrictions are imposed on the coefficients. An example is as follows: Y_{t+2}=a1Y_{t+1} + a2 Y_t + b x_t + e_t restriction a1+ a2 =1 Is there a function or a package that can estimate the coefficient of a model like this? I want to estimate the coefficients rather than test them. Thank you for your help Ahmad Abu Hammour --------------

Dear all, could someone please clarify me if this works, so as to implement lags and differences for example in y and in a independent x1 of a regression? model<-lm((diff(y), -i)~x1+lag(x1,-1), data=anydata) Thanks, a lot, Alexandre. ---------- In??cio da mensagem original ----------- De: r-help-bounces at stat.math.ethz.ch Para: r-help at stat.math.ethz.ch Cc:

Buenos días: Tengo un determinado código en R que en una máquina funciona correctamente, pero en otra devuelve un error del tipo "Error en as.formula(mymodel): Objeto ''mymodel'' no encontrado''" En mi código: mymodel <- "y ~ x1[, 1] + x2[, 1]" res <- lm (as.formula(mymodel)) ¿Qué puede estar motivando ese error?. La máquina en la que falla

hi, my linear model is y=c+a*x1+b*x2 i tried to found a, c, b by the use of: mymodel<-lm(y~1+x1+x2) where y, x1, x2 are 3 vectors with the same length the result is a=NA.so i want to know where is the problem. -- View this message in context: http://www.nabble.com/lm%28%29-function-tp17546079p17546079.html Sent from the R help mailing list archive at Nabble.com.

Dear list, In my previous email, the model I'd like to estimate is rea=a*st+b*mod+error, where a+b=1 and a,b>0. My apologies for the misunderstanding. Thanks for all your help, Jorge On Wed, Jul 23, 2008 at 3:35 PM, Jorge Ivan Velez <jorgeivanvelez@gmail.com> wrote: > > Dear list, > > I have a data set which looks like myDF (see below) and I'd like to >

Hi everyone, Can someone please help me in these questions?: 1)if I use crossvalidation with svm, do I have to use this equation to calculate RMSE?: mymodel <- svm(myformula,data=mydata,cross=10) sqrt(mean(mymodel$MSE)) But if I don’t use crossvalidation, I have to use the following to calculate RMSE: mymodel <- svm(myformula,data=mydata) mytest

I have dose response data analyzed with the package 'drc'. 'summary(mymodel)' prints my kinetic parameters. I want that text in an ASCII text file. I want to get exactly what I would get if I copied and pasted from the terminal window. I've read the documentation on data export to text files here:

This should be very easy, but alas, I'm very new to R. My end goal is to calculate p-values from arima(). Let's say I just ran this: > MyModel <- arima(y[1:58], order=c(1,0,0), xreg=MyData[1:58,7:14], > method="ML") > MyModel And I see: arima(x = y[1:58], order = c(1, 0, 0), xreg = MyData[1:58, 7:14], method = "ML") Coefficients: ar1

Hello again, I guess this is more of a Ruby question, but here goes anyways... I want to intercept all calls to non existent class methods and then call an existing class method with the non existent name. In code... class MyModel < ActiveRecord::Base def MyModel.method_missing(symbol, *args) MyModel.func(symbol.id2name) end def MyModel.func(func_called) puts "You

?sink On May 18, 2018 9:47:25 AM PDT, Ed Siefker <ebs15242 at gmail.com> wrote: >I have dose response data analyzed with the package 'drc'. >'summary(mymodel)' prints my kinetic parameters. I want >that text in an ASCII text file. I want to get exactly what I >would get if I copied and pasted from the terminal window. > >I've read the documentation on

I''m getting confused. If I have a view (.rhtml), I can use a field helper like: <%= text_field :mymodel, :attribute %> This uses the value from @mymodel.attribute Now, if I am rendering a template normally (i.e. not partial), I would expect @mymodel to be a model object defined as an instance variable in the controller. However, if I am rendering a partial template, I

Hi, I would like to extract x,y from mjca function in ca package to be able to draw other graph types. For example for simple mca function from MASS I can do the following: ------------- anydata.mca <- mca(anydata, nf = 2, abbrev = TRUE) x <- anydata.mca$rs[,1] #that is what I can 'extract' y <- anydata.mca$rs[,2] plot(anydata.mca) #standard plot library(lattice)

Hello! I am using package "VAR". I've fitted my model: mymodel<-VAR(mydata,myp,type="const") I can extract the Log Liklihood for THE WHOLE MODEL: logLik(mymodel) How could I calculate (other than manually) the corresponding Akaike Information Criterion (AIC)? I tried AIC - but it does not take mymodel: AIC(mymodel) # numeric(0) Thank you! -- Dimitri Liakhovitski