search for: rain

...R has given me. The code I have used, and some of the data, is as follows (once the database has been imported from SQL, and outliers removed). Z1sub <- Z1[, c(2, 5, 7,11, 12, 13, 15, 16)] colnames(Z1sub) <- c("temp", "hum", "wind", "press", "rain", "s.rad", "mean1", "sd1" ) attach(Z1sub) names(Z1sub) Model1d <- lm(mean1 ~ hum*wind*rain + I(hum^2) + I(wind^2) + I(rain^2) ) summary(Model1d) Call: lm(formula = mean1 ~ hum * wind * rain + I(hum^2) + I(wind^2) + I(rain^2)) Residuals: Min...

...two things from my data frame and was after some advice. For the example below i want to know. 1. How many unique Orders/Families and Genera there are per eco-name 2. How many incidences are there for each Order/Family and Genus there are per eco-region I have 650 econame. I.e for Biak-Numfoor rain forests there are 2 orders, 2 families and two genera. Also, Alismatales are represented once, Asparagales once etc etc. Thanks for any advice. Chris ECO_NAME Order Family Genus Biak-Numfoor rain forests Alismatales Araceae Homalomena Biak-Numfoor rain forests Asparagales...

I am using windows XP with R 2.8.1 I am generating a lattice plot of annual rain patterns using the following function: > xyplot(rain.stats$min+ rain.stats$max + rain.stats$ave ~ rain.stats$month |rain.stats$year, lty = 1, data = rain.stats, type = c("l","l", "l"), col = c("red", "blue", "green"), d...

Hi Gabor and everybody; Thanks Gabor, with the alternative step you've told me the problem is resolved. Comparing the two procedures: Extract from the source 'character' data: > rain$ts[2039:2046] [1] "25/03/2000 22:00:00 UTC" "25/03/2000 23:00:00 UTC" [3] "26/03/2000 00:00:00 UTC" "26/03/2000 01:00:00 UTC" [5] "26/03/2000 02:00:00 UTC" "26/03/2000 03:00:00 UTC" [7] "26/03/2000 04:00:00 UTC" "26/03/2000...

Hi, If I have the follow situation: A dependent variable (i.e. number of insects) that is affected by an independent variable (i.e. rain). The problem is that the measure of rain affect the population in other moment. So there exit a lag between the rain and the number of insects. Exist in R any tool to find what is this lag? Explain better. Suppose that I have a linear relationship between rain and insects. More rain make mor...

...d code followed by sessionInfo(). Thanks in advance for your help. I'm going to be doing a lot of text mining in the near future. I have an excellent O'Reilly book on regex's. What is the best reference for R's special treatment of these animals? Mark myCharVec <- c("[the rain in spain]", "(the rain in spain)") gsub('\\[*.\\]', '', myCharVec) #what I get # [1] "[the rain in spai" "(the rain in spain)" #what I want [1] "" "(the rain in spain)" > sessionInfo() R version 2.10.0 Under developmen...

See ?rpart.control. I get: > golf.rp = rpart(Outlook ~ ., golf, control=rpart.control(minsplit=1)) > golf.rp n= 14 node), split, n, loss, yval, (yprob) * denotes terminal node 1) root 14 9 rain (0.2857143 0.3571429 0.3571429) 2) Temperature< 71.5 6 2 rain (0.1666667 0.6666667 0.1666667) 4) Temperature< 64.5 1 0 overcast (1.0000000 0.0000000 0.0000000) * 5) Temperature>=64.5 5 1 rain (0.0000000 0.8000000 0.2000000) 10) Humidity>=75 3 0 rain (0.0000000 1...

Oats[Oats$Variety %in% c("Victory", "Golden Rain"),] or subset(Oats, Variety %in% c("Victory", "Golden Rain")) Stefano -----Messaggio originale----- Da: r-help-bounces at stat.math.ethz.ch [mailto:r-help-bounces at stat.math.ethz.ch]Per conto di elyakhlifi mustapha Inviato: luned? 23 aprile 2007 9.56 A: R-help at stat.m...

Dear anyone who knows more about R than me (so everyone). I have been bashing my head on the keyboard all day trying to do something with my table. I have some data, like so: yyyy-mm Rainfall(mm) 1 1977-02 17.4 2 1977-03 34.0 3 1977-04 26.2 4 1977-05 42.6 5 1977-06 58.6 6 1977-07 23.2 7 1977-08 26.8 8 1977-09 48.4 9 1977-10 47.0 10 1977-11 37.2 11 1977-12 15.0 12 1978-01 2.6 13...

Hi r-users,   I have this data below.  I would like to obtain the weekly rainfall sum.  That is I would like to find sum for day 1 to day 7, day 8 - day15, and so on.    year month day rain 1  1922     1   1  0.0 2  1922     1   2  0.0 3  1922     1   3  0.0 4  1922     1   4  0.0 5  1922     1   5  0.0 6  1922     1   6  0.0 7  1922     1   7  0.0 8  1922     1   8  6.6 9 ...

hello, I'd like know how to do to extract data from a frame for example how can I do to extract only the data where variety=victory or variety=golden rain thanks. > Oats Block Variety nitro yield 1 I Victory 0.0 111 2 I Victory 0.2 130 3 I Victory 0.4 157 4 I Victory 0.6 174 5 I Golden Rain 0.0 117 6 I Golden Rain 0.2 114 7 I Golden Rain 0.4 161 8 I Golden...

Hey Guys, I am quite new to R, am trying to plot a bar plot for my monthly rainfall data from 1920 to 1990. For every year I have the data from January to December. I have created this code: data<-read.table("C://R//Data//Rain.txt", header=TRUE) data<-data[,2:13] Rdata<-data[1,] for(i in 2:79){Rdata<-c(Rdata,data[i,])} times<-ts(Rdata,start=1920,frequ...

Hello, I love the work that is being done with compiz! I have made a patch for water.c that makes the time between rain-drops adjustable using gconf. This way you can leave the plug-in running and just have pleasant drops every now and then instead of the rainstorm that is hard-coded in now (which is still fun now and then!). The gconf value is in ms (1/1000 of a second). Suggestions for raindrop times: 10000 (ever...

Hello, I'm trying currently using barplot to summarize precipitation data. So far I've compiled total annual snow and rain accumulation in a table (attached). I've been successful at plotting it using the following code: fig=barplot(t(Annual_Precip_table), horiz=TRUE, space=0, col=c("grey70", "white"), axes=FALSE) . . The result is a stacked barplot with total annual rain represented by grey b...

...nts, ie, user count etc.. I swear I seen one before but I can't remember where. -=-=-=-=-=-=-=Ottawa Valley Weather-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-= Weather @ 2:20pm - Temp: 34.9?C - Humidity 29% - Wind: E @ 13 km/h Baro: 996kPa Steady - Vis: 8km - Sky: --- - Weather: Haze Hourly Rain: 0.00mm - Daily Rain: 0.00mm - Total Rain(May 28th): 18.00mm =-=-=-= Website: http://www.WeatherServer.net =-=-=-=-=-=-=-= Our Alert Lists: MTO-PEEL, MTO-TORONTO, MTO-YORK, NHC, OntarioDiscussion, SPC, USThunderStormWarnings, USTornadoWarnings, WxDispatch Signup at http://listserver.weatherserver....

...es in R, so some advice would be very appreciate. Here is a simple example similar to my real problem (R 2.6.0 for windows): Lets supose I have annual tree production measured in 9 trees during 3 years; the 9 trees are located in 3 different mountains (sites), and each tree receive different annual rainfall (different locations). I would like to known the parameters that explain the variability in production. The data would be something like: set.seed(111) mydat <- data.frame(tree= factor(rep(1:9,3)), year= gl(3,9,lab=2001:2003, orde=T), site= gl(3,3,27,lab=c("A","B","...

Oh sorry. An example say lots more than words. The data below, when submitted to twoord.plot return the mentioned error. Rain are bars and Salt lines, the bars appear and the error occurs with the salt data. Point Rain Salt Fev/03 365.6 13 Mar/03 235 18 Abr/03 115.1 18 Mai/03 47.4 18.75 Jun/03 112 15 Jul/03 156.8 17 Ago/03 66.1 15 Set/03 149.8 14 Out/03 167,1 11.5 Nov/03 269.3 17.5 Dez/03 283.7 Jan/04 415 1,3 Fev/04 322...

...nb) model1AIC.glm.nb<-stepAIC(model1.glm.nb) When I run it, it doesn't give me any significance. So I run a second model, with my autocorrelation table, and remove one more variable that does not autocorrelate with other (anything below .07) model2.glm.nb2 <- update(model1.glm.nb, ~ . - rain, data=climate.dat) summary(model2.glm.nb2) Now, when I run a 3rd. model to remove rain, it shows after summary that rain has been removed. Once I do a forth model rain is again in it. How can I remove more than one variable in one go? I was thinking doing something like this: model2.glm.nb2 <...

...s sync. Also, like for mail spools, large files that simply get appended to get the whole file recopied. Does anyone know of something that syncs an ext2/3 fs to another at the block level which result in less data transfer? -- --------------------------------------------------------------- Paul Raines email: raines@nmr.mgh.harvard.edu MGH/MIT/HMS Athinoula A. Martinos Center for Biomedical Imaging 149 (2301) 13th Street Charlestown, MA 02129 USA

...ook and extracted the variance and correlation components: >library(nlme) >data(Oats) >fm4OatsB <- lme(yield~nitro, data=Oats, random=list(Block=pdCompSymm(~Variety-1))) >VarCorr(fm4OatsB) Block = pdCompSymm(Variety - 1) Variance StdDev Corr VarietyGolden Rain 331.5271 18.20788 VarietyMarvellous 331.5271 18.20788 0.635 VarietyVictory 331.5271 18.20788 0.635 0.635 Residual 165.5585 12.86695 This is a compound symmetry variance-covariance structure. I then tried to find out the standard deviation and correlation ma...