similar to: R-squared value for linear regression passing through origin using lm()

Hi, I'm using a non linear model to fit experimental survival curves. This model describes the fraction of "still active" experiments as a function of time t as follows: f(t)=(1+exp(-etaD*cD)) / (1+exp(etaD(t-cD))) Moreover, when experiments are still active, they may change of state (from 0 to 1). But they may fall inactive before changing their state (their state still

I am wondering if there is any procedure to write a particular value in Percentage format, still maintaining it's numeric character. for example I want to write '.33' as '33%' regards, thanks in advance --------------------------------- [[alternative HTML version deleted]]

Hi, I'm analyzing experimental results where two different events ("T1" and "T2") can occur or not during an experiment. I made my experiments with one factor ("Substrate") with two levels ("Sand" and "Clay"). I would like to know wether or not "Substrate" affects the occurrence probability of the two events. Moreover, for each

Hi, I guess there have been discussions about this in the past and from what I understood hooking an R-function to facilitate automatic adjustment is problematic. So why not doing it like this: --- R-3.4.1/src/unix/sys-std.c 2017-03-24 00:03:59.000000000 +0100 +++ R-3.4.1/src/unix/sys-std.patched.c 2017-08-28 09:16:02.714204023 +0200 @@ -1005,6 +1005,9 @@ // introduced in

? stato filtrato un testo allegato il cui set di caratteri non era indicato... Nome: non disponibile Url: https://stat.ethz.ch/pipermail/r-help/attachments/20070917/f6d8d606/attachment.pl

Dear List, Is there any way to password-protect script files (either within R or otherwise)? platform i386-pc-mingw32 arch i386 os mingw32 system i386, mingw32 status major 2 minor 5.1 year 2007 month 06 day 27 svn rev 42083 language R version.string R version 2.5.1 (2007-06-27)

Hi, All, I am a beginner for R. Now I have installed R 2.5.1 in Window environment. After I run a program such as "gam" I would like to display a plot for the object. The following is an example. When I did this, only the last plot was presented on my screen. How can I get a plot before the last plot? I mean if the object has several plots how can I get those? "gam.object <-

I feel like an idiot posting this because every language I've ever seen has a string function that trims blanks off strings (off the front or back or both). Ideally, it would process whole data frames/matrices etc but I don't even see one that processes a single string. But I've searched and I don't even see that. There's a strtrim function but it does something completely

I have following syntax for putting a legend : legend("bottom", fill=c("red","blue"), legend=expression(p==0.30, p==0.50), bty="n") However what I want is that : the value "0.30" should be a value of a variable instead of a constant, so that I can put the name of this variable and in legend it's value will be displayed. Can anyone tell me how

Hi there, I've got a small question : is there any post-hoc test for Kruskal rank sum test integrated in R ? I know that the Nemenyi test is one of the post-hoc that can be used, but there's no (to by knowledge) R function for it. What should I do ? Thanks, Etienne

Is there an easy way to open an R script file in the R Editor found in more recent versions of R on Windows via a file association? I am looking for functionality akin to how the ".ssc" file extension works for S-Plus: upon double-clicking on a ".R" file, Rgui.exe starts up and loads the script file in the R Editor. As far as I can tell, Rgui.exe does not have a command line

Hi, Sorry for the naive question, but what exactly does the 'Adjusted R-squared' coefficient in the summary of linear model adjust for? Sample code: > x <- rnorm(15) > y <- rnorm(15) > lmr <- lm(y~x) > summary(lmr) Call: lm(formula = y ~ x) Residuals: Min 1Q Median 3Q Max -1.7828 -0.7379 -0.4485 0.7563 2.1570 Coefficients:

All, I am using a sum of squared differences in the objective function of an optimization problem I am doing and I have managed to speed it up using the outer function versus the nested for loops, but my suspicion is that the calculation could be done even quicker. Please see the code below for a simple example. If anyone can point out a faster way I would appreciate it greatly. Thanks, Roger

Hi, I have some problems when using AlgDesign->optFederov() generating designs. I have 6 variables, all factors. 3^2 and 4^4, I want to have a design that can take care of main effects and two interactions within 2 pair of variables v3-v4 and v5-v6, the following is the code ################ require(AlgDesign) set.seed(1) levels = c(v1=3,v2=3, v3=4,v4=4,v5=4,v6=4)

Hi, I've made a research about how to compare two regression line slopes (of y versus x for 2 groups, "group" being a factor ) using R. I knew the method based on the following statement : t = (b1 - b2) / sb1,b2 where b1 and b2 are the two slope coefficients and sb1,b2 the pooled standard error of the slope (b) which can be calculated in R this way: > df1 <-

Hi, sorry if this is explained somewhere but I didn't find anything. How can I use "Anova" from the car package to test a modell without between subject's factors? Suppose I have the following data mat.1 mat.2 mat.3 di ex 1 85 85 88 1 1 2 90 92 93 1 1 3 97 97 94 1 1 4 80 82 83 1 1 5 91 92 91 1 1 6 83 83

Oups, I sent the email by error, as I was still writing my reply… Spencer, Le 22-avr.-09 à 03:33, spencerg a écrit : > Is your first model a special case of the second with eta1 = 0? > If yes, what about using 2*log(likelihood ratio) being approximately > chi-square? Yes, the first model is a special case of the second with eta1=0… Could you give me more explanation about

On Fri, 1 Sep 2017, Ralf Goertz wrote: > Am Mon, 28 Aug 2017 09:33:31 +0200 > schrieb Ralf Goertz <r_goertz at web.de>: > > > Hello, me again > >> Hi, >> >> I guess there have been discussions about this in the past and from >> what I understood hooking an R-function to facilitate automatic >> adjustment is problematic. So why not doing it

On Thu, Nov 29, 2018 at 1:10 PM S Ellison <S.Ellison at lgcgroup.com> wrote: > > > > > plot(x=1:10, y=) > > > plot(x=1:10, y=, 10:1) > > > > > > In both cases, 'y=' is ignored. In the first, the plot is for y=NULL (so not > > 'missing' y) > > > In the second case, 10:1 is positionally matched to y despite the

On Sat, 15 Dec 2018, frederik at ofb.net wrote: > I agree with Steve and Achim that we should keep some examples with no > data frame. That's Objectively Simpler, whether or not it leads to > clutter in the wrong hands. As Steve points out, we have attach() > which is an excellent language feature - not to mention with(). Just for the record: Personally, I wouldn't recommend