similar to: Trouble when I call a function within another

Displaying 20 results from an estimated 900 matches similar to: "Trouble when I call a function within another"

2008 Sep 18
1
axis text size in scatter plot
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
2008 Sep 23
1
superscript text style in axis label in graph
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
2008 Sep 10
3
writing simple function through script
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)
2008 Sep 30
1
par() default values
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
2008 Oct 17
1
use of empty space in split graph device
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
2008 Oct 15
1
parameter assessment in differential equation
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.
2009 May 26
2
using lsoda() and nls() together
Thanks to Dieter Menne and Spencer Graves I started to get my way through lsoda() Now I need to use it in with nls() to assess parameters I have a go with a basic example dy/dt = K1*conc I try to assess the value of K1 from a simulated data set with a K1 close to 2. Here is (I think) the best code that I've done so far even though it crashes when I call nls()
2008 Oct 30
1
continue a loop after an error with confint
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
2009 Jul 07
1
curve from a formula with ggplot2
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
2009 Jul 09
1
ggplot2: geom_errorbarh()
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) +
2009 Jul 23
1
tweak a bit the legend wiht ggplot2
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:
2009 Apr 28
1
text in graph with expression
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
2010 Mar 20
1
grid lines aligned on each ticks with a log scale
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
2009 May 13
2
ode first step
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
2008 Sep 03
2
nls convergence trouble
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
2008 Oct 27
1
create list of data frames
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
2009 Aug 18
1
ggplot2: geom_smooth and legend
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
2010 Sep 10
3
ggplot bar geom: control the filling in the colour legend
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
2009 Dec 22
2
use of lm() and poly()
Hi all, I want to fit data called "metal" with a polynominal function as dP ~ a.0 + a.1 * U0 + a.2 * U0^2 + a.3 * U0^3 + a.4 * U0^4 The data set includes, the independant variable U0 and the dependant variable dP. I've seen that the combination of lm() and poly() can do that instead of using the nls() function. But I don't get how to interpret the results from the linear
2009 Sep 09
2
ggplot2: mixing colour and linetype in geom_line
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