[R] boot: how is t0 calculated?

Dear all, I am doing log-transformed bootstrap regression using: x=c(0.038,0.054,1.030,1.250,10.240,52.000,228.100,240.000,758.000,1502.600,0.140,9.170,280.000,1.000,0.150,0.388,20) y=c(3961.5,25987.5,526557,321094.5,642,23592715.5,4050, 228853521.1,320980392,42900,58435.5,13525240.5

[R] calculate RMSE using bootstrap

Dear all, I need to calculate confidence intervals of a and b and root mean square error (RMSE) of a power law given by Y = a X^b I calculate the confidence intervals by: reg.ln = function(tab.xy, indices){ tab.xy = storage[indices,] res.lm = lm(log(y)~log(x), data=tab.x

[R] fitting power growth

Dear all I am using the x and y vectors as defined below and want do to a power law regression: y = a x^b using > lm(log(y)~log(x)) gives reasonable values (b=1.23) but is not very popular due to biases of back-transformation from log to non-log values. Using > nls(y~a*x^b,start=list(a=1

[R] help with expression

Hi all, I still fail to plot an axis title with the following expression: plot(0,xlab=expression('(SOC [' * kgm^{-2} * '])' * ^{-2})) the xlab should look like: (SOC [kgm^2])^0.25 with an out bracket and a superscript. Thanks for your advice. Thomas [[alternative HTML version del

[R] determing font type in expression

Dear R-Maillist members, I would like to label the y-axis of a plot with a mathematical expression "sed.rate mm a-1", where the "-1" is printed subscript. I am using font.lab=15 for the x-axis. However, when plotting the y-axis with the expression-function it appears in a different font. Does

[R] creating expression with for-loop

Dear listmembers I would like to create an expression that looks like labl = expression(10^1,10^2,10^3,10^4,10^5) using a for-loop. However for (i in 1:5){ labl[i]=expression(10^i) } does not do the right thing. Does anybody knows help? Thanks Thomas

[R] filled.contour with log axis

Dear all, I would like to generate a filled.contour plot with log x and y axis, however using: filled.contour(as.line,log="xy") results in a warning message. Does anybody knos what to do? Thanks Thomas __ R-help@r-project.org mailing list https://

Re: [R] axis names in triangle.plot

ed row.names? > Please return the results of class(grain) and names(grain) > > Cheers, > PS: If you have questions related to ade4, you can use the adelist > (http://listes.univ-lyon1.fr/wws/info/adelist) > > > Thomas Hoffmann wrote: >> Hi folks, >> >> I am using

[R] axis names in triangle.plot

Hi folks, I am using a data.frame with sediment grain sizes: > grain sand silt clay OAT 10.03 56.77 18.25 OAT 10.40 57.40 17.94 WG1 50.03 20.68 12.57 WG1 43.20 25.69 13.41 WG1 33.89 31.10 14.48 WG2 2.84 62.81 20.79 WG2 2.79 60.46 19.16 WG2 16.27 33.04 6.48 WG2 1.39 57.90 9.13 WG3 4.

[R] fit to power law with inknown zero

Dear R-helpers I would like to do a fit of the form: y = a (x+c)**b, where a, b and c are unknown. Does anybody know how to do it? Thanks Thomas __ R-help@r-project.org mailing list https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read the post

[R] power law fit with unknown zero

Dear R-helpers I would like to do a fit of the form: y = a (x+c)**b, where a, b and c are unknown. Does anybody know how to do it? Thanks Thomas __ R-help@r-project.org mailing list https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read the post