Re: [R] How to generate SE for the proportion value using a randomization process in R?

2021-01-28 Thread Rui Barradas
Hello, Yes, sorry for my previous post, I had forgotten about boot.array. That's a much better solution for your problem. Rui Barradas Às 20:29 de 28/01/21, Marna Wagley escreveu: Thank you Rui, This is great. How about the following? SimilatedData<-boot.array(b, indices=T) seems it is givin

Re: [R] How to generate SE for the proportion value using a randomization process in R?

2021-01-28 Thread Marna Wagley
Thank you Rui, This is great. How about the following? SimilatedData<-boot.array(b, indices=T) seems it is giving the rows ID which are used in the calculation, isn't it? On Thu, Jan 28, 2021 at 12:21 PM Rui Barradas wrote: > Hello, > > I don't know why you would need to see the indices but

Re: [R] How to generate SE for the proportion value using a randomization process in R?

2021-01-28 Thread Rui Barradas
Hello, I don't know why you would need to see the indices but rewrite the function bootprop as bootprop_ind <- function(data, index){ d <- data[index, ] #sum(d[["BothTimes"]], na.rm = TRUE)/sum(d[["Time1"]], na.rm = TRUE) index } and call in the same way. It will now return a matrix of

Re: [R] random numbers with constraints

2021-01-28 Thread Abby Spurdle
I recognize the problems with global data. And my code could certainly be improved. However, I also note that the random numbers (ignoring transformations), need to be constant, while computing the rate. Otherwise, my algorithm wouldn't work well. As it is, rounding operations can cause "jumps".

Re: [R] How to generate SE for the proportion value using a randomization process in R?

2021-01-28 Thread Marna Wagley
Hi Rui, I am sorry for asking you several questions. In the given example, randomizations (reshuffle) were done 1000 times, and its 1000 proportion values (results) are stored and it can be seen using b$t; but I was wondering how the table was randomized (which rows have been missed/or repeated in

Re: [R] Regarding fitted value

2021-01-28 Thread Rui Barradas
Hello, From help('forecast::fitted.Arima'): h The number of steps to forecast ahead. So you have the default h = 1 step ahead forecast for your model. Hope this helps, Rui Barradas Às 12:13 de 28/01/21, Md. Moyazzem Hossain escreveu: Dear R-experts, I hope that all of you are doi

Re: [R] random numbers with constraints

2021-01-28 Thread Martin Maechler
> Abby Spurdle > on Thu, 28 Jan 2021 08:48:06 +1300 writes: > I note that there's a possibility of floating point errors. > If all values have one digit after the decimal point, you could replace: > qexp (p, rate) with round (qexp (p, rate), 1). > However, sometimes u

Re: [R] Citing

2021-01-28 Thread Bert Gunter
?citation 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, Jan 28, 2021 at 4:37 AM M Visconti wrote: > Hi there, > > Our statistician used R version 4.0.3

Re: [R] Citing

2021-01-28 Thread Patrick (Malone Quantitative)
When you start up R, one of the opening messages is: 'citation()' on how to cite R or R packages in publications. Since it sounds like you're not using R yourself, here are the results for R, but you'd have to go one by one through packages. R Core Team (2020). R: A language and environment for

[R] Citing

2021-01-28 Thread M Visconti
Hi there, Our statistician used R version 4.0.3 for a clinical project data analysis. I wanted to cite correctly using this below format (for a different example): NCSS 10 Statistical Software [2015]; NCSS, LLC, Kaysville, UT Can you please help me fill in the info after the semi colon for your

[R] Regarding fitted value

2021-01-28 Thread Md. Moyazzem Hossain
Dear R-experts, I hope that all of you are doing well. I got the filled value from the ARIMA model. I use the following working code. But I am not clear whether I got the fitted value for each *corresponding time* of the original data point like 2000, 2001, 2020 or get a *one-step-ahead* fitted

Re: [R] random numbers with constraints

2021-01-28 Thread Denis Francisci
Thanks again for your help, One digit after the decimal point is enough for my purposes; so, I can round the qexp function, even if possible errors in floating points are not a problem. Thank you very very much, Denis Il giorno mer 27 gen 2021 alle ore 20:48 Abby Spurdle ha scritto: > I note