similar to: VAR question

Hi all, I am running a bivariate regression with the following: p1=c(184,155,676,67,922,22,76,24,39) p2=c(1845,1483,2287,367,1693,488,435,1782,745) I1=c(1530,1505,2505,204,2285,269,1271,298,2023) I2=c(8238,6247,6150,2748,4361,5549,2657,3533,5415) R1=I1-p1 R2=I2-p2 x1=cbind(p1,R1) y1=cbind(p2,R2) fit1=lm(y1~-1+x1) summary(fit1) Response 2: Coefficients: Estimate Std. Error t value

Hi all, I'd like to fit a multivariate regression with the variance of the error term porportional to the predictors, like the WLS in the univariate case. y_1~x_1+x_2 y_2~x_1+x_2 var(y_1)=x_1*sigma_1^2 var(y_2)=x_2*sigma_2^2 cov(y_1,y_2)=sqrt(x_1*x_2)*sigma_12^2 How can I specify this in R? Is there a corresponding function to the univariate specification lm(y~x,weights=x)??

Hi all, I am going to sample two variables from two exponential distributions, but I want to specify a covariance structure between these two variables. Is there any way to do it in R? Or is there a "Multivariate Exponential" thing corresponding to the multivariate normal? Thanks in advance. Sincerely, Yanwei Zhang Department of Actuarial Research and Modeling Munich Re America Tel:

Hi all, Assume I have a random vector with four variables, i.e. A=(a,b,c,d). I am able to get the covariance matrix of vector A, but how can I get the covariance matrix of vector B=(a,c,b,d) by manipulating the corresponding covariance matrix of A? Thanks. Sincerely, Yanwei Zhang Department of Actuarial Research and Modeling Munich Re America Tel: 609-275-2176 Email:

Hi all, I have n independent variables A_1, A_2, A_3,......,A_n, and each with known variances var(A_1), var(A_2),..., but unknown mean. How can I get the approximation of the variance of the product of the variables using numerical computation, i.e. var(A_1*A_2*A_3*.....*A_n)? Thanks. Sincerely, Yanwei Zhang Department of Actuarial Research and Modeling Munich Re America Tel: 609-275-2176

Hi all, I have a 4 by 4 matrix, and I want to switch row 2 and row 3 first, then switch column 2 and column 3. Is there an easy way to do it? The following is a tedious way to get what I want. But I wonder if there is a way to simplify this. > a=matrix(rnorm(16),4,4) > a [,1] [,2] [,3] [,4] [1,] 0.33833811 -0.9422273 -0.06181611 -1.8346134 [2,] -0.68167996

Hi all, I met with a problem about the weighted least square regression. 1. I simulated a Normal vector (sim1) with mean 425906 and standard deviation 40000. 2. I simulated a second Normal vector with conditional mean b1*sim1, where b1 is just a number I specified, and variance proportional to sim1. Precisely, the standard deviation is sqrt(sim1)*50. 3. Then I run a WLS regression without the

Hi all, Is there an easy way to do cumulative matrix multipliation in R? What's the syntex? Thanks. Sincerely, Yanwei Zhang Department of Actuarial Research and Modeling Munich Re America Tel: 609-275-2176 Email: yzhang@munichreamerica.com<mailto:yzhang@munichreamerica.com> [[alternative HTML version deleted]]

Dear list, I have the following matrix. How can I make the following plot? 1. The x-axis has index 1:7, and the first column is plotted against index 1, second against 2, and so on. 2. I want the points from the left upper conner including the antidiagonal to be plotted with col=2, and the rest with col=3 [,1] [,2] [,3] [,4] [,5] [,6] [,7] [1,] 0.589 0.857 0.923 0.944 0.954 0.963

Hello, I ran AIC for some competing models I created. I get df and an AIC score from the AIC procedure. Can I use the models with the lowest AIC scores from this procedure to choose my 'best' models? If not, what else do I need to do (and know) and how can I do it in R to chose the 'best' models? Thank you kindly, Michael [[alternative HTML version deleted]]

Hello all, I looking at package dse or vars or mAr I know how to simulate a VAR(p) process, my problem is that most of those processes are unstable (not weakly stationary). Do anybody know how to generate a random VAR (or VARMA even better) process that is weakly stationary? Thanks -- View this message in context: http://r.789695.n4.nabble.com/simulating-stable-VAR-process-tp4261177p4261177.html

Hallo, as I see there is a new version for R available. Can anyone tell me how I can update my version 2.5.0 under Windows? The last times I just uninstalled the old version and installed the new one. Afterwards I had to install also all needed packages again. All in all it cost me half a day until my system works fine again. Is there a quicker option? If yes please tell me the commands. Thanks,

Dear R-experts, It is a quite stupid question but please help me. I am very confuced. I am able to import normal txt ant mat-files to R but unable to import .xls-file I do not understand the online help. Can please anyone send me the corresponding command lines? The .xls-file is attached. In my file we use commas for the decimal format (example: 0,712), changes might be needed. Thanks, Corinna

>From the R Studio downloads, look below the installers. This is off topic however. If there is no zipped, no exe, no installation required of R, then I thank you very much for your help and trolling. (BTW, I think my question was pretty clear, concise and specific, I appreciate that some of you tried to solve a problem related to what I have, but I have already reviewed all options, and what

I have a Mac OS X system. To deal with a long monthly electricity demand time-series I use the procedures TRAMO/SEATS with the MS-windows only Demetra programme and X12 under R resorting to the awkward - as far as the output is concerned - x12 R package running the relating Fortran code. I wonder if someone out there has attempted to translate TRAMO/SEATS and X12 into R native language? Ciao

When I do ARMA(2,2) using one lag of LCPIH data This is eview result > > *Dependent Variable: DLCPIH > **Method: Least Squares > **Date: 08/12/11 Time: 12:44 > **Sample (adjusted): 1970Q2 2010Q2 > **Included observations: 161 after adjustments > **Convergence achieved after 14 iterations > **MA Backcast: 1969Q4 1970Q1 > ** > **Variable Coefficient Std.

Hi R users Having both a faster CPU and more memory will boost computing power. I was wondering if only adding more memory (1GB -> 2GB) will significantly reduce R computation time? Taka, _________________________________________________________________ Get FREE company branded e-mail accounts and business Web site from Microsoft Office Live

Dear R People: Is it possible to run R on a netbook/mini, please? Thanks, Erin -- Erin Hodgess Associate Professor Department of Computer and Mathematical Sciences University of Houston - Downtown mailto: erinm.hodgess at gmail.com

I recently installed R 2.6.2 and am getting errors on startup that relate to svIDE being loaded by Tinn-R. Loading required package: tcltk Loading Tcl/Tk interface ... done Warning messages: 1: '\A' is an unrecognized escape in a character string 2: unrecognized escape removed from ";for Options\AutoIndent: 0=Off, 1=follow language scoping and 2=copy from previous line\n" 3: In

I'm using stableFit from the package fBasics to estimate the parameters of a truncated normal distribution (I'm interested in the parameters of the underlying normal distribution). It is correct to generalize this truncated normal distribution as a stable distribution ? Thanks David -- View this message in context: