Re: [R] Smallest Space Analysis (SSA) in R

2018-09-20 Thread Martin Maechler
> Eric Berger > on Thu, 20 Sep 2018 23:28:27 +0300 writes: > Hi Andrew, > I don't have any experience in this area but I was intrigued by your > question. Here is what I learned. > 1, A bit of poking around turned up a thread on stats.stackexchange that > mentions

[R] Fitting Production Curves

2018-09-20 Thread mikorym via R-help
Hi All, By a production curve I mean for example the output of a mine, peak oil production or the yield of a farm over time within the same season. It is this last example that we should take as the prototypical case. What I would like to do is to fit a curve that inherits qualities of the dis

Re: [R] How to vectorize this function

2018-09-20 Thread Lynette Chang
Here are the code and data file. I’m not sure if I put too much unrelated information here. My goal is to factor out volatilities from the data. I hope I can get sigV <- impVolC(callM, K, T, F, r), which has five vectors as input, and one vector as output. The length of all those six vectors a

Re: [R] tibble question with a mean

2018-09-20 Thread Erin Hodgess
David That's awesome! Thank you!!! Erin Hodgess, PhD mailto: erinm.hodg...@gmail.com On Thu, Sep 20, 2018 at 9:19 PM David L Carlson wrote: > > xt[, 2:3] %>% colMeans > y z > 2.500 -0.4401625 > > > xt[2] %>% colMeans > y > 2.5 > > t(xt[, 2]) %>% mean > [1] 2.5 > > ---

Re: [R] tibble question with a mean

2018-09-20 Thread David L Carlson
> xt[, 2:3] %>% colMeans y z 2.500 -0.4401625 > xt[2] %>% colMeans y 2.5 > t(xt[, 2]) %>% mean [1] 2.5 - David L. Carlson Department of Anthropology Texas A&M University -Original Message- From: R-help [mailto:r-help-boun...@r-project.o

Re: [R] [FORGED] Re: tibble question with a mean

2018-09-20 Thread Erin Hodgess
Thanks to all for the good suggestions!!! Erin Hodgess, PhD mailto: erinm.hodg...@gmail.com On Thu, Sep 20, 2018 at 7:38 PM Rolf Turner wrote: > > Please see inline below. > > n 09/21/2018 01:07 PM, Peter Langfelder wrote: > > I don't know tibble, so I'll do the same with a plain data frame: >

Re: [R] [FORGED] Re: tibble question with a mean

2018-09-20 Thread Rolf Turner
Please see inline below. n 09/21/2018 01:07 PM, Peter Langfelder wrote: I don't know tibble, so I'll do the same with a plain data frame: a = data.frame(x=LETTERS[1:4],y=1:4,z=rnorm(4),a=c("dog","cat","tree","ferret")) a x y z a 1 A 1 -0.08264865dog 2 B 2 0.32344426

Re: [R] tibble question with a mean

2018-09-20 Thread Peter Langfelder
I don't know tibble, so I'll do the same with a plain data frame: a = data.frame(x=LETTERS[1:4],y=1:4,z=rnorm(4),a=c("dog","cat","tree","ferret")) > a x y z a 1 A 1 -0.08264865dog 2 B 2 0.32344426cat 3 C 3 -0.80416061 tree 4 D 4 1.27052529 ferret > mean(a[2:3]) [1] NA

[R] tibble question with a mean

2018-09-20 Thread Erin Hodgess
Hello! Here is a toy tibble problem: xt <- tibble(x=LETTERS[1:4],y=1:4,z=rnorm(4),a=c("dog","cat","tree","ferret")) str(xt) Classes ‘tbl_df’, ‘tbl’ and 'data.frame': 4 obs. of 4 variables: $ x: chr "A" "B" "C" "D" $ y: int 1 2 3 4 $ z: num 0.3246 0.0504 0.339 0.4872 $ a: chr "dog" "cat"

Re: [R] How to vectorize this function

2018-09-20 Thread Jeff Newmiller
Sorry, misread your comment, I agree. 4/2 has one arithmetic operation, (1/2)*4 has two to accomplish the same calculation. On September 20, 2018 1:53:03 PM PDT, "MacQueen, Don" wrote: >You're asking me? > >I prefer >> 4/2 >[1] 2 > >not >> 1/2*4 >[1] 2 > >(I think that's what I said) > >And if

Re: [R] How to vectorize this function

2018-09-20 Thread MacQueen, Don via R-help
You're asking me? I prefer > 4/2 [1] 2 not > 1/2*4 [1] 2 (I think that's what I said) And if I did want to multiply by the reciprocal, which does happen from time to time, I'd certainly do it this way: (1/2)*4 -Don -- Don MacQueen Lawrence Livermore National Laboratory 7000 East Ave., L-6

Re: [R] Smallest Space Analysis (SSA) in R

2018-09-20 Thread Andrew
Hi Eric I will need to dig into this a bit deeper, but this looks like it might hold some promise. The web link you shared seems familiar - perhaps I came across it but not at the site you linked to. I will read the sources with interest. Thank you for bringing them to my attention. Regards,

Re: [R] Make Month variable be called “May”,“June” etc., instead of a numeric quantity 5,6, etc.

2018-09-20 Thread AbouEl-Makarim Aboueissa
Dear Jeff: thank you very much abou __ *AbouEl-Makarim Aboueissa, PhD* *Professor of Statistics* *Graduate Coordinator* *Department of Mathematics and Statistics* *University of Southern Maine* On Thu, Sep 20, 2018 at 1:37 PM Jeff Newmiller wrote: > airquality$Month

Re: [R] Smallest Space Analysis (SSA) in R

2018-09-20 Thread Eric Berger
Hi Andrew, I don't have any experience in this area but I was intrigued by your question. Here is what I learned. 1, A bit of poking around turned up a thread on stats.stackexchange that mentions that "smallest space analysis" (SSA) is a special case of "multidimensional scaling" (MDS). See the th

Re: [R] How to vectorize this function

2018-09-20 Thread Jeff Newmiller
re: your last comment... why do you prefer to multiply by the reciprocal? On September 20, 2018 10:56:22 AM PDT, "MacQueen, Don via R-help" wrote: >In addition to what the other said, if callM is a vector then an >expression of the form > if (callM <= call0) >is inappropriate. Objects inside t

Re: [R] How to vectorize this function

2018-09-20 Thread MacQueen, Don via R-help
In addition to what the other said, if callM is a vector then an expression of the form if (callM <= call0) is inappropriate. Objects inside the parentheses of if() should have length one. For example, > if (1:5 < 3) 'a' else 'b' [1] "a" Warning message: In if (1:5 < 3) "a" else "b" : th

Re: [R] Make Month variable be called “May”,“June” etc., instead of a numeric quantity 5,6, etc.

2018-09-20 Thread Jeff Newmiller
airquality$Month <- factor(airquality$Month, levels=1:12, labels=month.name ) On September 20, 2018 10:14:41 AM PDT, AbouEl-Makarim Aboueissa wrote: >Dear All: > >*Re:* How to make the Month variable be called “May”,“June”, "July", >"August", "September" instead of a numeric quantity (5,6,7,8,9)

[R] Make Month variable be called “May”,“June” etc., instead of a numeric quantity 5,6, etc.

2018-09-20 Thread AbouEl-Makarim Aboueissa
Dear All: *Re:* How to make the Month variable be called “May”,“June”, "July", "August", "September" instead of a numeric quantity (5,6,7,8,9) In the airquality data set, please see the code below; How to make the Month variable be called “May”,“June”, "July", "August", "September" instead of a

Re: [R] How to vectorize this function

2018-09-20 Thread Bert Gunter
Also: What package does polya() come from and "gamma" (as a numeric value)is undefined (it is a function). 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

Re: [R] How to vectorize this function

2018-09-20 Thread David L Carlson
Your function takes an argument "F" that is never used and uses an object "y" which is not defined. Give us some data to use for testing different approaches along with the answer you expect. It may be possible to use two ifelse() functions instead of the loop. -

[R] How to vectorize this function

2018-09-20 Thread Lynette Chang
Hello everyone, I’ve a function with five input argument and one output number. impVolC <- function(callM, K, T, F, r) I hope this function can take five vectors as input, then return one vector as output. My vectorization ran into problems with the nested if-else operation

Re: [R] about the series of numbers

2018-09-20 Thread PIKAL Petr
Hi > -Original Message- > From: R-help On Behalf Of Bogdan Tanasa > Sent: Thursday, September 20, 2018 9:26 AM > To: r-help > Subject: [R] about the series of numbers > > Dear all, > > if I may ask please a question that is likely very naive : > > shall I write in R > "1:9", it will gene

Re: [R] about the series of numbers

2018-09-20 Thread Albrecht Kauffmann
Dear Bogdan, for increments <>1, it needs "by". Try seq(.1,.9,.1) Best, Albrecht -- Albrecht Kauffmann alkau...@fastmail.fm Am Do, 20. Sep 2018, um 09:26, schrieb Bogdan Tanasa: > Dear all, > > if I may ask please a question that is likely very naive : > > shall I write in R > "1:9", it

[R] about the series of numbers

2018-09-20 Thread Bogdan Tanasa
Dear all, if I may ask please a question that is likely very naive : shall I write in R > "1:9", it will generate "1 2 3 4 5 6 7 8 9" shall I write > "0.1:0.9", why does it generate only 0.1 ? thank you ! -- bogdan [[alternative HTML version deleted]]