similar to: apply a function down each column

Hi! Does anybody know a string function that would calculate how many characters two strings share? I.e. ("Hello World","Hello Peter") would be 7. Thanks. Laetitia

Hi, I know that according to plot theory, one should not have two y-axis. It is just that in my field, we are often depending on plots where water depth is plotted in one y-axis and other parameters such as salinity on the other y-axis to see where the halocline (the depth where the two layers of water are not mixed) is located. I have many other examples. To publish, it is required of me to

Hi lattice users, I would like to add 5-year moving averages to my double y-plot. I have three factors needs to be plotted with moving averages in the same plot. One of these reads off y-axis 1 and two from y-axis 2. I have tried to use the rollmean function from the zoo-packages, but I fail in insering this into lattice (I am not an experienced lattice user). I want to keep the data points in

Hi, I've had a little problem for several weeks now. It is annoying and therefore I will ask for help now: When I write a script with several iterations, I make it write out a text file to save the data during the run. For example I write: if (i %% 25) write.table(output,"temporary_output.txt") Later on, when I read in this output and want to calculate things, R complains that

Hi, I am new to R. How can I get labels/different colours on the lines in my xy-plots. I plot "het" against "temp" conditional on "year" and "station". library(lattice) MyLines <- function(xi, yi, ...){ #Draws line in the panels while avoiding spaghetti-plots I <- order(xi) panel.lines(xi[I], yi[I], col = 1)} xyplot(het ~ temp | fstation,

Hi, I have confirmed temporal correlation problems in my data. Is there a possibility to use corCompSymm for a gamm()? I am an R-beginner. I have very short time series. There are three years and within each year, there are 10 weeks. he 10 weeks are the same every year and have not unique values, I seem not to be able to use AR-1 (I assume that I have too little data for autoregression

Hi, I am attempting to correct my models p-values due to temporal autocorrelations. It is not possible to model the correlation, so I have to make do with the p-value correction. I am feeling a bit thick here....I cannot get the autocorrelation values. What is the argument? My aim is to multiply the dependent variable autocorrelation with the independent variable autocorrelation and then

Hi, Dummy data script and scripts are attached below. I would like to change the plot to look like this:

Hi, I have used the mgcv library to generate a simple additive model. I want to know how to plot the function on the raw data with confidence intervals whan I have TWO variables in the model. I get it to work with one variable but not with two. I am on the limit for what I understand in R, so be gentle. I have read the help file on predict.gam, but did not get any help out of it. #My model:

Hi R people! I am looking for some suggestions writing an if-else function. The idea is to characterize different plots containing counts of variables (here parasites). If a plot has a count equal or higher than 4 for any parasite the function should return a 1 else a 0. Later I can loop the function over all plots. Here I have a little subset of my data: VariablePAR Plot1

Hi, Can anybody show me how to extract all columns in my dataset that are polymorphic? Or phrased in another way I would like to delete all columns that have no more than one letter in it (that are monomorphic). Thank you. Laetitia -------------- next part -------------- An embedded and charset-unspecified text was scrubbed... Name: seqCol.txt URL:

Hi, I want to plot means with standard deviations of Total Nitrogen (TN) across 4 stations (S1-S4) and over 3 years (2007-2009). I want this plot in one panel. I do not want medians (bwplot, boxplot). I have tried a few different packages and it seems that ggplots with plotmeans was the fastest (I am not extremely skilled in writing my own scripts). Unfortunately, there is no grouping

Hi! I need some help to finish my script. I have two tables that I combine randomly to produce a third table. This I do for hundreds of iterations. In the output file I get all the simulated tables after each other. It looks like this (in this case 3 iterations): output file: [[1]] [,1] [,2] [,3] [1,] "GM030005" "WI920024" "CCCC" [2,]

Dear Wu Gong and Peter Ehlers, thank you very much for your help debugging my script. Now I have a general following up question: Is there a straightforward way to rearrange the following dataset so that all first letters of each column will be combined in one column, all the second letters in a second column, all the third ones in a third column and so on, resulting in 7 columns, i.e. for

Hi, I have a data set with two continous variables that I want to plot MEANS (I am not intrerested in median values) on a double-y graph. I also have 2 factors. I want the factor combinations plotted in different panes. Dummy dataset: mydata <- data.frame(factor1 = factor(rep(LETTERS[1:3], each = 40)), factor2 = factor(rep(c(1:4), each = 10)), y1 =

Dear R users, I am trying to fit a glm with a distribution free family, link = log and variance = constant*mu. I guess I have to use the quasi family but the choices of variance are restricted to constant or mu or mu^2..., I don't know the way to choose the variance that I need, i.e. constant*mu. If you have any ideas or advice, please tell me. Thanks, Laetitia Mestdagh Laetitia Mestdagh

Help! Being a bit of a novice, please bear with me if this is a stupid question! I am trying to fit a saturated model to some count data that I have: model<-glm(COUNT~SP*LOC*COL*TIME*TREAT,poisson) but R keeps on crashing and coming up with (occasionally before crashing) an error that states: Algorithm did not converge in: (if(is.empty.model(mt)) glm.fit.null else glm.fit)(x = X, y = Y,

Hello everyone, How can I specify in tcltk file browser the initial directory to "MyComputer" in Windows where Drives and Partition are accessible? And just a little question if anyone knows, is there a way to use the function choose.files under windows to select a directory? Thanks a lot for your help. Laetitia.

Hi, I have R installed under a Mandrake linux system and I don't have shell utilities any more under my R console such as completion when writing a file path, back and forth in the history, bindkeys... Moreover when I quit R by saving, no .Rhistory file is created while the .Rdata is. I don't get how this work, I thought it was based on the user unix shell but it does not seem. Any

Hello, Is there a simple and fast function that returns a vector of the number of replications for each object of a vector ? For example : I have a vector of IDs : ids <- c( "ID1", "ID2", "ID2", "ID3", "ID3","ID3", "ID5") I want the function returns the following vector where each term is the number of replicates for the