similar to: How can i add a color bar with base package

hi, i have a vector like : x<-c(1,15,30,45,60,90,115) i know that step by step i have always more than 10 min(diff(x)) >=11 i want to add for each value a sequence of value:value+9 result should be : 1 2 3 4 5 6 7 8 9 10 15 16 17 18 19 20 21 22 23 24 30 31 (...) 39 45 46 (...) 54 60 61 etc.. how can i do this without a loop (i'm sure there is a "elegant" way like

Hi all I have the next data frame year STOD SLAGNR TAL TALT TALVEKT 1 2002 2120006 57 1 NA 1 2 1997 97030032 57 NA NA NA 3 1997 97030071 27 1 NA NA 4 1997 97030005 57 1 NA NA 5 1997 97020127 60 NA 1 NA 6 2001 1160025 27 1 NA 1 7 1998 98020069 60 1 NA NA 8

Hello, Using a fixed effects linear model (with lm), I can get exact p-values out of the AVOVA table, even if they are very small, eg. 1.0e-200. Using lme (linear mixed effects) from the nlme library, it appears that there is rounding of the p-values to zero, if the p-value is less than about 1.0e-16. Is there a way we can obtain the exact p-values from lme without rounding? used commands:

Dear R users, I'm coming back to R after while. I have a data frame with 200 columns, each column has a name. How to extract all columns to a new dataset, but the specified (by names) ones? I was playing with that for a little bit using the vector syntax but got several syntax errors. Thanks, Rado

What does this warning mean precisely? Is there any reason to care about it? Can I Avoid it by another way of programming? Thank you in advance.

Dear Sir/Madam, I am doing a project related to R. However, it is always difficult find some R functions. The R user guide seems not complete. Is there any free document about all R functions ? Who knows ? please help me. My email is myzhaogong@yahoo.com Thank you very much --------------------------------- [[alternative HTML version deleted]]

Un texte encapsul? et encod? dans un jeu de caract?res inconnu a ?t? nettoy?... Nom : non disponible Url : https://stat.ethz.ch/pipermail/r-help/attachments/20060609/ab46b028/attachment.pl

Un texte encapsul? et encod? dans un jeu de caract?res inconnu a ?t? nettoy?... Nom : non disponible Url : https://stat.ethz.ch/pipermail/r-help/attachments/20060609/84108484/attachment.pl

Hi, I am fairly new to GAM and started using package mgcv. I like the fact that optimal smoothing is automatically used (i.e. df are not determined a priori but calculated by the gam procedure). But the mgcv manual warns that p-level for the smooth can be underestimated when df are estimated by the model. Most of the time my p-levels are so small that even doubling them would not result

hi, i'm working with mgcv packages and specially gam. My exemple is: >test<-gam(B~s(pred1)+s(pred2)) >plot(test,pages=1) when ploting test, you can view pred1 vs s(pred1, edf[1] ) & pred2 vs s(pred2, edf[2] ) I would like to know if there is a way to access to those terms (s(pred1) & s(pred2)). Does someone know how? the purpose is to access to equation of smooths terms

I have just started using R for my PhD. I am importing my data from Excel via notepad into Word. Unfortunately, my data has many missing values. I have put '.' and this allowed me to import the data into R. However, I now want to interpolate these missing values. Please can someone give me some pointers as to the method/code I could use? Thankyou, Lillian.

hi, i'm using gam() function from package mgcv. if G is my gam object, then >SG=summary(G) Formula: y ~ +s(x0, k = 5) + s(x1) + s(x2, k = 3) Parametric coefficients: Estimate std. err. t ratio Pr(>|t|) (Intercept) 3.462e+07 1.965e+05 176.2 < 2.22e-16 Approximate significance of smooth terms: edf chi.sq p-value s(x0)

Dear R users, Could somebody please help me to find a way of comparing nonlinear, non-nested models in R, where the number of parameters is not necessarily different? Here is a sample (growth rates, y, as a function of internal substrate concentration, x): x <- c(0.52, 1.21, 1.45, 1.64, 1.89, 2.14, 2.47, 3.20, 4.47, 5.31, 6.48) y <- c(0.00, 0.35, 0.41, 0.49, 0.58, 0.61, 0.71, 0.83, 0.98,

Dear All: I would any appreciate any help with this plot I am struggling with. I have 4 estimates (95% CIs) I want to plot. I want the CI lines to be horizontal on each plotted point. I was trying to tweak some old codes (was for a vertical CI lines) into horizontal but not much dice. Many thanks in advance for your help. YA My working codes: x=c(1,1,1,1.1,1.1,1.1,2,2,2,2.1,2.1,2.1)

When concatenating mbox files like described here https://xaizek.github.io/2013-03-30/merge-mbox-mailboxes/. You will end up with an 'unsorted' mbox file. Is this going to be a problem esspecially when they are large >2GB's and new emails will be written to it? The email client nicely sorts the message from folder A "foldera 5 last" as last, but of course the mbox is

Hi there, I think I have found a small problem in the "/my/path/R-2.0.1/src/modules/X11/MakeFile" generation. During the configure step, I have specified a specific JPEG library by exporting CPPFLAGS variable. All compilation works well for individual files in the src/modules/X11/ directory, but when linking, the -ljpeg option doesn't work. I obtain the following message (in french

Hello, I'm trying to develop a box plot of time series data to look at the range in the data values over the entire period of record. My data initially starts out as a list of hourly data, and then I've been using this code to make this data into the final ts array. # Read in the station list stn.list <- read.csv("/home/kbennett/fews/stnlist3", as.is=T, header=F) # Read in

Hello, As I am relatively new to the R environment this question may be either a) Really simple to answer b) Or I am overlooking something relatively simple. I am trying to add a VarIdent structure to my gam model which is fitting smoothing functions to the time variables year and month for a particular species. When I try to add the varIdent weights to variable Month I get this error returned.

I would like to use lattice graphics to plot multiple functions (or groups or subpopulations) on the same plot region, using different line types "lty" or colors "col" to distinguish the functions (or groups). In traditional graphics, this seems straightforward: First plot all the data using 'type="n"', and subsequently execute a series of "points"

Dear list, I would like to plot points with two types of labels, one at the data point (the name of the point) and another offset a bit with another factor which is either of the two greek characters alpha or beta. I have tried to get the routine to plot a greek character with expression() or with substitute() and have not yet had any success. The following only plots the word in english in