search for: macqueen

...moving average data. To see this, break your work into two separate steps, like this: tnr.ma <- ma(dat3[1:28], order=3) TNR_moving_average <- forecast(tnr.ma, h=8) I think you will find that the warning comes from the second step. Print tnr.ma and you will see some NAs. -Don -- Don MacQueen Lawrence Livermore National Laboratory 7000 East Ave., L-627 Livermore, CA 94550 925-423-1062 Lab cell 925-724-7509 ? From: Bill Poling <Bill.Poling at zelis.com> Date: Friday, June 1, 2018 at 8:58 AM To: "MacQueen, Don" <macqueen1 at llnl.gov>, array R-help <r-help at r-...

...ue, say the average of the others? Also, I have 5 series of data I am working with using this script and of course each gave me that warning, but only the one series had the same 8 Point Forecast values, is that coincidental you think? Terrific of you to help, I really appreciate it. WHP From: MacQueen, Don [mailto:macqueen1 at llnl.gov] Sent: Friday, June 01, 2018 12:54 PM To: Bill Poling <Bill.Poling at zelis.com>; r-help (r-help at r-project.org) <r-help at r-project.org> Subject: Re: [R] Time-series moving average question You are right that there are no NAs in the practice data....

...https://cran.r-project.org/web/packages/forecast/forecast.pdf Since I created this practice data using the structure() function I am unsure why there would be NA?s as there are none apparently in the structure? No worries though, I am going to reach out to the package author. Cheers. WHP From: MacQueen, Don [mailto:macqueen1 at llnl.gov] Sent: Friday, June 01, 2018 11:24 AM To: Bill Poling <Bill.Poling at zelis.com>; r-help (r-help at r-project.org) <r-help at r-project.org> Subject: Re: [R] Time-series moving average question My guess would be that if you inspect the output from ma(...

...duced calculation time a lot by extracting the numeric columns of a data frame and working with them, then recombining them with the character columns. R?s performance working with data frames has improved since then, so I really don?t know if it would make a difference for your task. -Don -- Don MacQueen Lawrence Livermore National Laboratory 7000 East Ave., L-627 Livermore, CA 94550 925-423-1062 Lab cell 925-724-7509 From: Luca Meyer <lucam1968 at gmail.com> Date: Monday, April 30, 2018 at 8:08 AM To: Rui Barradas <ruipbarradas at sapo.pt> Cc: "MacQueen, Don" <macqueen1...

Is there a package or a command that does window aggregation like select sum(col1) over (partition by col2, col3 order by col4 rows between unbounded preceding and current row) as sum1 from table1 ; the above is Netezza syntax, but Postgre has same capability. Stephen B [[alternative HTML version deleted]]

...int, and jump to wherever the hypelink pointed. I don't remember then name of the package or function, and haven't been able to find it after some searching. Does anyone know the name, or have a suggestion to help track it down? Thanks -Don -- -------------------------------------- Don MacQueen Environmental Protection Department Lawrence Livermore National Laboratory Livermore, CA, USA 925-423-1062

Thanks. I thought lm() function is for linear model, such as the correlation below: Y= aX + b On Wed, Jun 14, 2017 at 5:25 PM, MacQueen, Don <macqueen1 at llnl.gov> wrote: > Start with the lm() function; i.e., see > > ?lm > > -Don > > -- > Don MacQueen > > Lawrence Livermore National Laboratory > 7000 East Ave., L-627 > Livermore, CA 94550 > 925-423-1062 > > > On 6/14/17, 3:4...

...ut from ma(dat3[1:28], order=3) you will find some NAs in it. And then forecast() doesn't like NAs. But I can't check, because I can't find the ma() and forecast() functions. I assume they come from some package you installed; it would be helpful to say which package. -Don -- Don MacQueen Lawrence Livermore National Laboratory 7000 East Ave., L-627 Livermore, CA 94550 925-423-1062 Lab cell 925-724-7509 ?On 6/1/18, 4:20 AM, "R-help on behalf of Bill Poling" <r-help-bounces at r-project.org on behalf of Bill.Poling at zelis.com> wrote: Good morning, I hope so...

Hi I have a vector like that readCsvFile$V1  [1]  30  31  32  33  34  35  36  37  38  39 310 311 312 313 314 315 316 317 318 [20] 319 320 321  20  21  22  23  24  25  26  27  28  29 210 211 212 213 214 215 [39] 216 217 218 219 220 221 222 223  40  41  42  43  44  45  46  47  48  49 410 [58] 411 412 413 414 415 416 417 418 419 420 421 and I am looking to find where the number 31 is located. So I

Does R have a comparison operator similar to the Like function, for example: a<-"Is a Fish" b<-"Fish" if(b in a){c<-TRUE} Michael R Howard Micron Technology Inc. Boise ID. Fab C Engineering Software (FCES) Software Engineer

...comment that regexpr might indeed be the time consuming piece of code. I might try to optimise this piece of code later on, but for the time being I am working on the following part of building a neural network to try indeed classifying some text. Again, thanks, Luca 2018-04-30 17:25 GMT+02:00 MacQueen, Don <macqueen1 at llnl.gov>: > Luca, > > > > If speed is important, you might improve performance by making d0 into a > true matrix, rather than a data frame (assuming d0 is indeed a data frame > at this point). Although data frames may look like matrices, they aren?t,...

...rectly, 2001 is prior to the current R support for locales and extended character sets. Using \265 is what I could find at that time to get a mu into my output. I came across this while checking some things; it?s not actually breaking my scripts, so I doubt it?s due to any recent change. -- Don MacQueen Lawrence Livermore National Laboratory 7000 East Ave., L-627 Livermore, CA 94550 925-423-1062

...er if you run data(tritest) x <- tritest$x y <- tritest$y vm <- voronoi.mosaic(x,y) plot(vm) par(new=T) plot(x,y,col='blue') it looks at it should. -zubin > > From: Prof Brian Ripley <ripley at stats.ox.ac.uk> > Date: 2006/07/19 Wed AM 03:15:51 EDT > To: Don MacQueen <macq at llnl.gov> > CC: zubin <binabina at bellsouth.net>, r-help at stat.math.ethz.ch > Subject: Re: [R] voronoi tessellations > > On Tue, 18 Jul 2006, Don MacQueen wrote: > > > I'll suggest going to the CRAN packages page and doing a search for "voron...

...6)) axis.Date(side=3, at=seq(min(LAI_simulation$Date), max(LAI_simulation$Date),by="year"), format="%Y", xlab="xxx",cex.axis=0.9) mtext("years",side=3,line=3,cex.lab=1.2) Ahmed Attia, Ph.D. Agronomist & Soil Scientist On Wed, Jul 5, 2017 at 4:52 PM, MacQueen, Don <macqueen1 at llnl.gov> wrote: > Not entirely sure what you really want the axis on the top to look like, but here's what I'd be likely to do: > > ## note, correction to format; in your example data year comes last > LAI_simulation$Date <- as.Date( LAI_simulation$D...

...; wrote: > Hi > > But X can be some function like - sin, cos, log, exp... > > Cheers > Petr > > > -----Original Message----- > > From: R-help [mailto:r-help-bounces at r-project.org] On Behalf Of lily li > > Sent: Thursday, June 15, 2017 1:28 AM > > To: MacQueen, Don <macqueen1 at llnl.gov> > > Cc: R mailing list <r-help at r-project.org> > > Subject: Re: [R] about fitting a regression line > > > > Thanks. I thought lm() function is for linear model, such as the > correlation > > below: > > Y= aX + b >...

I use kmeans to classify spectral events in high and low 1/3 octave bands: #Do cluster analysis CyclA<-data.frame(LlowA,LhghA) CntrA<-matrix(c(0.9,0.8,0.8,0.75,0.65,0.65), nrow = 3, ncol=2, byrow=TRUE) ClstA<-kmeans(CyclA,centers=CntrA,nstart=50,algorithm="MacQueen") This works well when the actual data shows 1,2 or 3 groups that are not "too close" in a cross plot. The MacQueen algorithm will give one or more empty groups which is what I want. However, there are cases when the groups are closer together, less compact or diffuse which leads t...

...quot;, newtst[i + 1])){ newtst_new[c(j, j + 1)] <- newtst[c(i, i + 1)] i <- i + 2 }else{ newtst_new[c(j, j + 1)] <- c(newtst[c(i)], "Fval: ") i <- i + 1 } j <- j + 2 } newtst_new which is also not very pretty. HTH Ulrik On Thu, 13 Jul 2017 at 16:48 MacQueen, Don <macqueen1 at llnl.gov> wrote: > Using Ulrik?s example data (and assuming I understand what is wanted), > here is what I would do: > > ex.dat <- c("FName: fname1", "Fval: Fval1.name1", "Fval: ", "FName: > fname2", "Fval: Fva...

...2.48" > fmt.signif(x[3],3) [1] "0.120" > > signif(x,3) [1] "2.50" "2.48" "0.120" For now, at least, I'm willing to assume that I never want the numbers formatted with exponential notation. -- -------------------------------------- Don MacQueen Environmental Protection Department Lawrence Livermore National Laboratory Livermore, CA, USA --------------------------------------

...called, file3 is deleted, file2 is renamed file3, file1 is renamed file2, and the name of file1 is returned to the user's script, which then writes to file1. I've done some searching (RSiteSearch) but haven't found anything. Thanks -Don -- -------------------------------------- Don MacQueen Environmental Protection Department Lawrence Livermore National Laboratory Livermore, CA, USA 925-423-1062

Start with the lm() function; i.e., see ?lm -Don -- Don MacQueen Lawrence Livermore National Laboratory 7000 East Ave., L-627 Livermore, CA 94550 925-423-1062 On 6/14/17, 3:40 PM, "R-help on behalf of lily li" <r-help-bounces at r-project.org on behalf of chocold12 at gmail.com> wrote: Hi R users, I have some data points (Xi, Yi)...