[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 well. I want

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

Re: [R] Computing Median for Subset of Data

2013-06-02 Thread Berend Hasselman
On 02-06-2013, at 08:08, Matt Stati mattst...@yahoo.com 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

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

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 janne...@henning.org wrote: mCaT_soma_AB __ R-help@r-project.org mailing list https://stat.ethz.ch/mailman/listinfo/r-help PLEASE

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

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 to

[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

[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 greatly

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

Re: [R] Multivariate EWMA covariance estimator?

2013-06-02 Thread Berend Hasselman
On 02-06-2013, at 15:17, Neuman Co neumanc...@gmail.com 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

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 package,

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] 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 rip...@stats.ox.ac.ukwrote: On 01/06/2013 22:44, Anthony Damico wrote: hope this helps.. :)

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

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

Re: [R] Multivariate EWMA covariance estimator?

2013-06-02 Thread Berend Hasselman
On 02-06-2013, at 19:03, Neuman Co neumanc...@gmail.com 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

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

Re: [R] Multivariate EWMA covariance estimator?

2013-06-02 Thread Berend Hasselman
On 02-06-2013, at 19:45, Neuman Co neumanc...@gmail.com 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

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, 2013

[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

[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

Re: [R] Scale package change comma defaults?

2013-06-02 Thread Ben Bolker
John Kane jrkrideau at 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

[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

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

[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

[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

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 matrix

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

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