Re: [R] Plot 2 different time series data with 2 y axes

2013-06-02 Thread Enrico Schumann
On Sun, 02 Jun 2013, Piyush Kumar writes: > Hello > > I am new to R and I have been trying to plot 2 different time series with 2 > y axes. I have been able to plot it so far with 2 y axis but the x axis is > not shifted. I mean, the x range for data1 from day1 and xrange for data2 > is from day2

Re: [R] 3.0.1 update and compiler package

2013-06-02 Thread Pascal Oettli
Hi, How did you upgraded your version of R? From source or from a Linux package? Regards, Pascal On 05/31/2013 11:33 AM, Christoph Knapp wrote: Hi, I recently updated to R 3.0.1. I'm running linux ubuntu 12.04. I realized that I have to update all the installed packages so I run > update.pack

Re: [R] 3.0.1 update and compiler package

2013-06-02 Thread Christoph Knapp
Hello, would it be possible to make a decision whether this forwarded to the list or declined. Its now 2 days. Thanks Christoph On 31/05/13 14:33, Christoph Knapp wrote: Hi, I recently updated to R 3.0.1. I'm running linux ubuntu 12.04. I realized that I have to update all the installed pac

Re: [R] X-axis in meandist {vegan}

2013-06-02 Thread Jari Oksanen
Dear Kumar Mainali, On 03/06/2013, at 00:01 AM, Kumar Mainali wrote: > I am using vegan package for estimating dissimilarity within and between > clusters, and plotting them. The following code gives me a dendrogram without > X-axis. How would I add X axis to the plot? Below "dataEnv" is a mat

[R] Conditional Forecast with MSBVAR

2013-06-02 Thread Devin Moore
Hi there! I'm trying to use the R package MSBVAR to generate a conditional forecast from a bayesian VAR. The condition that I am imposing for the forecast is observed values for all but one of the variables. I then want to forecast the path of that one variable, given observed values for all o

[R] Strange behaviour of R graphics copied to PowerPoint

2013-06-02 Thread Erling Johan Frøysa
Hello, I am using R to create graphics, especially to plot time series charts. These charts are then copied as metafiles (for best quality) to a PowerPoint presentation and then saved to PDF (via the "Save As" dialog"). Attached is two pictures. The first picture shows how my chart looks like in

Re: [R] Regularized Discriminant Analysis scores, anyone?

2013-06-02 Thread Matthew Fagan
Thank you Dr. Ligges, i very much appreciate the quick reply. i wondered if that was the case, based on the math as I (poorly) understood it. However i remain confused. page 107 from the "rrcov" package PDF makes me think I can derive LDA-style discriminant scores for a QDA: library(rrcov

[R] Problem with package 'mi' for multiple imputation

2013-06-02 Thread Sascha Sandragesan
Dear R-Community, I have just started using R (3.0.0) for my statistics and I find it the best programm for running my regressions. However, as I am still a junior user I soon came to the limit of my technical knowhow. Here's my situation: My dataset consists of 30 observations and 100 variab

[R] X-axis in meandist {vegan}

2013-06-02 Thread Kumar Mainali
I am using vegan package for estimating dissimilarity within and between clusters, and plotting them. The following code gives me a dendrogram without X-axis. How would I add X axis to the plot? Below "dataEnv" is a matrix of my environmental variables. dataFact$Source_by_Type has the levels that d

Re: [R] Scale package change comma defaults?

2013-06-02 Thread Ben Bolker
John Kane inbox.com> writes: > I was using the comma() funtion in the scales page and was > wondering how to chage a the defaults. comm(1000) gives me 1,000 > which is what I usually want but how would I change the output to > 1.000. Since the comma() function is just a wrapper for format():

[R] Scale package change comma defaults?

2013-06-02 Thread John Kane
I was using the comma() funtion in the scales page and was wondering how to chage a the defaults. comm(1000) gives me 1,000 which is what I usually want but how would I change the output to 1.000. I had thought that I could simply do comm(1000, big.mark = ".") but I am getting Error in form

[R] useR meetup group in Munich, Friday 7th: Web Service Frameworks with R

2013-06-02 Thread Markus Schmidberger
Dear all, I would like to invite Munich (Germany) area R users for our second meeting: 7th June 2013. The group is aimed to bring together practitioners from industry and academia in order to exchange knowledge and experience in solving data analysis & statistical problems by using R. More info

Re: [R] measuring distances between colours?

2013-06-02 Thread John Fox
Dear Kevin, When computer code is bug free, we'll probably all be out of business. Thank you for improving my original code. Best, John > -Original Message- > From: r-help-boun...@r-project.org [mailto:r-help-bounces@r- > project.org] On Behalf Of Kevin Wright > Sent: Sunday, June 02, 2

Re: [R] Multivariate EWMA covariance estimator?

2013-06-02 Thread Berend Hasselman
On 02-06-2013, at 19:45, Neuman Co wrote: > Again, a big thanks for your answer. > > On this webpage: > http://financetrainingcourse.com/education/2010/03/master-class-calculating-value-at-risk-var-final-steps/ > > I found, that I have to rescale by dividing the weights calculated in > Step B2

Re: [R] Multivariate EWMA covariance estimator?

2013-06-02 Thread Neuman Co
Again, a big thanks for your answer. On this webpage: http://financetrainingcourse.com/education/2010/03/master-class-calculating-value-at-risk-var-final-steps/ I found, that I have to rescale by dividing the weights calculated in Step B2 by 1-?n The "?" is the lambda, since the webpage cannot d

Re: [R] Multivariate EWMA covariance estimator?

2013-06-02 Thread Berend Hasselman
On 02-06-2013, at 19:03, Neuman Co wrote: > Thanks a lot for your answer, one more question: > I now use 100 values, so not infinity values. That means I cut some > values off, so the weights will not sum up to one. With which factor > do I have to multiply the (1-lambda)*summe2 to rescale it? S

Re: [R] Multivariate EWMA covariance estimator?

2013-06-02 Thread Neuman Co
Thanks a lot for your answer, one more question: I now use 100 values, so not infinity values. That means I cut some values off, so the weights will not sum up to one. With which factor do I have to multiply the (1-lambda)*summe2 to rescale it? So that I do not always underestimate the variance any

Re: [R] Official way to set/retrieve options in packages?

2013-06-02 Thread Bert Gunter
Have you read ?Startup which explains startup procedures, including reading of .Rprofile. I cannot speak for Brian Ripley, but I can tell you that I would prefer that packages did not mess with my .Rprofile files. OTOH, I have no objections if packages suggest that I set certain options in them

Re: [R] Official way to set/retrieve options in packages?

2013-06-02 Thread Jonathan Greenberg
What would be an example of setting, saving, and re-loading an option to a user's .Rprofile -- and would this be a no-no in a CRAN package? --j On Sat, Jun 1, 2013 at 4:57 PM, Prof Brian Ripley wrote: > On 01/06/2013 22:44, Anthony Damico wrote: > > hope this helps.. :) > > > > # define a

Re: [R] measuring distances between colours?

2013-06-02 Thread Kevin Wright
Sorry about the bug. How embarrassing. Especially because I've learned over the years to trust my gut feelings when something doesn't feel quite right, and when I was testing the function, I remember thinking "surely there a better matching named color than 'magenta'". Thanks for the fix. Kevin

Re: [R] Regularized Discriminant Analysis scores, anyone?

2013-06-02 Thread Uwe Ligges
On 02.06.2013 05:01, Matthew Fagan wrote: Hi all, I am attempting to do Regularized Discriminant Analysis (RDA) on a large dataset, and I want to extract the RDA discriminant score matrix. But the predict function in the "klaR" package, unlike the predict function for LDA in the "MASS" packa

Re: [R] Multivariate EWMA covariance estimator?

2013-06-02 Thread Berend Hasselman
On 02-06-2013, at 15:17, Neuman Co wrote: > Hi, > since I want to calculate the VaR of a portfolio consiting of 4 assets > (returns saved into "eonreturn","henkelreturn" and so on) I have to > estimate the covariance matrix. I do not want to take the rectangular > version with equal weights, but

Re: [R] Windows Vista computer will not use changes I made to R console file

2013-06-02 Thread Uwe Ligges
On 02.06.2013 13:49, Milan Bouchet-Valat wrote: Le samedi 01 juin 2013 à 15:28 -0400, Julie Royster a écrit : My newer laptop died so I installed R on an older Vista laptop. I did follow the online instructions about giving myself full control under the R properties security tab. I made change

[R] Robust GAM that covers Gaussian Distributions

2013-06-02 Thread Christos Giannoulis
Dear All, I was looking the r-archives and crantastic...for a package that has a robust approach to generalized additive models. I found two packages "robustgam" and "rgam" but their implemented functions cover only binomial and poisson distributions (pls correct me if I am wrong). I would greatl

[R] Multivariate EWMA covariance estimator?

2013-06-02 Thread Neuman Co
Hi, since I want to calculate the VaR of a portfolio consiting of 4 assets (returns saved into "eonreturn","henkelreturn" and so on) I have to estimate the covariance matrix. I do not want to take the rectangular version with equal weights, but the exponentially weighted moving average in a multiva

Re: [R] Windows Vista computer will not use changes I made to R console file

2013-06-02 Thread Milan Bouchet-Valat
Le samedi 01 juin 2013 à 15:28 -0400, Julie Royster a écrit : > My newer laptop died so I installed R on an older Vista laptop. > I did follow the online instructions about giving myself full control under > the R properties security tab. > I made changes in the R console file under the etc folder

Re: [R] metafor mixed-effect models: testing different predictor combinations for lowest AIC

2013-06-02 Thread Michael Dewey
At 22:12 01/06/2013, emhanlon wrote: Good afternoon, I am trying to use a mixed-effect model to find the best predictors of wildlife response (fledgling success, offspring growth rate, etc...) to human disturbance. I am using a Log Response Ratio as my measure of effect size along with study va

Re: [R] object not found

2013-06-02 Thread Jim Holtman
but that is not the variable that is not found. "mCaT_soma_AB" Sent from my iPad On Jun 2, 2013, at 4:12, Jannetta Steyn wrote: > mCaT_soma_AB __ R-help@r-project.org mailing list https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read the postin

Re: [R] object not found

2013-06-02 Thread Jannetta Steyn
Hi Jim It is defined and initialised: line 69 in the differential equation: dmCaT_soma_AB <- (mCaTx(v_soma) - mCaT_soma_AB)/taumCaT(v_soma) and initialised in line 83 dmCaT_soma_AB=0 However, l just noticed the mistake. It needs to be initialised without the d, i.e. mCaT_soma_AB=0 and not dmCaT_

Re: [R] Computing Median for Subset of Data

2013-06-02 Thread Berend Hasselman
On 02-06-2013, at 08:08, Matt Stati wrote: >> From my larger data set I created a subset of it by using: > > subset_1 <- subset(timeuse, IndepTrans = 1, Physical = 1) > > where my larger data set is "timeuse" and the smaller subset is "subset_1". > The subset was conditioned on "IndepTrans"

Re: [R] Computing Median for Subset of Data

2013-06-02 Thread peter dalgaard
On Jun 2, 2013, at 08:08 , Matt Stati wrote: >> From my larger data set I created a subset of it by using: > > subset_1 <- subset(timeuse, IndepTrans = 1, Physical = 1) That's not going to work. Try subset(timeuse, (IndepTrans == 1) & (Physical == 1)) (Whenever you do things like this, run s

[R] Computing Median for Subset of Data

2013-06-02 Thread Matt Stati
>From my larger data set I created a subset of it by using:  subset_1 <- subset(timeuse, IndepTrans = 1, Physical = 1) where my larger data set is "timeuse" and the smaller subset is "subset_1". The subset was conditioned on "IndepTrans" equaling "1" in the data and "Physical" equaling "1" as w