similar to: qplot, stacked area, own colourscheme

Dear List, yes, me again trying to work with qplot ;-) I would like to make several single plots within a loop, like this (simplified and so on...): trials <- c("A","B","C") mycolours <- ("wheat","darkolivegreen","lightgreen",

Dear list users, I have trouble with covariance analysis. I measured nitrate concentrations in the soil (NO3) and the percentage of legumes (LEG, continuous), affected by 2 different CO2 concentrations (CO2, discrete). I suspect that CO2 has an effect on LEG and NO3, but also that LEG has an effect on NO3, so this is the formula I wrote to test this: NO3 ~ CO2 + LEG + CO2:LEG Will LEG be

Hi everybody, Well, I really hoped that until now i know how to do a selection ..... but obviously i don't. So risking your laughing at me - here you go. I have a table with lots of vegetation data, and one column specifies in which region the data is. There are only 2 regions, region 1 and region 3. I want to select all data that belong to region 1 and after do some xtabs on it. If i do

Dear all, I'm getting an error in one of the stock examples in the 'mvpart' package. I tried: require(mvpart) data(spider) fit3 <- rpart(gdist(spider[,1:12],meth="bray",full=TRUE,sq=TRUE)~water+twigs+reft+herbs+moss+sand,spider,method="dist") #directly from ?rpart summary(fit3) ...which returned the following: Error in apply(formatg(yval, digits - 3), 1,

Dear all I'm not sure if I did the right analysis for my specific split splot design. We are studying biomass increase with different CO2 concentrations with four different functional plant groups (e.g. grasses, herbs, broad-leafed trees and conifers). Of each functional plant group we have four species. The design is orthogonal. The design is: Blocks: 2 (climate chambers, called

Dear List, I have got error message when I run the R-code. Can anyone has a suggestion? v.code <- df.bm7[,c(10:31)]; v.code[1:3,] names(v.code) CM = v.code # variable binomial code sim.sp <- function(data,CM,n,N) { C <- matrix(rep(NA,N),ncol=1) for(i in 1:N) { j <- n xx <- which(colSums(CM[j,])==1) V <- names(xx) V <- paste(V,

Dear all, I was wondering if there is a simple way to avoid printing the multiple cross-validation automatic output to the console of recursive partitionning functions like rpart or mvpart. For example... > data(spider) > mvpart(data.matrix(spider[,1:12])~herbs+reft+moss+sand+twigs+water,spider,xv="1se",xvmult=100) *X-Val rep : 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15

Dear R-group, I am trying desperately to get this Tukey test working. Its my first time here so I hope my question is not too stupid, but I couldn't find anything helpful in the help or in the forum. I am analysing a dataset of treated grasses. I want to find out, if the grasses from different Provenances react differently. In the aov test I found a significance for the combination Treatment

Dear List, I'm using gam in a multiple imputation framework -- specifying the knot locations, and saving the results of multiple models, each of which is fit with slightly different data (because some of it is predicted when missing). In MI, coefficients from multiple models are averaged, as are variance-covariance matrices. VCV's get an additional correction to account for how

Hello all, turns out i'm having a bad R week. I am at my wits end with a function that I am trying to write. When I run the lines of code outside of a function, I get the desired output. When I wrap the lines of code into a function it doesn't work as expected. Not sure what is going on here. I suspected that the syntax of the if statement with the for loop was the culprit, but when I

Hello. I am sorry for posting a (seemingly) simple question, but I have just spent 2 hours trying to find the answer, without success. I want to make a histogram with conditioning on a factor, using Trellis graphics. However, I do not want any colours (only black and white) either in the histograms or in the strip. There must be some simple argument but I can’t find it. Here is my code so

Dear all, Does any one knows why should I get the following error message, when trying to do a simple data.frame?? DataF<-data.frame(Subject,BiomR,Spp,Capas,Litter,Herbs,LitterD,MaxCanH,DDifS p,DSSp,Slope, CanDens,NearestSp) Erro em data.frame(Subject, BiomR, Spp, Capas, Litter, Herbs, LitterD, : arguments imply differing number of rows: 202, 0 The data I am using

I have a simple barplot that looks fine in the R graphics device window. However when I export it to png I am losing the x-label. It must be an obvious problem but I cannot see it. Trying to resize the plot does not seem to help. Code is below. Any help gratefully received. ###### Start Code############ Groups <-c(21.8,45, 43, 17.2, 8.3, 18) names(Groups) <- c("Exeter",

Someone was asking how to do a 16 category piechart in OpenOffice Calc and it appears that it can not be done (which we, probably, should be happy about) but I thought that I'd try it in ggplot2. It works but I then thought I'd like to make the colours more distinctive but fro some reason I don't seem to be able to use manually assigned colours. Can anyone suggest where I'm

Dear all, Suppose I have a linear mixed-effects model (from the package nlme) with nested random effects (see below); how would I present the results from the random effects part in a publication? Specifically, I?d like to know: (1) What is the total variance of the random effects at each level? (2) How can I test the significance of the variance components? (3) Is there something like an

Dear List, I'm getting an error in mgcv, and I can't figure out where it comes from. The setup is the following: I've got a fitted GAM object called "MI", and a vector of "prediction data" (with default values for predictors). I feed this into predict.gam(object, newdata = whatever) via the following function: makepred = function(varstochange,val){ for

Bill Venables expected howls of protest about deprecation of underscore for assign. Here's my one-and-only howl on this point: (1) I don't like <- as an assignment operator, because it creates an ambiguity with comparisons: I can write "x+1", "x-1", "x<1" but not "x<-1". However, I guess it's a bit much to expect the entire R/S

Thank you for taking time to do this, Andrew. But, of course, it will be too late for me. I’ve just returned from exile, where I went on a spirit quest of sorts. Except that, on this quest, I was obliged to keep my distance until I had found a way to export and import all users, groups and group membership from my old samba domain to my new one. I updated schema to support Kerio Connect using

Dear R community, I am analyzing plant communities with the function mvpart, using a dissimilarit matrix as input. The matrix is calculated with the funtion gdist. fit <- mvpart(gdist (ba12[,18:29], meth="maximum", full=TRUE, sq=F) ~ beers + slope_dem + elev_dem+ plc_dem + pr_curv+ +curv+max_depth+doc_rocks+ abandon+land_use+ca_old, data=ba12, xv="p") This

In VIM one can easily change colours with ":set backgorund=dark". This doesn't actually change the background, but rather uses a colour scheme that is designed for a dark background. Is there any quick command like this for bash? I don't want to edit the whole config file if there is a quick way to get a better colour scheme. Thanks. -- Dotan Cohen http://gibberish.co.il