Re: [R] create new

[Bert Gunter]

?ifelse

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 County" comic strip )


On Fri, Mar 24, 2017 at 8:56 PM, Val  wrote:
> Hi all,
>
>
> I have several variables in a group and one group  contains three
> variables. Sample of data ( Year, x1, x3 and x2)
>
> mydat <- read.table(header=TRUE, text=' Year  x1  x3  x2
> Year1  10  120
> Year2   0  150
> Year3   0   020
> Year4  25   0   12
> Year5  15  25   12
> Year6   0  16   14
> Year7   0  100')
>
> I want create another variable( x4) based on the following condition.
>
> if x1  > 0  then x4 = x1; regardless of  x2 and x3 values.
> if x1  = 0  and x2  > 0 then x4 = x2;
> if x1  = 0 and  x2  = 0 then x4 = x3
>
> The desired output looks like as follows
> Yearx1  x3  x2   x4
> Year1  10  120  10
> Year20  150  15
> Year300   20  20
> Year4  250   12  25
> Year5  15  25   12  15
> Year60  16   14  14
> Year70  10 0  10
>
> Thank you in advance
>
> __
> R-help@r-project.org mailing list -- To UNSUBSCRIBE and more, see
> https://stat.ethz.ch/mailman/listinfo/r-help
> PLEASE do read the posting guide http://www.R-project.org/posting-guide.html
> and provide commented, minimal, self-contained, reproducible code.

__
R-help@r-project.org mailing list -- To UNSUBSCRIBE and more, see
https://stat.ethz.ch/mailman/listinfo/r-help
PLEASE do read the posting guide http://www.R-project.org/posting-guide.html
and provide commented, minimal, self-contained, reproducible code.

[R] create new

[Val]

Hi all,


I have several variables in a group and one group  contains three
variables. Sample of data ( Year, x1, x3 and x2)

mydat <- read.table(header=TRUE, text=' Year  x1  x3  x2
Year1  10  120
Year2   0  150
Year3   0   020
Year4  25   0   12
Year5  15  25   12
Year6   0  16   14
Year7   0  100')

I want create another variable( x4) based on the following condition.

if x1  > 0  then x4 = x1; regardless of  x2 and x3 values.
if x1  = 0  and x2  > 0 then x4 = x2;
if x1  = 0 and  x2  = 0 then x4 = x3

The desired output looks like as follows
Yearx1  x3  x2   x4
Year1  10  120  10
Year20  150  15
Year300   20  20
Year4  250   12  25
Year5  15  25   12  15
Year60  16   14  14
Year70  10 0  10

Thank you in advance

__
R-help@r-project.org mailing list -- To UNSUBSCRIBE and more, see
https://stat.ethz.ch/mailman/listinfo/r-help
PLEASE do read the posting guide http://www.R-project.org/posting-guide.html
and provide commented, minimal, self-contained, reproducible code.

Re: [R] change the value of vector based on the string of its name

[David Winsemius]

> On Mar 24, 2017, at 10:12 AM, Qill Bel  wrote:
> 
> Dear R-users,
> 
> Imagine I have 3 vectors:
> *proc1<-0*
> *proc2<-0*
> *cpd<-c("proc1","proc2")*

The `cpd` vector has only character values in it, and has no relation to either 
the `proc1` or `proc2` vectors. If you want it refer to  them (or contain their 
values),  then omit the quotes. The values will not retain their names in this 
instance. The R names `proc`` and `proc2` will be looked up and their values 
places in an unnamed numeric vector. You could however make it so with:

cpd<-c("proc1"=proc1,"proc2"=proc2)

#The quotes are not needed and so this would be equivalent"

cpd<-c( proc1=proc1, proc2=proc2)

There is a `get`-function, but teaching you how to use it to solve the problem 
is likely to be a disservice since it would allow you to continue using R as a 
macro-processor.

> 

> How can I change the value of proc1 to 1, based on vector cpd only?
> I tried with *as.factor(cpd[1])<-1*, but it produces an error.

Doesn't make much sense to me to use assignment to a factor in this instance, 
but that's probably because we come from different programming backgrounds. 
Neither your code nor mine would have produced a factor vector. You should say 
in natural language what it is that you want to happen.

(Perhaps) Try this:

cpd[1] <- 1

Or:

cpd['proc1'] <- 1


> 
> Any idea how could I achieve that?
> 
>   [[alternative HTML version deleted]]

And learn to post in plain text. Reading the Posting Guide will also be a good 
idea. As would be studying a bit longer the "Introduction to R" that is shipped 
with every installation.


> 
> __
> R-help@r-project.org mailing list -- To UNSUBSCRIBE and more, see
> https://stat.ethz.ch/mailman/listinfo/r-help
> PLEASE do read the posting guide http://www.R-project.org/posting-guide.html
> and provide commented, minimal, self-contained, reproducible code.

David Winsemius
Alameda, CA, USA

__
R-help@r-project.org mailing list -- To UNSUBSCRIBE and more, see
https://stat.ethz.ch/mailman/listinfo/r-help
PLEASE do read the posting guide http://www.R-project.org/posting-guide.html
and provide commented, minimal, self-contained, reproducible code.

[R] Unzip multiple layer of compressed (zip/tar) file

[Mohammad Tanvir Ahamed via R-help]

Hi, 

In a folder , there are multiple compressed (zip+tar) files. Each of compressed 
file can have multiple layer of  compression (compressed file inside the 
compressed file and so on).

Can anyone suggest, is there any way to decompress all the files in a single 
directory ?   
 
Tanvir Ahamed 
Göteborg, Sweden  |  mashra...@yahoo.com 

__
R-help@r-project.org mailing list -- To UNSUBSCRIBE and more, see
https://stat.ethz.ch/mailman/listinfo/r-help
PLEASE do read the posting guide http://www.R-project.org/posting-guide.html
and provide commented, minimal, self-contained, reproducible code.

[R] change the value of vector based on the string of its name

[Qill Bel]

Dear R-users,

Imagine I have 3 vectors:
*proc1<-0*
*proc2<-0*
*cpd<-c("proc1","proc2")*

How can I change the value of proc1 to 1, based on vector cpd only?
I tried with *as.factor(cpd[1])<-1*, but it produces an error.

Any idea how could I achieve that?

[[alternative HTML version deleted]]

__
R-help@r-project.org mailing list -- To UNSUBSCRIBE and more, see
https://stat.ethz.ch/mailman/listinfo/r-help
PLEASE do read the posting guide http://www.R-project.org/posting-guide.html
and provide commented, minimal, self-contained, reproducible code.

Re: [R] Default value of Numerals argument in read.table function R 3.3.2 on Mac

[Bert Gunter]

... and I should have added:

It is a specific instance of how ?match.arg works .

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 County" comic strip )


On Fri, Mar 24, 2017 at 1:17 PM, Bert Gunter  wrote:
> Follow the "type.convert" link in the numerals section of ?read.table.
> That's what it's there for.
>
> -- 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 County" comic strip )
>
>
> On Fri, Mar 24, 2017 at 11:55 AM, Ramnik Bansal  
> wrote:
>> Hi,
>>
>> The help file for function read.table mentions the default value for the
>> argument numerals  as
>> c("allow.loss", "warn.loss", "no.loss")
>>
>> How are the three values used as default values ?
>>
>> Thanks
>> Ramnik
>>
>> [[alternative HTML version deleted]]
>>
>> __
>> R-help@r-project.org mailing list -- To UNSUBSCRIBE and more, see
>> https://stat.ethz.ch/mailman/listinfo/r-help
>> PLEASE do read the posting guide http://www.R-project.org/posting-guide.html
>> and provide commented, minimal, self-contained, reproducible code.

__
R-help@r-project.org mailing list -- To UNSUBSCRIBE and more, see
https://stat.ethz.ch/mailman/listinfo/r-help
PLEASE do read the posting guide http://www.R-project.org/posting-guide.html
and provide commented, minimal, self-contained, reproducible code.

Re: [R] Default value of Numerals argument in read.table function R 3.3.2 on Mac

[Bert Gunter]

Follow the "type.convert" link in the numerals section of ?read.table.
That's what it's there for.

-- 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 County" comic strip )


On Fri, Mar 24, 2017 at 11:55 AM, Ramnik Bansal  wrote:
> Hi,
>
> The help file for function read.table mentions the default value for the
> argument numerals  as
> c("allow.loss", "warn.loss", "no.loss")
>
> How are the three values used as default values ?
>
> Thanks
> Ramnik
>
> [[alternative HTML version deleted]]
>
> __
> R-help@r-project.org mailing list -- To UNSUBSCRIBE and more, see
> https://stat.ethz.ch/mailman/listinfo/r-help
> PLEASE do read the posting guide http://www.R-project.org/posting-guide.html
> and provide commented, minimal, self-contained, reproducible code.

__
R-help@r-project.org mailing list -- To UNSUBSCRIBE and more, see
https://stat.ethz.ch/mailman/listinfo/r-help
PLEASE do read the posting guide http://www.R-project.org/posting-guide.html
and provide commented, minimal, self-contained, reproducible code.

Re: [R] Default value of Numerals argument in read.table function R 3.3.2 on Mac

[Sarah Goslee]

Those are the available options, not the defaults. If you follow the
chain to ?type.convert as the help for read.table suggests, you will
learn that "allow.loss" is the default value. The help for
type.convert also explains what each option means.

Sarah

On Fri, Mar 24, 2017 at 2:55 PM, Ramnik Bansal  wrote:
> Hi,
>
> The help file for function read.table mentions the default value for the
> argument numerals  as
> c("allow.loss", "warn.loss", "no.loss")
>
> How are the three values used as default values ?
>
> Thanks
> Ramnik
>

-- 
Sarah Goslee
http://www.functionaldiversity.org

__
R-help@r-project.org mailing list -- To UNSUBSCRIBE and more, see
https://stat.ethz.ch/mailman/listinfo/r-help
PLEASE do read the posting guide http://www.R-project.org/posting-guide.html
and provide commented, minimal, self-contained, reproducible code.

[R] Default value of Numerals argument in read.table function R 3.3.2 on Mac

[Ramnik Bansal]

Hi,

The help file for function read.table mentions the default value for the
argument numerals  as
c("allow.loss", "warn.loss", "no.loss")

How are the three values used as default values ?

Thanks
Ramnik

[[alternative HTML version deleted]]

__
R-help@r-project.org mailing list -- To UNSUBSCRIBE and more, see
https://stat.ethz.ch/mailman/listinfo/r-help
PLEASE do read the posting guide http://www.R-project.org/posting-guide.html
and provide commented, minimal, self-contained, reproducible code.

Re: [R] Setting up a .Rprofile file

[Henrik Bengtsson]

R for Windows is a bit peculiar where it locates your .Rprofile file,
or rather what it consider to be your home directory.  If you look
from within R, the file you do want to create / edit is:

> f <- normalizePath("~/.Rprofile", mustWork = FALSE)
> f
[1] "C:\\Users\\joe\\Documents\\.Rprofile"

For instance, to change what your R prompt looks you can add this from
within R as:

> cat('options(prompt = "> R ")\n', file = f, append = TRUE)

Hope this helps

Henrik

On Fri, Mar 24, 2017 at 2:36 AM, Bruce Ratner PhD  wrote:
> Henrico:
> Thanks for quick reply.
> However, one last question:
> If I want to change working directory, and put setwd() in the Rprofile file, 
> logically R will not know where the be work directory is, correct?
>
> So, should I install R in my preferred working directory?
>
> Thanks again, in advance.
> Bruce
>
>
> __
> Bruce Ratner PhD
> The Significant Statistician™
> (516) 791-3544
> Statistical Predictive Analytics -- www.DMSTAT1.com
> Machine-Learning Data Mining -- www.GenIQ.net
>
>
>
>> On Mar 24, 2017, at 3:48 AM, Enrico Schumann  wrote:
>>
>> On Thu, 23 Mar 2017, Bruce Ratner PhD writes:
>>
>>> Hi R'ers:
>>> I would like to setting up a .Rprofile file with
>>> setwd("C:/R_WorkDir")
>>> set.seed(12345)
>>> options (prompt "> R ")
>>>
>>> ---
>>> Can you help providing the code or instructive link,
>>> I've find many links, but I can't figure it out?
>>>
>>> Thanks.
>>> Bruce
>>>
>>
>> Quoting from ?Startup:
>>
>> ,
>> | [...] unless ‘--no-init-file’ was given, R searches
>> | for a user profile, a file of R code.  The path of
>> | this file can be specified by the ‘R_PROFILE_USER’
>> | environment variable (and tilde expansion will be
>> | performed).  If this is unset, a file called
>> | ‘.Rprofile’ is searched for in the current directory
>> | or in the user's home directory (in that order).  The
>> | user profile file is sourced into the workspace.
>> `
>>
>> --
>> Enrico Schumann
>> Lucerne, Switzerland
>> http://enricoschumann.net
>>
>
> __
> R-help@r-project.org mailing list -- To UNSUBSCRIBE and more, see
> https://stat.ethz.ch/mailman/listinfo/r-help
> PLEASE do read the posting guide http://www.R-project.org/posting-guide.html
> and provide commented, minimal, self-contained, reproducible code.

__
R-help@r-project.org mailing list -- To UNSUBSCRIBE and more, see
https://stat.ethz.ch/mailman/listinfo/r-help
PLEASE do read the posting guide http://www.R-project.org/posting-guide.html
and provide commented, minimal, self-contained, reproducible code.

Re: [R] Setting up a .Rprofile file

[Bert Gunter]

Inline.

Cheers,

Bert

On Fri, Mar 24, 2017 at 2:36 AM, Bruce Ratner PhD  wrote:
> Henrico:
> Thanks for quick reply.
> However, one last question:
> If I want to change working directory, and put setwd() in the Rprofile file, 
> logically R will not know where the be work directory is, correct?

No. See ?setwd . It requires a dir argument and will therefore give an
error if one is not supplied. Nothing to do with "knowing" what the
current working directory is.
>
> So, should I install R in my preferred working directory?

No (imho, anyway). Your working directory may change over time,
depending on how you organize things of course.

Please re-read the advice you have been given as well as ?startup.

Cheers,
Bert


>
> Thanks again, in advance.
> Bruce
>
>
> __
> Bruce Ratner PhD
> The Significant Statistician™
> (516) 791-3544
> Statistical Predictive Analytics -- www.DMSTAT1.com
> Machine-Learning Data Mining -- www.GenIQ.net
>
>
>
>> On Mar 24, 2017, at 3:48 AM, Enrico Schumann  wrote:
>>
>> On Thu, 23 Mar 2017, Bruce Ratner PhD writes:
>>
>>> Hi R'ers:
>>> I would like to setting up a .Rprofile file with
>>> setwd("C:/R_WorkDir")
>>> set.seed(12345)
>>> options (prompt "> R ")
>>>
>>> ---
>>> Can you help providing the code or instructive link,
>>> I've find many links, but I can't figure it out?
>>>
>>> Thanks.
>>> Bruce
>>>
>>
>> Quoting from ?Startup:
>>
>> ,
>> | [...] unless ‘--no-init-file’ was given, R searches
>> | for a user profile, a file of R code.  The path of
>> | this file can be specified by the ‘R_PROFILE_USER’
>> | environment variable (and tilde expansion will be
>> | performed).  If this is unset, a file called
>> | ‘.Rprofile’ is searched for in the current directory
>> | or in the user's home directory (in that order).  The
>> | user profile file is sourced into the workspace.
>> `
>>
>> --
>> Enrico Schumann
>> Lucerne, Switzerland
>> http://enricoschumann.net
>>
>
> __
> R-help@r-project.org mailing list -- To UNSUBSCRIBE and more, see
> https://stat.ethz.ch/mailman/listinfo/r-help
> PLEASE do read the posting guide http://www.R-project.org/posting-guide.html
> and provide commented, minimal, self-contained, reproducible code.

__
R-help@r-project.org mailing list -- To UNSUBSCRIBE and more, see
https://stat.ethz.ch/mailman/listinfo/r-help
PLEASE do read the posting guide http://www.R-project.org/posting-guide.html
and provide commented, minimal, self-contained, reproducible code.

Re: [R] Using betareg package to fit beta mixture with given initial parameters

[Michael Dayan]

Dear Achim,

Thank you for your help, this is exactly what I needed. Your help and
didactic example is very much appreciated.

Best wishes,

Michael

On Wed, Mar 22, 2017 at 11:34 AM, Achim Zeileis 
wrote:

> On Wed, 22 Mar 2017, Michael Dayan wrote:
>
> The method of setting the initial values given lambda, alpha1, etc. should
>> not depend on the exact values of lambda, alpha1, etc. in my situation,
>> i.e. it does not depend on my data.
>>
>
> Presently, flexmix() that betamix() is built on cannot take the parameters
> directly for initialization. However, it is possible to pass a matrix with
> initial 'cluster' probabilities. This can be easily generated using dbeta().
>
> For a concrete example consider the data generated in this discussion on
> SO:
>
> http://stats.stackexchange.com/questions/114959/mixture-of-b
> eta-distributions-full-example
>
> Using that data with random starting values requires 42 iterations until
> convergence:
>
> set.seed(0)
> m1 <- betamix(y ~ 1 | 1, data = d, k = 2)
> m1
>
> ## Call:
> ## betamix(formula = y ~ 1 | 1, data = d, k = 2)
> ## ## Cluster sizes:
> ##   1   2
> ##  50 100
> ## ## convergence after 42 iterations
>
> Instead we could initialize with the posterior probabilities obtained at
> the observed data (d$y), the true alpha/beta parameters (10; 30 vs. 30; 10)
> and the true cluster proportions (2/3 vs. 1/3):
>
> p <- cbind(2/3 * dbeta(d$y, 10, 30), 1/3 * dbeta(d$y, 30, 10))
> p <- p/rowSums(p)
>
> This converges after only 2 iterations:
>
> set.seed(0)
> m2 <- betamix(y ~ 1 | 1, data = d, k = 2, cluster = p)
> m2
>
> ## Call:
> ## betamix(formula = y ~ 1 | 1, data = d, k = 2, cluster = p)
> ## ## Cluster sizes:
> ##   1   2
> ## 100  50
> ## ## convergence after 2 iterations
>
> Up to label switching and small numerical differences, the parameter
> estimates agree. (Of course, these are on the mu/phi scale and not
> alpha/beta as explained in the SO post linked above.)
>
> coef(m1)
> ##(Intercept) (phi)_(Intercept)
> ## Comp.11.196286  3.867808
> ## Comp.2   -1.096487  3.898976
> coef(m2)
> ##(Intercept) (phi)_(Intercept)
> ## Comp.1   -1.096487  3.898976
> ## Comp.21.196286  3.867808
>
>
>
> On Mar 22, 2017 04:30, "David Winsemius"  wrote:
>>
>>
>> On Mar 21, 2017, at 5:04 AM, Michael Dayan 
>>>
>> wrote:
>>
>>>
>>> Hi,
>>>
>>> I would like to fit a mixture of two beta distributions with parameters
>>> (alpha1, beta1) for the first component, (alpha2, beta2) for the second
>>> component, and lambda for the mixing parameter. I also would like to set
>>> a
>>> maximum of 200 iterations and a tolerance of 1e-08.
>>>
>>> My question is: how to use the betareg package to run the fit with
>>> initial
>>> values for the parameters alpha1, beta1, alpha2, beta2 and lambda? I saw
>>>
>> in
>>
>>> the documentation that I would need to use the 'start' option of the
>>> betareg function, with start described as "an optional vector with
>>>
>> starting
>>
>>> values for all parameters (including phi)". However I could not find how
>>>
>> to
>>
>>> define this list given my alpha1, beta1, alpha2, beta2 and lambda.
>>>
>>> The current code I have is:
>>> mydata$y <- 
>>> bmix <- betamix(y ~ 1 | 1, data = mydata, k = 2, fstol = 1e-08, maxit =
>>> 200, fsmaxit = 200)
>>>
>>>
>>> And I suspect I would need to do something along the lines:
>>>
>>> initial.vals <- c(?, ?, ?, ?, ?)
>>> bmix <- betamix(y ~ 1 | 1, data = mydata, k = 2, fstol = 1e-08, maxit =
>>> 200, fsmaxit = 200, control=betareg.control(start=initial.vals)))
>>>
>>> But I do not know what to use for initial.vals.
>>>
>>
>> If there were sensitivity to data, then wouldn't  that depend on your
>> (unprovided) data?
>>
>>
>>
>>> Best wishes,
>>>
>>> Michael
>>>
>>>   [[alternative HTML version deleted]]
>>>
>>> __
>>> R-help@r-project.org mailing list -- To UNSUBSCRIBE and more, see
>>> https://stat.ethz.ch/mailman/listinfo/r-help
>>> PLEASE do read the posting guide http://www.R-project.org/
>>>
>> posting-guide.html
>>
>>> and provide commented, minimal, self-contained, reproducible code.
>>>
>>
>> David Winsemius
>> Alameda, CA, USA
>>
>> [[alternative HTML version deleted]]
>>
>> __
>> R-help@r-project.org mailing list -- To UNSUBSCRIBE and more, see
>> https://stat.ethz.ch/mailman/listinfo/r-help
>> PLEASE do read the posting guide http://www.R-project.org/posti
>> ng-guide.html
>> and provide commented, minimal, self-contained, reproducible code.
>>
>>

[[alternative HTML version deleted]]

__
R-help@r-project.org mailing list -- To UNSUBSCRIBE and more, see
https://stat.ethz.ch/mailman/listinfo/r-help
PLEASE do read the posting guide http://www.R-project.org/posting-guide.html
and provide commented, minimal, self-contained, 

Re: [R] Setting up a .Rprofile file

[Bruce Ratner PhD]

Henrico:
Thanks for quick reply.
However, one last question:
If I want to change working directory, and put setwd() in the Rprofile file, 
logically R will not know where the be work directory is, correct?

So, should I install R in my preferred working directory?

Thanks again, in advance. 
Bruce


__
Bruce Ratner PhD
The Significant Statistician™
(516) 791-3544
Statistical Predictive Analytics -- www.DMSTAT1.com
Machine-Learning Data Mining -- www.GenIQ.net



> On Mar 24, 2017, at 3:48 AM, Enrico Schumann  wrote:
> 
> On Thu, 23 Mar 2017, Bruce Ratner PhD writes:
> 
>> Hi R'ers:
>> I would like to setting up a .Rprofile file with
>> setwd("C:/R_WorkDir")
>> set.seed(12345)
>> options (prompt "> R ")
>> 
>> ---
>> Can you help providing the code or instructive link,
>> I've find many links, but I can't figure it out?
>> 
>> Thanks.
>> Bruce 
>> 
> 
> Quoting from ?Startup:
> 
> ,
> | [...] unless ‘--no-init-file’ was given, R searches
> | for a user profile, a file of R code.  The path of
> | this file can be specified by the ‘R_PROFILE_USER’
> | environment variable (and tilde expansion will be
> | performed).  If this is unset, a file called
> | ‘.Rprofile’ is searched for in the current directory
> | or in the user's home directory (in that order).  The
> | user profile file is sourced into the workspace.
> `
> 
> -- 
> Enrico Schumann
> Lucerne, Switzerland
> http://enricoschumann.net
> 

__
R-help@r-project.org mailing list -- To UNSUBSCRIBE and more, see
https://stat.ethz.ch/mailman/listinfo/r-help
PLEASE do read the posting guide http://www.R-project.org/posting-guide.html
and provide commented, minimal, self-contained, reproducible code.

Re: [R] How to change parameter values as a function to time with the package "deSolve"

[Marine Regis]

Thanks very much David for your answer. Sorry ! here is the code without the 
error:


library(deSolve)
param <- c(a = 0.1, b = 1)
yini <- c(alpha0 = 6, beta0 = 2)

mod <- function(times, yini, param) {

  with(as.list(c(yini, param)), {

gamma0 <- ifelse(times %in% seq(0,10,1), 5, 0)

## print(gamma0)

dalpha0 <- - a*alpha0 + gamma0
dbeta0 <- a*alpha0 - b*beta0
return(list(c(dalpha0, dbeta0)))

  })}

times <- seq(from = 0, to = 10, by = 1/24)
out <- ode(func = mod, times = times, y = yini, parms = param)
plot(out, lwd = 2, xlab = "day")

Thanks in advance for your help!
Marine




De : David Winsemius 
Envoyé : vendredi 24 mars 2017 06:01
À : Marine Regis
Cc : r-help@r-project.org
Objet : Re: [R] How to change parameter values as a function to time with the 
package "deSolve"


> On Mar 23, 2017, at 2:59 PM, Marine Regis  wrote:
>
> Hello,
>
> I am trying to solve an ODE in R using deSolve. With the following code, I 
> expected the parameter �gamma0� takes the values 5 at time step 0, 1, 2, 3, 
> 4, 5, 6, 7, 8, 9 and 10, and 0 otherwise. However, the print(gamma0) shows 
> that �gamma0� stays at 0.
>
> Here is my ODE:
>
> param <- c(a = 0.1, b = 1)
> yini <- c(alpha0 = 0, beta0 = 0)
>
> mod <- function(times, yini, param) {
>
>  with(as.list(c(yini, parameters)), {
>
>  gamma0 <- ifelse(times %in% seq(0,10,1), 5, 0)
>
>  # print(gamma0)
>
>  dalpha0 <- - a*alpha0 + gamma0
>  dbeta0 <- a*alpha0 - b*beta0
>  return(list(c(dalpha0, dbeta0)))
>
>  })}
>
> times <- seq(from = 0, to = 10, by = 1/24)
> out <- ode(func = mod, times = times, y = yini, parms = param)
> plot(out, lwd = 2, xlab = "day")
>
>
> What am I doing wrong? Thanks in advance for your help!

When I execute that code I get:

Error in as.list(c(yini, parameters)) : object 'parameters' not found

(The mistake seems fairly obvious when you look are the parameter list for your 
`mod` function.)

Additionally: It's not generally good practice to use `with` inside functions, 
but that's not the problem here.

--
David.

> Marine
>
>
>[[alternative HTML version deleted]]
>
> __
> R-help@r-project.org mailing list -- To UNSUBSCRIBE and more, see
> https://stat.ethz.ch/mailman/listinfo/r-help
R-help Info Page - Homepage - SfS – Seminar for 
Statistics
stat.ethz.ch
The main R mailing list, for announcements about the development of R and the 
availability of new code, questions and answers about problems and solutions 
using R ...



> PLEASE do read the posting guide http://www.R-project.org/posting-guide.html
> and provide commented, minimal, self-contained, reproducible code.

David Winsemius
Alameda, CA, USA


[[alternative HTML version deleted]]

__
R-help@r-project.org mailing list -- To UNSUBSCRIBE and more, see
https://stat.ethz.ch/mailman/listinfo/r-help
PLEASE do read the posting guide http://www.R-project.org/posting-guide.html
and provide commented, minimal, self-contained, reproducible code.

Re: [R] Setting up a .Rprofile file

[Jeff Newmiller]

OP is on some variant of MSWindows, which doesn't use a Bash shell, so the 
syntax of that example is a bit foreign, Enrico. 

I would say that the concept of using .Rprofile to change directories is 
ill-advised, since having ALL of your R work in one place is not sustainable. 
However you can make one or more clickable icons (shortcuts) on your desktop or 
start menu with the Startup Directory property set to the directory you want to 
work in. You can modify a shortcut made by the installer, or  hold the Alt key 
down while dragging the RGui.exe icon to some convenient place, and use the 
right-click popup menu to change Properties.

I used to create a .RData file in the directory I wanted to work in and double 
click that to get R started in that directory, but the side effects of loading 
objects into the R global environment can really confuse things and create 
puzzling behaviour if that file gets changed. These days I use RStudio and 
double click on a .Rproj (project) file to get RStudio going with the directory 
set to where that file is without mucking with the global environment. 
-- 
Sent from my phone. Please excuse my brevity.

On March 24, 2017 3:03:34 AM PDT, Enrico Schumann  
wrote:
>On Fri, 24 Mar 2017, Bruce Ratner PhD writes:
>
>> Henrico:
>> Thanks for quick reply.
>> However, one last question:
>> If I want to change working directory, and put setwd() in the
>Rprofile
>> file, logically R will not know where the be work directory is,
>> correct?
>>
>> So, should I install R in my preferred working directory?
>>
>> Thanks again, in advance. 
>> Bruce
>>
>
>Hm, I think I don't understand what you mean here. Where R
>is installed and where it is run are (usually) not the same
>places.
>
>Let $ be a shell prompt, and let > be the R
>prompt. Then:
>
>  $ cd /tmp/
>  $ R -q
>  > getwd()
>  [1] "/tmp"
>  > q("no")
>  $ cd ~/Documents/
>  $ R -q
>  > getwd()
>  [1] "~/Documents"
>
>Now I write into my "~./Rprofile" file:
>
>  setwd("~/Downloads")
>
>$ cd /tmp/
>$ R -q
>> getwd()
>[1] "~/Downloads"
>
>
>But maybe I am completely misunderstanding what you
>mean
>
>Kind regards
> Enrico
>
>
>>
>>
>>> On Mar 24, 2017, at 3:48 AM, Enrico Schumann 
>wrote:
>>> 
>>> On Thu, 23 Mar 2017, Bruce Ratner PhD writes:
>>> 
 Hi R'ers:
 I would like to setting up a .Rprofile file with
 setwd("C:/R_WorkDir")
 set.seed(12345)
 options (prompt "> R ")
 
 ---
 Can you help providing the code or instructive link,
 I've find many links, but I can't figure it out?
 
 Thanks.
 Bruce 
 
>>> 
>>> Quoting from ?Startup:
>>> 
>>> ,
>>> | [...] unless ‘--no-init-file’ was given, R searches
>>> | for a user profile, a file of R code.  The path of
>>> | this file can be specified by the ‘R_PROFILE_USER’
>>> | environment variable (and tilde expansion will be
>>> | performed).  If this is unset, a file called
>>> | ‘.Rprofile’ is searched for in the current directory
>>> | or in the user's home directory (in that order).  The
>>> | user profile file is sourced into the workspace.
>>> `
>>> 
>
>-- 
>Enrico Schumann
>Lucerne, Switzerland
>http://enricoschumann.net
>
>__
>R-help@r-project.org mailing list -- To UNSUBSCRIBE and more, see
>https://stat.ethz.ch/mailman/listinfo/r-help
>PLEASE do read the posting guide
>http://www.R-project.org/posting-guide.html
>and provide commented, minimal, self-contained, reproducible code.

__
R-help@r-project.org mailing list -- To UNSUBSCRIBE and more, see
https://stat.ethz.ch/mailman/listinfo/r-help
PLEASE do read the posting guide http://www.R-project.org/posting-guide.html
and provide commented, minimal, self-contained, reproducible code.

Re: [R] Error: could not find function "ap"

[PIKAL Petr]

Hi

Did you define function ap?

Cheers
Petr

> -Original Message-
> From: R-help [mailto:r-help-boun...@r-project.org] On Behalf Of Frederic
> Ntirenganya
> Sent: Friday, March 24, 2017 1:18 PM
> To: r-help@r-project.org
> Subject: [R] Error: could not find function "ap"
>
> Dear All,
>
> I hope you are doing well.
>
> I am new in using evapotranspiration packages and got the following error:
>
> Error: could not find function "ap"
>
> The code I am using and data are:
>
> dat1<-read.csv("/home/fredo/Documents/Meteo Data/Meteo Rwa
> data.csv",header=T,na.string="")
> dat1$u.daily<-dat1$Wind.Speed*0.51
> dat1$Date<-as.Date(paste(dat1$Year,dat1$Month,dat1$Day,sep="-"))
> #dat1$DOY<-dayOfYear(dat1$Date)
> dat1$Longitude<-30.11
> dat1$Latitude<--1.95
> dat1$Elevation<-1490
> dat1$RHmax.daily<-dat1$RHmax.daily
> dat1$RHmin.daily<-dat1$RHmin.daily
> dat1$Tmax.daily<-dat1$Tmax.daily
> dat1$Tmin.daily<-dat1$Tmin.daily
> #dat1<-subset(select=-Wind.direction.measured.in.Degrees)
> Sunshine<-dat1$n.daily
> lat<-as.numeric(dat1$Latitude)
> lon<-as.numeric(dat1$Longitude)
> days<-dat1$Date
> require(zoo)
> A <- 0.21
> B <- 0.57
> dat1$RS.daily<-ap(days,lat,lon,A,B,Sunshine,extraT=NULL)
> View(dat1)
>
>
>
> dput(head(dat1))
> structure(list(Year = c(1984L, 1984L, 1984L, 1984L, 1984L, 1984L ), Month =
> c(1L, 1L, 1L, 1L, 1L, 1L), Day = 1:6, Wind.Speed = c(5L, 4L, 4L, 3L, 5L, 6L), 
> n.daily
> = c(6.3, 4.8, 0.6, 8.2, 7.3, 1.7 ), Tmax.daily = c(27.4, 26.3, 22.9, 27.7, 
> 28.5, 25.5),
> Tmin.daily = c(14.5, 16, 14.4, 14.8, 16.6, 15.4), RHmax.daily = c(100L, 95L, 
> 97L,
> 100L, 97L, 99L), RHmin.daily = c(45L, 62L, 72L, 55L, 54L, 63L ), Station.Name 
> =
> structure(c(1L, 1L, 1L, 1L, 1L, 1L), .Label = "KIGALI AERO", class = 
> "factor"),
> Elevation = c(1490, 1490, 1490, 1490, 1490, 1490), Longitude = c(30.11,
> 30.11, 30.11, 30.11, 30.11, 30.11), Latitude = c(-1.95, -1.95,
> -1.95, -1.95, -1.95, -1.95), u.daily = c(2.57222, 2.057776,
> 2.057776, 1.543332, 2.57222, 3.086664), Date = structure(c(5113,
> 5114, 5115, 5116, 5117, 5118), class = "Date")), .Names = c("Year", 
> "Month",
> "Day", "Wind.Speed", "n.daily", "Tmax.daily", "Tmin.daily", "RHmax.daily",
> "RHmin.daily", "Station.Name", "Elevation", "Longitude", "Latitude",
> "u.daily", "Date"), row.names = c(NA, 6L), class = "data.frame")
>
>
> I appreciate your help in advance!
>
> Best regards,
> Fredo.
>
>
>
>
>
>
> Frederic Ntirenganya
> Maseno University,
> African Maths Initiative,
> Kenya.
> Mobile:(+254)718492836
> Email: fr...@aims.ac.za
> https://sites.google.com/a/aims.ac.za/fredo/
>
>   [[alternative HTML version deleted]]
>
> __
> R-help@r-project.org mailing list -- To UNSUBSCRIBE and more, see
> https://stat.ethz.ch/mailman/listinfo/r-help
> PLEASE do read the posting guide http://www.R-project.org/posting-
> guide.html
> and provide commented, minimal, self-contained, reproducible code.


Tento e-mail a jakékoliv k němu připojené dokumenty jsou důvěrné a jsou určeny 
pouze jeho adresátům.
Jestliže jste obdržel(a) tento e-mail omylem, informujte laskavě neprodleně 
jeho odesílatele. Obsah tohoto emailu i s přílohami a jeho kopie vymažte ze 
svého systému.
Nejste-li zamýšleným adresátem tohoto emailu, nejste oprávněni tento email 
jakkoliv užívat, rozšiřovat, kopírovat či zveřejňovat.
Odesílatel e-mailu neodpovídá za eventuální škodu způsobenou modifikacemi či 
zpožděním přenosu e-mailu.

V případě, že je tento e-mail součástí obchodního jednání:
- vyhrazuje si odesílatel právo ukončit kdykoliv jednání o uzavření smlouvy, a 
to z jakéhokoliv důvodu i bez uvedení důvodu.
- a obsahuje-li nabídku, je adresát oprávněn nabídku bezodkladně přijmout; 
Odesílatel tohoto e-mailu (nabídky) vylučuje přijetí nabídky ze strany příjemce 
s dodatkem či odchylkou.
- trvá odesílatel na tom, že příslušná smlouva je uzavřena teprve výslovným 
dosažením shody na všech jejích náležitostech.
- odesílatel tohoto emailu informuje, že není oprávněn uzavírat za společnost 
žádné smlouvy s výjimkou případů, kdy k tomu byl písemně zmocněn nebo písemně 
pověřen a takové pověření nebo plná moc byly adresátovi tohoto emailu případně 
osobě, kterou adresát zastupuje, předloženy nebo jejich existence je adresátovi 
či osobě jím zastoupené známá.

This e-mail and any documents attached to it may be confidential and are 
intended only for its intended recipients.
If you received this e-mail by mistake, please immediately inform its sender. 
Delete the contents of this e-mail with all attachments and its copies from 
your system.
If you are not the intended recipient of this e-mail, you are not authorized to 
use, disseminate, copy or disclose this e-mail in any manner.
The sender of this e-mail shall not be liable for any possible damage caused by 
modifications of the e-mail or by delay with transfer of the email.

In case that this e-mail forms part of business dealings:
- 

[R] Error: could not find function "ap"

[Frederic Ntirenganya]

Dear All,

I hope you are doing well.

I am new in using evapotranspiration packages and got the following error:

Error: could not find function "ap"

The code I am using and data are:

dat1<-read.csv("/home/fredo/Documents/Meteo Data/Meteo Rwa
data.csv",header=T,na.string="")
dat1$u.daily<-dat1$Wind.Speed*0.51
dat1$Date<-as.Date(paste(dat1$Year,dat1$Month,dat1$Day,sep="-"))
#dat1$DOY<-dayOfYear(dat1$Date)
dat1$Longitude<-30.11
dat1$Latitude<--1.95
dat1$Elevation<-1490
dat1$RHmax.daily<-dat1$RHmax.daily
dat1$RHmin.daily<-dat1$RHmin.daily
dat1$Tmax.daily<-dat1$Tmax.daily
dat1$Tmin.daily<-dat1$Tmin.daily
#dat1<-subset(select=-Wind.direction.measured.in.Degrees)
Sunshine<-dat1$n.daily
lat<-as.numeric(dat1$Latitude)
lon<-as.numeric(dat1$Longitude)
days<-dat1$Date
require(zoo)
A <- 0.21
B <- 0.57
dat1$RS.daily<-ap(days,lat,lon,A,B,Sunshine,extraT=NULL)
View(dat1)



dput(head(dat1))
structure(list(Year = c(1984L, 1984L, 1984L, 1984L, 1984L, 1984L
), Month = c(1L, 1L, 1L, 1L, 1L, 1L), Day = 1:6, Wind.Speed = c(5L,
4L, 4L, 3L, 5L, 6L), n.daily = c(6.3, 4.8, 0.6, 8.2, 7.3, 1.7
), Tmax.daily = c(27.4, 26.3, 22.9, 27.7, 28.5, 25.5), Tmin.daily = c(14.5,
16, 14.4, 14.8, 16.6, 15.4), RHmax.daily = c(100L, 95L, 97L,
100L, 97L, 99L), RHmin.daily = c(45L, 62L, 72L, 55L, 54L, 63L
), Station.Name = structure(c(1L, 1L, 1L, 1L, 1L, 1L), .Label = "KIGALI
AERO", class = "factor"),
Elevation = c(1490, 1490, 1490, 1490, 1490, 1490), Longitude = c(30.11,
30.11, 30.11, 30.11, 30.11, 30.11), Latitude = c(-1.95, -1.95,
-1.95, -1.95, -1.95, -1.95), u.daily = c(2.57222, 2.057776,
2.057776, 1.543332, 2.57222, 3.086664), Date = structure(c(5113,
5114, 5115, 5116, 5117, 5118), class = "Date")), .Names = c("Year",
"Month", "Day", "Wind.Speed", "n.daily", "Tmax.daily", "Tmin.daily",
"RHmax.daily", "RHmin.daily", "Station.Name", "Elevation", "Longitude",
"Latitude", "u.daily", "Date"), row.names = c(NA, 6L), class = "data.frame")


I appreciate your help in advance!

Best regards,
Fredo.






Frederic Ntirenganya
Maseno University,
African Maths Initiative,
Kenya.
Mobile:(+254)718492836
Email: fr...@aims.ac.za
https://sites.google.com/a/aims.ac.za/fredo/

[[alternative HTML version deleted]]

__
R-help@r-project.org mailing list -- To UNSUBSCRIBE and more, see
https://stat.ethz.ch/mailman/listinfo/r-help
PLEASE do read the posting guide http://www.R-project.org/posting-guide.html
and provide commented, minimal, self-contained, reproducible code.

Re: [R] Setting up a .Rprofile file

[Enrico Schumann]

On Fri, 24 Mar 2017, Bruce Ratner PhD writes:

> Henrico:
> Thanks for quick reply.
> However, one last question:
> If I want to change working directory, and put setwd() in the Rprofile
> file, logically R will not know where the be work directory is,
> correct?
>
> So, should I install R in my preferred working directory?
>
> Thanks again, in advance. 
> Bruce
>

Hm, I think I don't understand what you mean here. Where R
is installed and where it is run are (usually) not the same
places.

Let $ be a shell prompt, and let > be the R
prompt. Then:

  $ cd /tmp/
  $ R -q
  > getwd()
  [1] "/tmp"
  > q("no")
  $ cd ~/Documents/
  $ R -q
  > getwd()
  [1] "~/Documents"

Now I write into my "~./Rprofile" file:

  setwd("~/Downloads")

$ cd /tmp/
$ R -q
> getwd()
[1] "~/Downloads"


But maybe I am completely misunderstanding what you
mean

Kind regards
 Enrico


>
>
>> On Mar 24, 2017, at 3:48 AM, Enrico Schumann  wrote:
>> 
>> On Thu, 23 Mar 2017, Bruce Ratner PhD writes:
>> 
>>> Hi R'ers:
>>> I would like to setting up a .Rprofile file with
>>> setwd("C:/R_WorkDir")
>>> set.seed(12345)
>>> options (prompt "> R ")
>>> 
>>> ---
>>> Can you help providing the code or instructive link,
>>> I've find many links, but I can't figure it out?
>>> 
>>> Thanks.
>>> Bruce 
>>> 
>> 
>> Quoting from ?Startup:
>> 
>> ,
>> | [...] unless ‘--no-init-file’ was given, R searches
>> | for a user profile, a file of R code.  The path of
>> | this file can be specified by the ‘R_PROFILE_USER’
>> | environment variable (and tilde expansion will be
>> | performed).  If this is unset, a file called
>> | ‘.Rprofile’ is searched for in the current directory
>> | or in the user's home directory (in that order).  The
>> | user profile file is sourced into the workspace.
>> `
>> 

-- 
Enrico Schumann
Lucerne, Switzerland
http://enricoschumann.net

__
R-help@r-project.org mailing list -- To UNSUBSCRIBE and more, see
https://stat.ethz.ch/mailman/listinfo/r-help
PLEASE do read the posting guide http://www.R-project.org/posting-guide.html
and provide commented, minimal, self-contained, reproducible code.

Re: [R] Modeling Time Series with Missing Observations

[Berend Hasselman]

There is also a package imputeTS, which may be able to do what you want.
It has a nice Introduction vignette  and also appears to have nice plot 
functions

Berend Hasselman


> On 24 Mar 2017, at 02:11, David Winsemius  wrote:
> 
> There's also an irts-package "irregular time series"
> 
> Sent from my iPhone
> 
>> On Mar 23, 2017, at 3:42 PM, John C Frain  wrote:
>> 
>> Google "arima missing data r"  will bring up several references including
>> http://stats.stackexchange.com/questions/104565/how-to-use-auto-arima-to-impute-missing-values.
>> There are several other useful results in that search.
>> 
>> John C Frain
>> 3 Aranleigh Park
>> Rathfarnham
>> Dublin 14
>> Ireland
>> www.tcd.ie/Economics/staff/frainj/home.html
>> mailto:fra...@tcd.ie
>> mailto:fra...@gmail.com
>> 
>>> On 23 March 2017 at 15:42, Jeff Newmiller  wrote:
>>> 
>>> Even the most basic introduction to R discusses the use of NA for missing
>>> data. Injecting values that could be mistaken for actual readings is a
>>> dangerous approach. You can use the merge function to introduce missing
>>> rows into zoo objects or data frames.
>>> --
>>> Sent from my phone. Please excuse my brevity.
>>> 
>>> On March 23, 2017 8:22:47 AM PDT, Paul Bernal 
>>> wrote:
 Dear all,
 
 Hope you are doing well. I am trying to model the historical number of
 transits of a particular market segment, but the problem is that I have
 missing data.
 
 I am working with monthly data, so I have 12 observations per year (in
 general). The problem is that, when I bring the data from the database,
 the
 following happens, for example:
 
 January-2000, Feb-2000, Apr-2000, Jun 2000 (I have missing
 observations)
 
 when I am supposed to have the sequence January-2000, Feb-2000,
 Mar-2000,
 Apr-2000, May-2000, Jun-2000, etc.
 
 How can I model a time series when there are missing months? I was
 planning
 making up fictional or fake observations with a value of 1 to fill in
 the
 gaps but not sure if this is a reasonable approach.
 
 Any help and/or guidance will be greatly appreciated,
 
 Best regards,
 
 Paul
 
 [[alternative HTML version deleted]]
 
 __
 R-help@r-project.org mailing list -- To UNSUBSCRIBE and more, see
 https://stat.ethz.ch/mailman/listinfo/r-help
 PLEASE do read the posting guide
 http://www.R-project.org/posting-guide.html
 and provide commented, minimal, self-contained, reproducible code.
>>> 
>>> __
>>> R-help@r-project.org mailing list -- To UNSUBSCRIBE and more, see
>>> https://stat.ethz.ch/mailman/listinfo/r-help
>>> PLEASE do read the posting guide http://www.R-project.org/
>>> posting-guide.html
>>> and provide commented, minimal, self-contained, reproducible code.
>>> 
>> 
>>   [[alternative HTML version deleted]]
>> 
>> __
>> R-help@r-project.org mailing list -- To UNSUBSCRIBE and more, see
>> https://stat.ethz.ch/mailman/listinfo/r-help
>> PLEASE do read the posting guide http://www.R-project.org/posting-guide.html
>> and provide commented, minimal, self-contained, reproducible code.
> 
> __
> R-help@r-project.org mailing list -- To UNSUBSCRIBE and more, see
> https://stat.ethz.ch/mailman/listinfo/r-help
> PLEASE do read the posting guide http://www.R-project.org/posting-guide.html
> and provide commented, minimal, self-contained, reproducible code.

__
R-help@r-project.org mailing list -- To UNSUBSCRIBE and more, see
https://stat.ethz.ch/mailman/listinfo/r-help
PLEASE do read the posting guide http://www.R-project.org/posting-guide.html
and provide commented, minimal, self-contained, reproducible code.

Re: [R] Setting up a .Rprofile file

[Enrico Schumann]

On Thu, 23 Mar 2017, Bruce Ratner PhD writes:

> Hi R'ers:
> I would like to setting up a .Rprofile file with
> setwd("C:/R_WorkDir")
> set.seed(12345)
> options (prompt "> R ")
>
> ---
> Can you help providing the code or instructive link,
> I've find many links, but I can't figure it out?
>
> Thanks.
> Bruce 
>

Quoting from ?Startup:

,
| [...] unless ‘--no-init-file’ was given, R searches
| for a user profile, a file of R code.  The path of
| this file can be specified by the ‘R_PROFILE_USER’
| environment variable (and tilde expansion will be
| performed).  If this is unset, a file called
| ‘.Rprofile’ is searched for in the current directory
| or in the user's home directory (in that order).  The
| user profile file is sourced into the workspace.
`

-- 
Enrico Schumann
Lucerne, Switzerland
http://enricoschumann.net

__
R-help@r-project.org mailing list -- To UNSUBSCRIBE and more, see
https://stat.ethz.ch/mailman/listinfo/r-help
PLEASE do read the posting guide http://www.R-project.org/posting-guide.html
and provide commented, minimal, self-contained, reproducible code.

[R] Setting up a .Rprofile file

[Bruce Ratner PhD]

Hi R'ers:
I would like to setting up a .Rprofile file with
setwd("C:/R_WorkDir")
set.seed(12345)
options (prompt "> R ")

---
Can you help providing the code or instructive link, I've find many links, but 
I can't figure it out?

Thanks.
Bruce 

__
Bruce Ratner PhD
The Significant Statistician™
(516) 791-3544
Statistical Predictive Analytics -- www.DMSTAT1.com
Machine-Learning Data Mining -- www.GenIQ.net

__
R-help@r-project.org mailing list -- To UNSUBSCRIBE and more, see
https://stat.ethz.ch/mailman/listinfo/r-help
PLEASE do read the posting guide http://www.R-project.org/posting-guide.html
and provide commented, minimal, self-contained, reproducible code.

Re: [R] (no subject)

[Bert Gunter]

Please read and follow the directions at the bottom of your email.

-- 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 County" comic strip )


On Thu, Mar 23, 2017 at 9:55 AM, munevver kaya  wrote:
> Dear Sir/Madam,
> I do not want to receive any email by r-help. Can you delete me in your
> mail list? Thank you.
>
> Sincerely,
> Munevver Basman
>
> [[alternative HTML version deleted]]
>
> __
> R-help@r-project.org mailing list -- To UNSUBSCRIBE and more, see
> https://stat.ethz.ch/mailman/listinfo/r-help
> PLEASE do read the posting guide http://www.R-project.org/posting-guide.html
> and provide commented, minimal, self-contained, reproducible code.

__
R-help@r-project.org mailing list -- To UNSUBSCRIBE and more, see
https://stat.ethz.ch/mailman/listinfo/r-help
PLEASE do read the posting guide http://www.R-project.org/posting-guide.html
and provide commented, minimal, self-contained, reproducible code.