Re: [R] how to get the 'starting directory'

2016-03-01 Thread PIKAL Petr
Hi Not for me. I do not use Rstudio nor GIT and I use R mainly for interactive statistics. The closest I came to R programming is a set of my functions transformed to a package. So you has to wait for others more experienced in R programming. Maybe you could try to change your second file with

Re: [R] path + name of the running R file

2016-03-01 Thread PIKAL Petr
Hi What do you mean by R file? When you start R session you can get actual working directory by mydir <- getwd() and info about files in this directory by list.files() Cheers Petr > -Original Message- > From: R-help [mailto:r-help-boun...@r-project.org] On Behalf Of MAURICE > Jean

Re: [R] how to get the 'starting directory'

2016-03-01 Thread MAURICE Jean - externe
Thanks Petr, ALT + CTRL + R is a shortcut in RSTUDIO to RUN ALL As I shall work with another developer, I installed GIT (a source manager) but now R 'bombs'. I have 3 versions of my R sources (in fact a file with the 'master' and one with all the functions) : one project within Git, one project

Re: [R] how to get the 'starting directory'

2016-03-01 Thread PIKAL Petr
Hi I am not an expert in making R scripts but what you see probably means that your R session starts in users\jm\mydocuments directory which is set as your working directory. AFAIK there are basically 2 directories - one with location of R program D:\programy\R-devel\bin\i386\Rgui.exe and one

Re: [R] Help in R code

2016-03-01 Thread PIKAL Petr
Well, so again, my first answer was maybe too elaborated. 0. Send your replies also to help list. 1. Do not post in HTML. 2. Send some data by using output from dput(yourobject) 3. Prepare desired output and again show it by dput(yourresult) I still do not understand your desired output. In first

[R] path + name of the running R file

2016-03-01 Thread MAURICE Jean - externe
Hi, Is it possible to get the path and name of the running R file ? Jean Ce message et toutes les pièces jointes (ci-après le 'Message') sont établis à l'intention exclusive des destinataires et les informations qui y figurent sont strictement confidentielles. Toute utilisation de ce Message

Re: [R] help in maximum likelihood

2016-03-01 Thread ProfJCNash
It's useful to add "print.level=2" inside your call to find that there's essentially nothing wrong. Rvmmin doesn't give the msg and numDeriv gives a similar (BUT NOT EXACTLY THE SAME!) hessian estimate. It's almost always worthwhile turning on the diagnostic printing when doing optimization, even

[R] help in maximum likelihood

2016-03-01 Thread Alaa Sindi
Hi all, what is wrong with this code? I am trying to estimate the model parameters by maximizing the likelihood function and I am getting this warning Warning message: In nlm(fn, p = c(-50, 20), hessian = TRUE) : NA/Inf replaced by maximum positive value x <- c(1.6907, 1.7242, 1.7552, 1

[R] ggplot2-theme(panel.border)

2016-03-01 Thread Hoji, Akihiko
A quick question. How come data points disappear when “theme(panel.border=element_rect()) is used and how could I get those points back on the plot ? Thanks in advance. AH __ R-help@r-project.org mailing list -- To UNSUBSCRIBE and more, see https:/

Re: [R] removing factor values in the main data frame

2016-03-01 Thread peter dalgaard
> On 01 Mar 2016, at 18:33 , David L Carlson wrote: > > You need to learn how to send emails in plain text since html gets mangled on > r-help. See your message below. If I understand your question, it has to do > with what happens to factor levels when you subset your data. Subsetting a > fa

Re: [R] FW: Multivariate ARIMA

2016-03-01 Thread Thomas Lofaro
Hi Paul, Thanks so much for reaching out. I read through the link but I am not having any luck. It seems the first part of the link refers to a single variable ARIMA model with no exogenous variable. Then the link starts discussing forecasting. Anyway, I have two columns of data in my matr

Re: [R] Help in R code

2016-03-01 Thread Deepak Aggarwal
Hahaha this is not homework..i am just trying to explore r vs sas Sent from my iPhone > On 01-Mar-2016, at 9:36 PM, Bert Gunter wrote: > > ... but if this is homework (it looks like it) ask your teachers for > help, as there is a no homework policy on this list. > > Cheers, > Bert > > > Bert

Re: [R] plotting spline term when frailty term is included

2016-03-01 Thread David Winsemius
> On Mar 1, 2016, at 2:59 AM, Conor Edward Little wrote: > > Dear colleagues, > > I'd very much appreciate your help in resolving a problem that I'm having > with plotting a spline term. > > I have a Cox PH model including a smoothing spline and a frailty term as > follows: > > fit<-coxph

Re: [R] Copula Regression

2016-03-01 Thread Martin Maechler
> Janmaat, John > on Tue, 1 Mar 2016 04:53:43 + writes: > Well, seem to have solved own problem. > The article "Enjoy the joy of copulas: with a package copula" has a nice appendix that describes precisely how to set up what I am looking for. The only minor issue is tha

Re: [R] Fw: removing factor values in the main data frame

2016-03-01 Thread Jeff Newmiller
I advocate not converting to factor in the first place. Delay that until you won't be wanting specific levels to be accounted for. -- Sent from my phone. Please excuse my brevity. On March 1, 2016 9:33:31 AM PST, David L Carlson wrote: >You need to learn how to send emails in plain text since

Re: [R] Fw: removing factor values in the main data frame

2016-03-01 Thread David L Carlson
You need to learn how to send emails in plain text since html gets mangled on r-help. See your message below. If I understand your question, it has to do with what happens to factor levels when you subset your data. Subsetting a factor does not remove empty factor levels. This is documented on t

Re: [R] removing factor values in the main data frame

2016-03-01 Thread Ulrik Stervbo
Without more information it is hard to see where things go wrong. To judge from your text, I am not sure you do things correct. You are doing it along the lines of: str(mydata) submydata <- subset(mydata, ... ) submydata <- droplevels(submydata) str(submydata) right? On Tue, 1 Mar 2016 at 17:45

Re: [R] removing factor values in the main data frame

2016-03-01 Thread Ulrik Stervbo
Hi, I think droplevels(df) is what you are looking for. Best wishes, Ulrik On Tue, 1 Mar 2016 at 17:33 hoda rahmati via R-help wrote: > Hi all,I have the following main data frame:(mydata)$ TE : num 40 > 40 20 20 20 20 20 20 20 40 ...$ TR : num 49 49 28 28 28 28 28 28 28 > 49 .

Re: [R] removing factor values in the main data frame

2016-03-01 Thread hoda rahmati via R-help
I also tried droplevels but after getting str(submydata) again I see no changes in COUNTRY  On Tuesday, March 1, 2016 8:38 AM, Ulrik Stervbo wrote: Hi, I think droplevels(df) is what you are looking for. Best wishes, Ulrik On Tue, 1 Mar 2016 at 17:33 hoda rahmati via R-help wrote:

Re: [R] source() is 'bombing'

2016-03-01 Thread Ben Tupper
Hi, I'm a little fuzzy on what a bomb is in this context, but the source() function expects the input argument to be either a connection object or character string. See ... ?source Try enclosing you filename with quotes. Cheers, Ben > On Mar 1, 2016, at 11:29 AM, MAURICE Jean - externe >

[R] Fw: removing factor values in the main data frame

2016-03-01 Thread hoda rahmati via R-help
Hi all,I have the following main data frame:(mydata)        $ TE : num 40 40 20 20 20 20 20 20 20 40 ...        $ TR : num 49 49 28 28 28 28 28 28 28 49 ...         $ COUNTRY : Factor w/ 27 levels "","AU","BA","BE",..: 8 8 8 8 8 ...among the COUNTRY I just need US and AU,first I get a subs

[R] removing factor values in the main data frame

2016-03-01 Thread hoda rahmati via R-help
Hi all,I have the following main data frame:(mydata)        $ TE : num 40 40 20 20 20 20 20 20 20 40 ...        $ TR : num 49 49 28 28 28 28 28 28 28 49 ...        $ COUNTRY : Factor w/ 27 levels "","AU","BA","BE",..: 8 8 8 8 8 ...among the COUNTRY I just need US and AU,first I get a subset to c

[R] source() is 'bombing'

2016-03-01 Thread MAURICE Jean - externe
Hi, This morning, we transfered 2 R files into Git because someone else is going to work with me. A file is the 'main script' and the second contains the functions. I had problem with the encoding of this second file : French accent were replaced by question mark. Without Git, I could 'reopen wi

Re: [R] Help with SpaceCAP

2016-03-01 Thread Bert Gunter
**Always post back to the list**, which I have cc'ed, not just to me. There is obviously some problem in your data files. As I have no idea what their structure is or should be, I have no clue. Cheers, Bert Bert Gunter "The trouble with having an open mind is that people keep coming along and

Re: [R] Help in R code

2016-03-01 Thread Bert Gunter
... but if this is homework (it looks like it) ask your teachers for help, as there is a no homework policy on this list. Cheers, Bert Bert Gunter "The trouble with having an open mind is that people keep coming along and sticking things into it." -- Opus (aka Berkeley Breathed in his "Bloom Co

Re: [R] Heteroscedastic data due to responses at one factor level being too low to measure and assumed to be 0?

2016-03-01 Thread Bert Gunter
This is a list about the R (statistical) programming language, not a statistical advice list. That would be (among others) stats.stackexchange.com, to which you can try posting. However, recommendation: As you stated, you are well out of your statistical depth here, and I would strongly advise tha

Re: [R] R package for bootstrapping a mixed-design (between-within) MANOVA

2016-03-01 Thread Bert Gunter
R is a programming LANGUAGE. You can always write your own code (it is not clear exactly how you would bootstrap, as this is not an iid error structure), and there are many online courses and tutorials to help you do so. See here for some suggestions: https://www.rstudio.com/resources/training/onl

Re: [R] Help in R code

2016-03-01 Thread PIKAL Petr
Hi See in line > -Original Message- > From: R-help [mailto:r-help-boun...@r-project.org] On Behalf Of deepak > aggarwal > Sent: Tuesday, March 01, 2016 1:30 PM > To: r-help@r-project.org > Subject: [R] Help in R code > > Hi , > > Seeking your help in coding following requirement in R. > >

[R] Heteroscedastic data due to responses at one factor level being too low to measure and assumed to be 0?

2016-03-01 Thread Dr. Leigh S. Sutherland
Hello all, I have a 3-factor experimental design. The response is the load to break an adhesive joint. I have a question about Heteroscedastic data ('non-uniform variability', Google told me that's what it's called - I'm not a statistician, just an engineer) *In essence:* - At one level of a

[R] how to get the 'starting directory'

2016-03-01 Thread MAURICE Jean - externe
Hi, it's very difficult for me to understand the fullrefman.pdf document !!! I open a R script that is in the Git\reposit\gesdyn\jm directory. I 'type' alt + ctrl + R to run it. I would like to get 'Git\reposit\gesdyn\jm' in a variable ... Getwd() gives something like users\jm\mydocuments Jean

Re: [R] legend for vectorplot in rasterVis

2016-03-01 Thread Adrienne Wootten
All, Thanks anyway folks, but I'm going to call myself a bonehead and move on now that I've found it. The key.arrow argument in vectorplot will do what I need to make a scale legend. Thanks all! Adrienne On Tue, Mar 1, 2016 at 9:14 AM, Adrienne Wootten wrote: > Jim, > > Thanks! Interesting

Re: [R] legend for vectorplot in rasterVis

2016-03-01 Thread Adrienne Wootten
Jim, Thanks! Interestingly when working with lengthKey it gives me an error that plot.new hasn't been called after the plot has been created with vectorplot. Really bizarre to me. A On Tue, Mar 1, 2016 at 5:22 AM, Jim Lemon wrote: > Hi Adrienne, > I'm not sure if this will help, but lengthKe

[R] Help in R code

2016-03-01 Thread deepak aggarwal
Hi , Seeking your help in coding following requirement in R. Vector 1 has sentences for ex. vector 1="he is a nice human being","he is smart and fast in his work" vector 2 has keywords found in each sentence seperated by comma for ex. vector2 for vector 1 =nice,being vector 2 for vector 1 =smar

[R] RSNNS neural network

2016-03-01 Thread jake88
I am new to R and neural networks . So I trained and predicted an elman network like so :  require ( RSNNS )  mydata = read.csv("mydata.csv",header = TRUE)  mydata.train = mydata[1000:2000,]  mydata.test = mydata[800:999,]  fit <- elman ( mydata.train[,2:10],mydata.train[,1], size =100       lear

[R] plotting spline term when frailty term is included

2016-03-01 Thread Conor Edward Little
Dear colleagues, I'd very much appreciate your help in resolving a problem that I'm having with plotting a spline term. I have a Cox PH model including a smoothing spline and a frailty term as follows: fit<-coxph(Surv(start,end,exit) ~ x + pspline(z) + frailty(a)) When I run a model withou

Re: [R] Color of points in legend() ignored if plotting to PNG

2016-03-01 Thread Helmut Schütz
Hi Boris, Jim, Petr, and David, I stand corrected; works as desired. In order to come up with an example I oversimplified my original code which contained lwd=0 in the legend. lwd=0 indeed switches the 'borders' of points off. Case closed. Helmut -- Ing. Helmut Schuetz BEBAC - Consultancy S

Re: [R] legend for vectorplot in rasterVis

2016-03-01 Thread Jim Lemon
Hi Adrienne, I'm not sure if this will help, but lengthKey in the plotrix package will display a scale showing the relationship of vector length to whatever numeric value is being displayed. However, you do have to sort of the scaling manually. Jim On Tue, Mar 1, 2016 at 7:30 AM, Adrienne Wootte

Re: [R] divide polygon shapefile into 3 equal areas

2016-03-01 Thread Shane Carey
This is super, thanks!! However, I cannot read in my shapefile. I am using readShapeSpatial and the error I receive is: Error in getinfo.shape(fn): Error opening SHP file. The projection is Irish Transverse Mercator!! Thanks in advance On Mon, Feb 29, 2016 at 6:35 PM, Barry Rowlingson < b.rowlin

Re: [R] Version 3.2.3: package not available error with https

2016-03-01 Thread Loris Bennett
Hi Bjørn-Helge, Bjørn-Helge Mevik writes: > Loris Bennett writes: > >> It seems that R needs libcurl 7.28.0, but my platform (Scientific Linux >> 6.7) only provides version 7.19.7. > > We got "bit" by this when upgrading to 3.2.2. If you cannot upgrade > libcurl on your machine(s), you can put

Re: [R] Version 3.2.3: package not available error with https

2016-03-01 Thread Bjørn-Helge Mevik
Loris Bennett writes: > It seems that R needs libcurl 7.28.0, but my platform (Scientific Linux > 6.7) only provides version 7.19.7. We got "bit" by this when upgrading to 3.2.2. If you cannot upgrade libcurl on your machine(s), you can put local({ options(useHTTPS = FALSE) }) in the Rpro