similar to: use of lm() and poly()

Hi all, I'm using 'expression' in the axes labels in graphs to get the greek letters (ylab), as follows: plot( layer_1$U0 ,layer_1$dP/(H1.Theo) #,ylim=c(0,100) ,xlim=c(0,1.1) #,main="Pressure Drop", ,xlab=expression(U[0]~(Nm/s)) ,ylab=expression(Delta~P~(Pa/m)) ,col=2 ,pch=16 ,cex=1.5) I just need to know how to avoid the use of a space (~) between the greek

Hi all, I try to write a simple function in a script. The script is as follows yo<-function(Xdata) { n<-length(Xdata[,1]) Lgm<-nls(formula=LgmFormula, data=Xdata, start=list(a=1500,b=0.1),weights=Xdata$Qe) return(Lgm) } After the execution of the script, when I call the function yo on data called NC60.DATA I get an error. #yo(NC60.DATA) Erreur dans eval(expr, envir, enclos)

Hi all, Is it possible to change the filling of the squares used to represent the colour legend in a bar plot with ggplot? in this example, fillings are raven black, I'd like them white. ggplot(diamonds, aes(clarity, colour = cut)) + geom_bar() Regards -- ------------- Benoit Boulinguiez Ph.D student Ecole de Chimie de Rennes (ENSCR) Bureau 1.20 Equipe CIP UMR CNRS 6226 "Sciences

Hi all, Is that possible to remove the grey colour in the legend key that goes with the geom_smooth? In my case it doesn't ease the reading of the legend. http://www.4shared.com/file/125864977/e10644f8/desorb.html Cordialement / Regards ------------------------------------------- Benoit Boulinguiez Ecole de Chimie de Rennes (ENSCR) Bureau 1.20 Equipe CIP UMR CNRS 6226 "Sciences

Hi all, I've got a list (Reg3Lst) with 1000 nls regression results in it. I'd like to get the confidence interval of the parameters obtained with the nonlinear regressions. Thus I've used this: for (i in 1:1000) { foo<-list(foo,confint(Reg3Lst[[i]])) } For some regressions the confidence interval is not estimated because of a singular gradient result. Then it

Hi all, I'm smoothly transferring my lattice graphs to ggplot2 graphs, but I'm stuck on representing a curve from a formula. I'm looking for the equivalent of curve() in ggplot2, Hadley Wickham mentions geom_curve, but as far as I've seen in the help it doesn't exist. My need is to plot a regular scatter plot of experimental data (easy to do actually) and then add the

Hi all, quick question: is the optional command "width" effective in the geom_errorbarh() layer of ggplot? Cause I can't get it works on this graph http://www.4shared.com/file/116919103/93488d88/iso_2PrsH.html pdf(file = "iso_2PrsH.pdf", width = 7, height = 7) NC60.iso.graph<-ggplot( NC60.DATA ,aes(Ce,Qe)) + geom_point(col=MaCouleur1, size=4) +

Hi all, I've this ggplot2 graph (http://www.4shared.com/file/120101043/f3e9350/isotherm.html). I'd like to tweak the legend for the adsorbent part. I'd like to have full square of the colors instead of the full circles. I guess that I've to work on the scale_color_hue or fill somewhere but I haven't found my way yet. Here is my code to get the graph:

Hi all, I need to realize nonlinear regression on a thousand data sets. I guess the lapply function would help me on that thus I'd like to create a list of data frames, each data frame containing the data as follows: Ce Qe 1 1.849147 0.1958672 2 10.054250 0.5771036 3 18.077246 0.7718514 4 27.576468 0.8079606 5 35.146862 0.8500489 6 43.245078 0.8366673 7 51.745760 0.8879672

Hi all, for publication purpose I must provide a graph with grid lines aligned with each tick mark. Thing is that the graph has a log scale on the x-axis. I looked at the grid() and par() documentation but still don't figure out how to get it done. simple example plot(seq(1,9), log="x", panel.first=grid()) how to get lines at x = 2 and 5 on the graph? Of course the option

Hi, Quite simple question I guess. Where could we specify the size of the text (axis,legend) for a scatter plot graphic made with plot()? 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 35700 Rennes, France Tel 33 (0)2 23

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

Hi all, I try to represent a multiple curve graphic where the x-axis is the temperature and the different y-axes are the different X (X22,X43,X44...) some X corresponds to the same molecule (22 and 44 are for CO2 for instance) so I use the same colour for them. I wanna mix the linetype with the colour to be able to visually see the difference between X43 and X45 The best I have done up to now

Hi, Parameters assessment in R with nls doesn't work, though it works fine with MS Excel with the internal solver :( I use nls in R to determine two parameters (a,b) from experimental data. m V C0 Ce Qe 1 0.0911 0.0021740 3987.581 27.11637 94.51206 2 0.0911 0.0021740 3987.581 27.41915 94.50484 3 0.0911 0.0021740 3987.581 27.89362

Hi all, I try to assess the parameters (K1,K2) of a model that describes the adsorption of a molecule onto on adsorbent. equation: dq/dt = K1*C*(qm-q)-K2*q I know the value of 'qm' and I experimentally measure the variables 'q', 'C', and the time 't'. t C q 1 0 144.05047 0.0000000 2 565 99.71492 0.1105625 3 988 74.99426

Hi all, Is that possible to get back to the par() default values easily? I played a bit with par() 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 35700 Rennes, France Tel 33 (0)2 23 23 80 83 Fax 33 (0)2 23 23 81 20

Hi all, I divide the graphical device into 4 ---mfrow=c(2,2)--- to plot 3 boxplots graph. Thus I get an empty space for the fourth graph. I d'like to place a legend in this space. legend(xxxx) creates a legend in the last graph (the 3rd) I'd like a legend in the empty space. How may I do it? #graph parameter par( fin=c(6,6), #dim graph mai=c(1,1,0.5,0.5) #marges graph

Dear R users, I seek for a more elegant manner to manipulate my data that that I produced so far. Data is in a data frame THC515.DATA -sample at the end of this mail- of 5 variables, a string and 4 continuous numeric variables. I need to get the mean of two variables "Q" and "T" at the levels of a third one "dP", individually for each level of "n". I

Hi, I'd like to know whether R is capable to assess parameters in a model describing the kinetic of a pollutant adsorption onto activated carbon. A common relation is for instance the Adam-Bohart-Thomas' one: dx/dt = K1 * (qm-x)*C - K2x where {K1,K2} are the unknown paramters and {qm,C} are known parameters Of course I get experimental data sets of measured x as a function of time.

Hello there, I am using nlme() to fit a logistic mixed effect model on our data. The outcome variable is binary. I got the error when I wanted to add a group factor to my model. My initial model is as below: model.a <- nlme(response~ 1/(1 + exp( -intercept- u0 - slope*TIME - u1*TIME)), + fixed=intercept+slope~1, random= u0+u1~1 |studentID,