search for: mmol

...code which look like this ------------------------------------------------------------- gof <- benthic.flux(ID="Gulf of Finland", meas.conc=conc, bw.conc=bw.conc, time=times, substance=expression(DIC~(mmol~m^{-3})) ) ------------------------------------------------------------- I get the output by Sweave in my pdf file, like this: --------------------------------------------------------------- > gof <- benthic.flux(ID = "Gulf of Finland", meas.conc = conc, + bw...

...quot;black",0.4) ,width=1) + geom_errorbarh( aes(xmax = NC60.DATA$Ce+NC60.DATA$sdCe ,xmin=NC60.DATA$Ce-NC60.DATA$sdCe) ,colour=alpha("black",0.4) ,width=1) + geom_line(data=NC60.Res4.curve ,aes(x,y) ,size=1 ,colour=alpha("black",0.5)) + xlab(C[e]~(mmol/m^3)) + ylab(q[e]~(mmol/m^3)) print(NC60.iso.graph) dev.off() Regards/Cordialement ------------- Benoit Boulinguiez Ph.D student Ecole de Chimie de Rennes (ENSCR) Bureau 1.20 Equipe CIP UMR CNRS 6226 "Sciences Chimiques de Rennes" Avenue du Général Leclerc CS 50837 35708 Rennes C...

...the real data, let me just post a simulation that shows a contourplot reasonably similar to what I am seeing. TC.ran <- exp(rnorm(400,1.5,.3)) HDL.ran <- exp(rnorm(400,.4,.3) ) f1<-kde2d(HDL.ran,TC.ran,n=25,lims=c(0,4,2,10)) contour(f1$x,f1$y,f1$z,ylim=c(0,8),xlim=c(0,3),ylab="TC mmol/L", xlab="HDL mmol/L") lines(f1$x,5*f1$x) # iso-ratio lines lines(f1$x,4*f1$x) lines(f1$x,3*f1$x) Two questions: Is there a 2d density estimation function that has provision for probability weights (or inverse sampling probabilities)? I seem to remember a discussion on...

Hi all, I'd like to use superscript or subscript in the axis label of a graph. Is that possible in R? For instance in xlab="Ce (mmol/m3)", I'd need the "3" of cubic meter in superscript. Regards/Cordialement ------------- Benoit Boulinguiez Ph.D Ecole de Chimie de Rennes (ENSCR) Bureau 1.20 Equipe CIP UMR CNRS 6226 "Sciences Chimiques de Rennes" Campus de Beaulieu, 263 Avenue du Général Leclerc 3...

...le line? Cheers, David y.bar.new = c(30, 15) se.ybar.new = c(2,3) error.bars(y.bar.new, se.ybar.new, c("Control (n=18)", "CKD (n=18)")) error.bars<-function(yv,z,nn){ xv<-barplot(yv,ylim=c(0,(max(yv)+max(z))),names=nn,ylab="Total five hour potassium excretion (mmol)", sub= "Figure 1B: Hour 1-5 potassium excretion") g=(max(xv)-min(xv))/50 for (i in 1:length(xv)) { lines(c(xv[i],xv[i]),c(yv[i]+z[i],yv[i]-z[i])) lines(c(xv[i]-g,xv[i]+g),c(yv[i]+z[i], yv[i]+z[i])) lines(c(xv[i]-g,xv[i]+g),c(yv[i]-z[i], yv[i]-z[i])) }} > sessionInfo() R version...

...} params1 = c(um = 1, #maximal phytoplankton growth rate, d-1 Kn = 0.2, #nitrogen half-saturation constant uM/L Kw = 0.04, #light attenuation due to water, unit: m^-1 Kc = 0.03, #light attenuation by phytoplankton, #unit: m^2 (mmol N)^-1 alpha = 0.025, Im = 0.6, #zooplankton maximal ingestion rate, unit: /d Kp = 1, #zooplankton half-saturation for ingestion, unit: mmol N/m^3 eps = .3, #Growth efficiency of zooplankton delz = 5,...

...ose) ncol <- length(names(glucose)) xrange <- range(1, ncol) yrange <- range(0, max(!is.na(glucose[1:ncol]))) nrecs <- nrow(glucose) colors <- rainbow(nrecs) linetype <- c(1:nrecs) plot(xrange, yrange, type="n", xlab="Measurement moment",  ylab="Glucose (mmol/liter)") for (i in 1: nrecs) {   daily_values <- glucose[i,]   lines(daily_values, type="b", lwd=1.5, lty=linetype[i], col=colors[i]) } So I want to loop over the file and add the plot lines one by one (day by day), but it seems that something's not yet right because nothing...

Dear list, I get the following message while computing truehist in R 1.4.1 on Redhat Linux 7.1: > truehist(lsk$Pox, nbins = "FD" , prob = TRUE, xlab = "Pox [mmol/kg]") Error in switch(casefold(nbins), scott = nclass.scott(data), "freedman-diaconis" = , : couldn't find function "nclass.fd" Maybe the "nclass.fd" should be changed into "nclass.FD"? But I could not find it in the directories. Regar...

...do not know how to pass on these parameters as a vector this productes and error rGrow = c(1.0, 1.3, 1.7) # /day, growth rate of prey rMort = 0.2 , # /day, mortality rate of predator assEff = 0.5, # -, assimilation efficiency K = 10) # mmol/m3, carrying capacity yini <- c(Prey = c(1, 2,3), Predator = 2) times <- seq(0, 200, by = 1) out <- ode(yini, times, LVmod, pars) summary(out) Thanks Frank

...("black",0.4)    ,width=1) + geom_errorbarh(   aes(xmax = NC60.DATA$Ce+NC60.DATA$sdCe    ,xmin=NC60.DATA$Ce-NC60.DATA$sdCe)    ,colour=alpha("black",0.4)    ,width=1) + geom_line(data=NC60.Res4.curve   ,aes(x,y)   ,size=1   ,colour=alpha("black",0.5)) + xlab(C[e]~(mmol/m^3)) + ylab(q[e]~(mmol/m^3)) print(NC60.iso.graph) dev.off() Regards/Cordialement ------------- Benoit Boulinguiez Ph.D student Ecole de Chimie de Rennes (ENSCR) Bureau 1.20 Equipe CIP UMR CNRS 6226 "Sciences Chimiques de Rennes" Avenue du Général Leclerc CS 50837 35708 Rennes CED...