search for: metstat

...d a function that I use quit frequently. I was wondering if there is away to save this function in an existing library so I can call it by using the function name once the library is loaded? Thanks, Doug -- --------------------------------- Douglas M. Hultstrand, MS Senior Hydrometeorologist Metstat, Inc. Windsor, Colorado voice: 970.686.1253 email: dmhultst at metstat.com web: http://www.metstat.com

...tions are appreciated. Thanks, Doug # Example for (i in 1:2) { x=c(1,2,3,4) y=c(10,20,30,40) G <- paste("name", i, sep="") G[i] <- data.frame(x,y) } -- --------------------------------- Douglas M. Hultstrand, MS Senior Hydrometeorologist Metstat, Inc. Windsor, Colorado voice: 720.771.5840 email: dmhultst at metstat.com web: http://www.metstat.com

...i,12],data[i,14],data[i,16],data[i,18],data[i,20],data[i,22]) print(freq_ppt) return_value <- approx(freq_ppt,yrs,xout=data2[i,3]) print(return_value) } Thanks for your help, Doug -- --------------------------------- Douglas M. Hultstrand, MS Senior Hydrometeorologist Metstat, Inc. Windsor, Colorado voice: 970.686.1253 email: dmhultst at metstat.com web: http://www.metstat.com

...end("topleft", legend=leg.txt, col=c("black","yellow","red","orange","blue","black"), pch = c(16,16,16,16,1,1), cex=0.85 ) Thanks, Doug -- --------------------------------- Douglas M. Hultstrand, MS Senior Hydrometeorologist Metstat, Inc. Windsor, Colorado voice: 970.686.1253 email: dmhultst at metstat.com web: http://www.metstat.com --------------------------------- -------------- next part -------------- A non-text attachment was scrubbed... Name: ZR_lewis_pre_20100115_1400.png Type: image/png Size: 31386 bytes Desc: not av...

...i in 1:length(d[,1])) { avg[i] = mean(d[i,1:4], na.rm=TRUE) } # Final Output wanted. st1 st2 st3 st4 avg 1 1 2 5 6 3.50 2 2 5 5 5 4.25 3 3 6 NA 7 5.33 4 4 7 7 8 6.50 -- --------------------------------- Douglas M. Hultstrand, MS Senior Hydrometeorologist Metstat, Inc. Windsor, Colorado voice: 720.771.5840 email: dmhultst at metstat.com web: http://www.metstat.com

...38.95 -95.67 1.000 5 8.553 > d.new lat.add lon.add PPT.add Z.add sd_dif 1 37.08 -95.57 0.508 106 -999.000 2 38.85 -94.73 1.000 5 -999.000 Thanks for all the help, Doug -- --------------------------------- Douglas M. Hultstrand, MS Senior Hydrometeorologist Metstat, Inc. Windsor, Colorado voice: 970.686.1253 email: dmhultst at metstat.com web: http://www.metstat.com

...00824 22,25,0.04258 22,50,0.057 22,100,0.101 23,1,2.10328 23,5,2.02755 23,10,1.93808 23,25,1.78408 23,50,1.67407 23,100,1.568 24,1,3.20842 24,5,3.09228 24,10,2.95452 24,25,2.71661 24,50,2.54607 24,100,2.38108 -- --------------------------------- Douglas M. Hultstrand, MS Senior Hydrometeorologist Metstat, Inc. Windsor, Colorado voice: 970.686.1253 email: dmhultst at metstat.com web: http://www.metstat.com

...quot;)"), xlab="Maximum Average Depth of Percipitation (inches)", scales=list(y=list(log=TRUE)), ylim=c(1,1000000), type='l', auto.key=list(space='right')) dev.off() -- --------------------------------- Douglas M. Hultstrand, MS Senior Hydrometeorologist Metstat, Inc. Windsor, Colorado voice: 970.686.1253 email: dmhultst at metstat.com web: http://www.metstat.com

...t;, side=1, line=2) mtext("Incremental Precipitation (in)", side=2, line=2.5) mtext("Accumulated Precipitation (in)", side=4, line=2) dev.off() Thanks for all your help, Doug --------------------------------- Douglas M. Hultstrand, MS Senior Hydrometeorologist Metstat, Inc. Windsor, Colorado voice: 970.686.1253 email: dmhultst at metstat.com web: http://www.metstat.com --------------------------------- -------------- next part -------------- A non-text attachment was scrubbed... Name: mass_curve_33.355-111.325_1147_Zone1.png Type: image/png Size: 36251 bytes De...

...later be called for plotting the data. Below is a clip of the for loop I am using, any suggestions? k = 1 for (i in 1:length(stats$hour)) { "x(i)" = dataset[k,(3:15)] k = k+1 } Thanks, Doug -- --------------------------------- Douglas M. Hultstrand, MS Senior Hydrometeorologist Metstat, Inc. Windsor, Colorado voice: 970.686.1253 email: dmhultst at metstat.com web: http://www.metstat.com

...file of data (two observations). Any thoughts would be very helpful. # Read in data file data <- read.table("file.txt", skip=4, header=F) temp <- data$V4 elev <- data$V3 Thank you, Doug -- --------------------------------- Douglas M. Hultstrand, MS Senior Hydrometeorologist Metstat, Inc. Windsor, Colorado voice: 720.771.5840 email: dmhultst at metstat.com web: http://www.metstat.com --------------------------------- -------------- next part -------------- An embedded and charset-unspecified text was scrubbed... Name: file.txt URL: <https://stat.ethz.ch/pipermail/r-help/at...

...0,] d.2 <- d.1[d.1$t < 2.0,] How can I make d.1 and d.2 into one step? I have tried several different methods such as example below. d.1 <- d[d$t >= 1.0, && d$t < 2.0,] Thanks, Doug -- --------------------------------- Douglas M. Hultstrand, MS Senior Hydrometeorologist Metstat, Inc. Windsor, Colorado voice: 970.686.1253 email: dmhultst at metstat.com web: http://www.metstat.com

...00010.00 00006.91,0000020.00 00005.90,0000050.00 00005.02,0000100.00 00004.09,0000200.00 00003.55,0000300.00 00002.96,0000500.00 00002.27,0001000.00 00001.64,0002000.00 00000.82,0005000.00 00000.77,0005360.00 -- --------------------------------- Douglas M. Hultstrand, MS Senior Hydrometeorologist Metstat, Inc. Windsor, Colorado voice: 970.686.1253 email: dmhultst at metstat.com web: http://www.metstat.com --------------------------------- -------------- next part -------------- An embedded and charset-unspecified text was scrubbed... Name: da_zone1_15hr_1166.stdsummary URL: <https://stat.ethz.c...

...ate the plot/.png. Example: lm_mod <- lm(data1~data2) intercept = 48.54 slope = 0.4856 I want to display equation on a plot as y=0.4856x+48.54 Any thoughts or suggestions are appreciated. Thanks, Doug -- --------------------------------- Douglas M. Hultstrand, MS Senior Hydrometeorologist Metstat, Inc. Windsor, Colorado voice: 970.686.1253 email: dmhultst at metstat.com web: http://www.metstat.com

...& P[i]<=1) { w[i] = 88 } else { w[i] = 77} } d.f <- data.frame(Z,P,w) d.f Z P w 1 25 0 0 2 23 1 1 3 28 2 2 4 29 3 3 5 30 4 4 6 31 5 5 7 28 0 0 8 29 1 1 9 22 5 77 -- --------------------------------- Douglas M. Hultstrand, MS Senior Hydrometeorologist Metstat, Inc. Windsor, Colorado voice: 970.686.1253 email: dmhultst at metstat.com web: http://www.metstat.com

...? Xi (location) = 35.396 Alpha (scale) = 1.726 Kappa (shape) = 0.397 Year, Return value, CI_u, CL_l 2, 36.0, 5, 37.3, 10, 38.0, 25, 38.5, 50, 38.8, 100, 39.0, Thank you for your help/suggestions, Doug -- --------------------------------- Douglas M. Hultstrand, MS Senior Hydrometeorologist Metstat, Inc. Windsor, Colorado voice: 970.686.1253 email: dmhultst at metstat.com web: http://www.metstat.com

...to square grids r2 = raster(r) res(r2) = min(res(r)) res(r2) r2 <- resample(r, r2, method='bilinear') writeRaster(r2,file="test_grid24.asc", format="ascii") Thanks, Doug -- --------------------------------- Douglas M. Hultstrand, MS Senior Hydrometeorologist Metstat, Inc. Windsor, Colorado voice: 720.771.5840 email: dmhultst at metstat.com web: http://www.metstat.com ---------------------------------

...ded packages are in /tmp/RtmpILmmXw/downloaded_packages Warning message: In install.packages("Cairo") : installation of package 'Cairo' had non-zero exit status > Thanks, Doug -- --------------------------------- Douglas M. Hultstrand, MS Senior Hydrometeorologist Metstat, Inc. Windsor, Colorado voice: 970.686.1253 email: dmhultst at metstat.com web: http://www.metstat.com