Re: [R] to R Core T: mle function in 32bits not respecting the constrain

2018-05-28 Thread francesc badia roca
The have set an epsilon > 0 to the point where the assimptota is, to compute the hessian at the lower it is required to evaluate at the left of the lower, and this is the point, because I want the estimate that usually is very close to the lower. I want to set a lower bound in which neither the est

Re: [R] help on the R package Softmaxreg

2018-05-28 Thread Vy Ong
Many thanks, Bert & Rolf. I have read the package manual and its github code to try understanding the things. Bests, V.O. On Mon, May 28, 2018 at 2:42 AM, Bert Gunter wrote: > ... who can be found by R's ?maintainer function. > > A web search on "R softmaxreg package" also seemed to bring up rel

Re: [R] Difficulty in writing R code for one pool dynamic model

2018-05-28 Thread David Winsemius
> On May 28, 2018, at 5:02 PM, Kebebe Ergano via R-help > wrote: > > > Hi everyone, > > I was trying to mode the following exercise using R. > > The question: Set up a one pool model using numericintegration. The model > will run from time 1 to time 30 using a time step of 1.The pool (A)

Re: [R] How to generate a conditional dummy in R?

2018-05-28 Thread Jim Lemon
Hi Faradj, What a problem! I think I have worked it out, but only because the result is the one you said you wanted. # the sample data frame is named fkdf Y2Xby3<-function(x) { nrows<-dim(x)[1] X<-rep(0,nrows) for(i in 1:(nrows-2)) { if(!is.na(x$Y[i])) { if(x$Y[i] == 1 && any(is.na(x$Y[(i+

[R] Difficulty in writing R code for one pool dynamic model

2018-05-28 Thread Kebebe Ergano via R-help
Hi everyone, I was trying to mode the following exercise using R. The question: Set up a one pool model using numericintegration. The model will run from time 1 to time 30 using a time step of 1.The pool (A) will be fed by flux "inA" at a rate of 5 units per hour anddrained by flux "outA" at

Re: [R] How to generate a conditional dummy in R?

2018-05-28 Thread Bert Gunter
I was not able to decipher your meaning, and your failure to garner a response so far may mean that others may have had similar difficultes. You might therefore try providing providing a **small reproducible** example of X's and Y's that show what you have and what you want to end up with to clarif

Re: [R] to R Core T: mle function in 32bits not respecting the constrain

2018-05-28 Thread Berry, Charles
> On May 27, 2018, at 10:31 PM, francesc badia roca wrote: > > I have an issue using mle in versions of 32 bits. > > I am writing a package which I want to submit to the CRAN. > When doing the check, there is an example that has an error running in the > 32 bits version. > > The problem comes

[R] How to generate a conditional dummy in R?

2018-05-28 Thread Faradj Koliev
Hi everyone, I am trying to generate a conditional dummy variable ”X" with the following rules set X=1 if Y is =1, two years prior to the NA. [0,0,NA]. For example, if the pattern for Y is 0,0,NA then the X variable is =0 for all the two years prior to the NA. If the pattern for Y is 0,1

Re: [R] help on the R package Softmaxreg

2018-05-28 Thread Bert Gunter
... who can be found by R's ?maintainer function. A web search on "R softmaxreg package" also seemed to bring up relevant info. -- 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 "Bloo

Re: [R] help on the R package Softmaxreg

2018-05-28 Thread Rolf Turner
On 28/05/18 18:13, Vy Ong wrote: Hello, Are there anyone knowing about the Softmaxreg R package? I suspect that the maintainer of softmaxreg (please note the *lower case* "s") "knows about" this package! cheers, Rolf Turner -- Technical Editor ANZJS Department of Statistics University of

[R] help on the R package Softmaxreg

2018-05-28 Thread Vy Ong
Hello, Are there anyone knowing about the Softmaxreg R package? I try to understand the mathematics behind its models, by its command softmaxreg, especially the L2 regularized version. Is it regularized on the neural network layers or on the softmax regression layer? Are there any material/article

Re: [R] TukeyHSD for multiple response

2018-05-28 Thread Sergio Ferreira Cardoso
Michael, Thank you very much for your answer. I finally tried lsmeans to compare what I wanted. I'll follow your advice and explore the CDA. It's probably a better solution to assess what I want. Best, Sérgio. - Mensagem original - > De: "Michael Friendly" > Para: "Sergio Ferreira Ca