similar to: Creating Vectors

Dear All, I'm trying to save a dataframe using write.table command. It works, but when I retrieved, there's an error message as shown below: > write.table(soil.dat,file="C:/soil.rdata") > load("C:/soil.rdata") Error: bad restore file magic number (file may be corrupted)-- no data loaded I can figure out the error message. Any assistance to solve the

Hi I have data like 1, A, 24, The Red House 2, A, 25, King's Home, by the Sea ... I'd like to read this in as three variables. I first tried temp <-read.csv(addresses, sep = "," ) it worked but line 2 was broken after King's Home, and by the Sea as placed in another line. and so i eneded up with more number of rows than in the data. when i tried temp

Dear R-Help, I am trying to find a function that will give me the significance of the correlation of 2 variables (in the same dimension arrays) correcting for serial autocorrelation. How can I view the function cor.test's code? I would like to know a lot more detail about the function than written in the documentation at http://finzi.psych.upenn.edu/R/library/stats/html/cor.test.html to

Dear list, When compiling the R-2.3.0 on FC4 x86_64, I got the following errors: make[3]: Entering directory `/project/scratch3/ligroup/wuming/src/R-2.3.0/src/main' gcc -Wl,--export-dynamic -L/usr/local/lib64 -o R.bin Rmain.o CConverters.o CommandLineArgs.o Rdynload.o Renviron.o RNG.o apply.o arithmetic.o apse.o array.o attrib.o base.o bind.o builtin.o character.o coerce.o colors.o complex.o

Hello, I'm looking for a way to plot different data series on one graph. I have a series of hourly rainfall and quarterly flow measurements (i.e. 4 times an hour) of a catchment. The rainfall should be plotted in bars, the flow as a line. Both on the same X axe (time) but with different Y axes. The problem is the plot() function does not support add=TRUE... Furthermore I'm not sure

Tengo una consulta, alguien tiene algun script en R para hacer graficos de anomalias estandarizadas (precipitacion, temperatura, etc) y que me permita pintar de rojo anomalias positivas y azules las negativas, ​Muchas gracias de antemano​ -- Luis Alfonso LOPEZ ALVAREZ MSc- Meteorologia - Ing. Sistemas PROFESIONAL ESPECIALIZADO (E) METEOROLOGO - IDEAM e-mail: llopez@ideam.gov.co

Hello, I am developing GlobalISel for MIPS. I have a few questions and observations about defining legality of generic instruction and also possible combining of instructions and artifacts in pre/post legalizer combiner or elsewhere (e.g. in some sort of instruction-select patterns). I look at legality as "If generic instruction can be selected into machine instruction, it is legal".

Hi According to vorbis I spec, there are several configurations. There are 2 formats of floor, 3 formats of residue and 3 look up types of VQ codebook, etc. I think to get all kinds of vorbis file with these configurations but what I have got from network only support some of the formats (Ex. floor1, residue 1 & 2, VQ lookup type0 & 1). I downloaded several freeware vorbis encoder

I am new to R and am looking for a book that can help in learning to program in R. I have looked at the R website suggested books but I am still not sure which book best suite my needs. I am interesting in programming, data manipulation not statistics. Any suggestions? Raphael

Dear all, When analyzing data from a climate change experiment using linear mixed-effects models, I recently came across a situation where: - the summary(model) showed a significant difference between the levels of a two-level factor, - while the anova(model) showed no significance for that factor (see below). My question now is: Is the anova.lme() approach correct for that model? And why does

Hi, Can anyone help to write a more elegant version of my code? I am sure this can be put into a loop but I am having trouble creating the objects b1,b2,b3,...,etc. b1 <- rigamma(50,1,1) theta1 <- rgamma(50,0.5,(1/b1)) sim1 <- rpois(50,theta1) b2 <- rigamma(50,1,1) theta2 <- rgamma(50,0.5,(1/b2)) sim2 <- rpois(50,theta2) b3 <- rigamma(50,1,1) theta3 <-

Hi, Could someone tell me if this is the correct model syntax for the following dataset: lme(height~treatment+genotype+treatment*genotype,drought,random=~genotyp e) The dataset has two factors: one fixed - treatment, and one random - genotype. I need to test the effect of both factors to identify their significance. There are multiple (but not equal) replicates at each level of genotype (the

I have been trying to do tukey's test to no avail. This is what I have and the error produced: pen.df = data.frame(blend, treatment, y) source("tukey.1.r") tukey.1(aov.pen, pen.df) Error in tukey.1(aov.pen, pen.df) : the model must be two-way This is a two-way design. Therefore I am confused. Can anyone help? Raphael > pen.df blend treatment y 1 1 A 89 2

Hi @all, I am biting my nails with the following problem: I have a data set of daily rainfall measurements for the last 20 years. What I want to do is calculate the daily cumulative sum of rainfall but only for every year which means that the cumulative sum has to be reset each year. After the calculations I want to plot each year of cumulative rainfall as a separate line in one graph preferably

I've been trying to split a space delimited string with double-quotes in R for some time but without success. An example of a string is as follows: /rainfall snowfall "Channel storage" "Rivulet storage"/ It's important for us because these are column headings that must match the subsequent data. Here is some code I've been trying: str <- 'rainfall

Dear ALL, I am trying to run Tinn R for the first time with R. I have already selected the SDI option. When ever I run the code I get the following error: > source(.trPaths[4]) Error in source(.trPaths[4]) : object '.trPaths' not found Can any one help? Raphael

Hi, For my master thesis I have 24 micro-plots on which I did measurements during 3 months. The measurements were: - Rainfall and runoff events throughout 3monts (runoff being dependant on the rainfall, a coefficient (%) has been made per rainfall event and per 3 months) - Soil texture (3 different textures were differentiated) - Slope (3 classes of slopes) - Stoniness (one time measurement)

hi i try to use arima and holtwinter to predict drought from 1895-2006 but i cannot read whole period of time and i try to do the exponent fitting, but it comes out as the coordinate x-y error i send the source code and data to take a look if anyone can help me, i am really new in R thank u so much jam

Dear All, I am a research student in environment. I have only little programming knowledge. I am currently doing the last project about rainfall impact on ground water quality in an area. It happens that I have to use R to read rainfall data (3 dimension) from ASC file (*.asc), and then write them into one NCDF file (*.nc). I have been working very hard on study R, but I

Fwiw the engine behind geom_raster needs explicit observation-per-row form for input (with no structural normalization), so conversion to points is perfectly proper here, albeit confusing in context. (It's closer to what graphics devices actually use ultimately, but the expansion is laid out very early in ggplot2 because there's no standard for intermediate forms.) Cheers, Mike On Wed,