similar to: write values of points beside points in plot

Hi there, I trying to solve this problem for the whole day not going anywhere, so I really hope maybe somebody can help me in this community... I've got an object coefficient2 which I want to plot in differerent ways, with colors and labels added to the points, but somehow there seems to be a problem if a value is NA within the independent variable, resulting in false labels and false colors

Hi, how do I create a string of the comma-separated content of a vector? I've got the vector i with several numeric values as content: >str(i) num 99 and want to create a SQL statement to look like the following where the part '(2, 4, 6, 7)' should be the content of the vector i: select * from [biomass_data$] where site_no in (2, 4, 6, 7) Here my approach (which doesn't

Hi there, I've got a database or rather spreadsheet with several columns and rows. For one column named sites I want to loop through all possible values and retrieve all data out of the database where site = x and write it into an object named 'sitex_data'. Somehow I'm really missing something as I'm not able to create these sitex_data objects with the database values,

Hi there, sorry for troubling everybody once again, I've got a problem rewriting Sarkar's function for rewriting the tick locations in a logaritmic way (s. http://lmdvr.r-forge.r-project.org/code/Chapter08.R): His example works for log 2 but I need log e (natural logarithm). My problem is that if I replace 2 with "e" (using paste()), I get the error message that the location

Hi there, I've got a question which is really trivial for sure but still I have to ask as I'm not making any progress solving it by myself (please be patient with an undergraduate student): I've got a linear model (lm and lmer fitted with method="ML"). Now I want to compare the coefficients (slope, intercept, not the random effects) of both models with a given value (e.g.

hi all, after fitting a multiple linear regression model <- lm(y ~ a + b+ c+d) i wanted to plot diagnostics plot(model) but get the error message Error in object$coefficients : $ operator is invalid for atomic vectors. which does not make a lot of sense, since there is no "$" - i am working with an attached dataset. can anyone help, please?? thanks a lot, kat

Hello, I have a dataframe containing various time series (not time series objects though!)with hourly time steps. I?d like to perform ccf for I need to know the correlation factors for different lags. Here is an example: x<-as.POSIXct(c("2008-12-25 16:00:00", "2008-12-25 17:00:00", "2008-12-25 18:00:00", "2008-12-25 19:00:00", "2008-12-25

hi, i am trying to work with the survey package in order to apply survey design weights. the data set i am using - ess - contains missing values. my question: when using svytable(~variable1+variable2, design=my.svydesign.object, na.rm=T) the resulting crosstable contains all missings although i defined the na's as such. any ideas??? thanks a lot, katharina

Hi, how can I use rm() on objects named like: paste("site",i,"_data",sep="") while looping through i? I tried rm(paste("site",i,"_data",sep="")) but I get the error that rm() must contain names or text strings which is confusing me as I thought paste() would create something like that...? Thanks, Katharina -- Time flies

Dear all, I would like to produce a map with information about the patenting activity in German districts, by coloring districts with different degrees of patenting activity in different colors. I work with the packages maptools, maps and spdep. The map data is read from an external .shp file (+ the corresponding .shx and .dbf files). Plotting a map with the IDs or the patenting indicator itself

Hello, I have a dataframe assigning various scores on around 20 variables to a list of countries. The scores are rated on a scale of (D, C, B, A) and there are also some not rated ones (NR) and others are left blank (NA). I now wanted to transfer the scores into numeric values (such as NR=0, D=25, C=50, B=27, A=100), while also leaving the NA in the dataframe. Can I transform the whole

Hi all! I have lots of functions called in the following pattern 'NameOfFunctionNumber' where the name always stays the same and the number varies from 1 to 98. Another function which I run in advance returns the number of the function which has to be called next. Now I want to combine 'NameOfFunction' with the 'Number' returned so that i can call the desired function. I

Hi, I've got the following problem which I cannot think of a solution right now: if got a lattice xyplot in black and white and a grouping variable with many (more than 8 values) and I plot it as regression lines (type="r"), just like this one (not reproducable but that's I guess not the point here): xyplot(log(AGWB) ~ log(BM_roots), data=sub_agwb_data, groups=species,

hi, is there a way to make r open the help file in a separate window? i'm quite new to linux(kubuntu) and used to work with r under windows... thanks, katharina

Hi, there´s something that really bothers me about R and after hours and hours of internet research, I´m still stuck with the same problem: I installed R and it is in Spanish, as the system on my work. I would really like it in English and there seems to be no way to change that!?! I´m doing an abroad internship, so my Spanish is far from perfect and my tutorial is in English. Is there any

Hello, I'm trying to use na.spline (package zoo) to fill some missing data in a time series. this works fine, however, if I apply the 'maxgap' argument, I always get the error: <------ Error in na.spline.vec(x., coredata(object.), xout = xout., ...) : attempt to apply non-function ------> I couldn't find a similar error for this case in the mailing lists and zoo vignette,

note: whole function is below- I am sure I am doing something silly. when I use it like USGS(input="precipitation") it is choking on the precip.1 <- subset(DF, precipitation!="NA") b <- ddply(precip.1$precipitation, .(precip.1$gauge_name), cumsum) DF.precip <- precip.1 DF.precip$precipitation <- b$.data part, but runs fine outside of the function: days=7

Hi there, I'm a bit confused concerning the axis tck setting in the lattice package as the ticks on left sided axis aren't drawn at all with the following setting: dados <- data.frame(varsep = factor(rep(1:2,10)), i = runif(20)) library(lattice) my.theme <- list( axis.components = list(left = list(tck = 1, pad1 = 1, pad2 = 2), top = list(tck = 0, pad1

I am working with data from the USGS with data every 30 minutes from 4/27/2011 to 8/25/2011. I am having trouble with setting the frequency. My R script is below: > shavers=read.csv("shavers.csv") > names(shavers) [1] "agency_cd" "site_no" "datetime" "tz_cd" "Temp" [6] "X04_00010_cd"

Dear R-helpers, I try to perform glm's with negative binomial distributed data. So I use the MASS library and the commands: model_1 = glm.nb(response ~ y1 + y2 + ...+ yi, data = data.frame) and predict(model_1, newdata = data.frame) So far, I think everything should be ok. But when I want to perform a glm with a subset of the data, I run into an error message as soon as I want to predict