similar to: AMRAtoMA

Hello R Users! I have a question about the output of ARMAtoMA when used to calculate the variance of a model.? I have a mixed model of the form ARMA(1,1).? The actual model takes the form: X(t) = 0.75X(t-12) + a(t) - 0.4a(t-1) Given that gamma(0) takes the form [(1 + theta^2 - 2*theta*phi)/(1-phi^2)]*sigma(a), I would expect a process variance of 4.02*sigma(a) when I substitute 0.75 for phi and

Hello all- I would appreciate any guidance that can be provided. I am new to R and am using it exclusively in a statistics program I am undertaking that mainly references Minitab. My focus is on data modeling and further more multivariate data analysis as much of my work in involves chemical measurements from custom sensors using all sorts of transduction methods. I am looking for a

Hello R users wiser than I - I am trying to produce a boxplot with quantiles defined by the type 6 algorithm used mainly by Minitab and SPSS for a comparison study. I found how to compare the results in tabular form using the quantile(x, type = 6) function, but am stuck on how to show the results in a summary boxplot form. Would someone be able the help with this situation? Many thanks, Matt

Dear friends - I am trying to understand fda and working with the 2009 book from Springer. I run the scripts directly from a new installed fda library and include sessionInfo() below. This is from the script fdarm-ch09.R - I seem to have got not fda_2.2.8 but 2.2.7 - is that the problem? Where do I get .8 as I just reinstalled? Best wishes Troels Ring, Nephrology Aalborg, Denmark >

Hi R users, Someone knows how to replace Infinite value by zero. I have a vector with some Inf value and I want to substitute these values by zero to get the mean of the components of the vector. Any idea? Many thanks, Marlene. [[alternative HTML version deleted]]

How can I use a function with some additional input parameters in constrOptim? For example, something like fr <- function(x,a) { ## Rosenbrock Banana function x1 <- x[1] x2 <- x[2] a * (x2 - x1 * x1)^2 + (1 - x1)^2 } where the optimum is to be found w.r.t. x. Calling optim(c(-1.2,1), fr, NULL, a=100) works as expected, but I fail to provide the a=100 in the constrained case:

Dear R-user Suppose I have the following data y=c(2,1,5,8,11,3,1,7,50,21,33,7,60) x=data.frame(y) for(i in 4:nrow(x)) x[i,] =var(x[i-3:i-1,]) I'm trying to get a new variable with the variance of the 3 previous values (just an example) and with NA in the three first positions. I know that my for() is wrong but I'm not able to find my error. Any idea? Thanks, Marlene.

Hi R-users I'm writing a code to run a fuction but I found an error that I can't fix. I reproduced the error with a simple example. The correct answer is k but I can't fill my s matrix. What I'm doing wrong? s<-matrix(data=NA,nrow=1,ncol=5 ) s for(i in 1:5) { k=sqrt(i) s[,i]<-k[i] print(k) } s Thanks in advance, Marlene. [[alternative

Full_Name: Madeleine Yeh Version: 1.9.1 OS: AIX 5.2 Submission from: (NULL) (151.121.225.1) After compiling R-1.9.1 on AIX 5.2 using the IBM cc compiler, I ran the checks. One of them failed. Here is the output from running the check solo. root@svweb:/fsapps/test/build/R/1.9.1/R-1.9.1/tests/Examples: ># ../../bin/R --vanilla < stats-Ex.R R : Copyright 2004, The R

Hi, I need to repor the index of a min value of each row in a matrix, but I don't know how to do that when I have more than one min value. Here is my example > dat <- matrix(c(5.4,4.8,5.6,4.8,NA,4.4,4.6,3.4,NA,NA,4,2.4,NA,NA,NA,2),byrow=TRUE,ncol=4) > dat [,1] [,2] [,3] [,4] [1,] 5.4 4.8 5.6 4.8 [2,] NA 4.4 4.6 3.4 [3,] NA NA 4.0 2.4 [4,] NA NA NA 2.0

Hello, I am relatively new to R and have a rudimentary question, I think. How does one truncate the number of digits displayed after the decimal when viewing the results of analyses? My apologies if this question has been answered previously, I was not able to find references very easily. Thank you in advance, Matt [[alternative HTML version deleted]]

Hi, I am translating some Gauss code to R. Gauss has an interesting way of handling constraints. Observe the following code snipplet: e1 = x[.,23] .eq 0; @ remove obs with Regular Hours = 0 @ e2 = x[.,12] .gt 1; @ remove obs with non-regular work status @ e3 = x[.,4] .lt 15; @ remove obs with agricultural and mining industry code (< 15)@ esum = e1 + e2 + e3; e = esum .gt 0; @

I have been trying some different sample rate and bitrate combinations to get a feel for how CELT behaves. I then encoded it with a couple of different frame sizes and sample rates. There were some small differences, but nothing horrible. So, I decided to run with 24KHz sample rate with 16 bit samples. This seemed like a reasonable tradeoff against the quality degradation that my folks

Hi, i try to install, but getting no success. What this mean cannot find -lz , any necessary lib left in my SuSE Installation? Many thanks, Christian linux:/usr/lib/R/bin # R CMD INSTALL RMySQL_0.5-2.tar.gz * Installing *source* package 'RMySQL' ... creating cache ./config.cache checking how to run the C preprocessor... cc -E checking for mysql_init in -lmysqlclient... yes checking

Dear List, I just tried installing R 1.1.1 on SuSE 7.0 and encountered some problems. Using the precompiled binaries, I had a failed dependency: R was looking for "libblas.so.2" and apparently didn't find it. I then tried compiling from sources, which failed with the message: "/usr/bin/ginstall: zero.so No such file or directory" (or something to that effect). I'm

Queridos miembros de la lista, Me gustaría anunciar que el paquete ''SCperf " ya está disponible en CRAN. El paquete implementa diferentes modelos de inventario, el efecto látigo y otras variables de desempeño de la cadena de suministro. Esta es la primera versión del paquete y también es mi primer paquete para R. Sugerencias, informes de errores y otros comentarios son

Dear R users, I'm running a SVR in package e1071 but I did not able to calculate the parameters w and b of the regression. I don't know how to do that and if it is possible to do it with this package. Someone have some idea. Any help would be much appreciated. Marlene [[alternative HTML version deleted]]

I have a data frame (df) with colums x, y and z. e.g. df <- data.frame(x = sample(4), y = sample(4), z = sample(4)) I can extract column z by: df$z or df[3] I can also extract columns x,y by: df[1:2] or by df[-3]. Is it possible to extract x,y columns in a "symbolic" fashion i.e. by equivalent of df[-z] (which is illegal) ??? Or alternativeley, is there an equivalent of

Hi all, I''m trying to get text on to the screen in response to pressing a button. This code doesn''t work, and I don''t know why. My puts command says that the window is shown, but it''s not displayed. Also, it seems to me that there should be an easier way to get text on a screen in response to pressing a button. Any hints? Thanks! Marlene require

Dear Professor, I am preparing for a Ph.D in semiparametric regression at Cairo university in Egypt. Referring to R package KernGPLM, I obtained R version 2.4.1 but I did not find package KernGPLM. Please, help me how can I obtain this package. Thanks in advance. Name: Magda Haggag E-mail: magdahaggag@yahoo.com Address: 27, Notrdam Desion st., Gleem, Alexandria, Egypt.