Hello,
I'm trying to do some non-linear regression with 2 cell types and 4 tissue
type treatments using selfStart models
Following Ritz and Streibig (2009), I wrote the following routines:
##Selfstart
expDecayAndConstantInflowModel <- function(Tb0, time, aL, aN, T0){
exp(-time*aL)*(T0*aL+(-
Hi,
I've been trying to work with the gnls() function in the "nlme" package. My
decision to use gnls() was so that I could fit varPower and such to some of
the data. However, in working with a small dataset, I've found that the
results given by gnls() don't seem to make any sense and they di
Hi Kingsford,
That's exactly what I was looking for!
Thanks.
Mike
On Thu, 22 Oct 2009, Kingsford Jones wrote:
help(gnls, pack=nlme)
hth,
Kingsford Jones
On Thu, Oct 22, 2009 at 4:36 PM, Michael A. Gilchrist wrote:
Hello,
I've been fitting a random effects model using nl
Hello,
I've been fitting a random effects model using nlme to some data, but I am
discovering that the variation in my random effect is very small. As a result,
I would like to replace it as a fixed effect (i.e. essentially fit the same
model but with no random effect).
As I understand it I
On Sat, 17 Oct 2009, Douglas Bates wrote:
On Fri, Oct 16, 2009 at 7:09 PM, Michael A. Gilchrist wrote:
Hello,
I'm having problems creating and using a selfStart model with nlme.
Briefly, I've defined the model, a selfStart object, and then combined them
to make a selfStart.def
Hello,
I'm having problems creating and using a selfStart model with nlme. Briefly,
I've defined the model, a selfStart object, and then combined them to make a
selfStart.default model.
If I apply getInitial to the selfStart model, I get results. However, if I try
usint it with nls or nlsL
i, 16 Oct 2009, Dieter Menne wrote:
Michael A. Gilchrist wrote:
-
nlme(Count ~ quad.PBMC.model(aL, aN, T0),
+ data = tissueData,
+ weights = varConstPower(form =~ Count),
+ start = list( fixed = c(rep(1000, 8), -2, -2) ),
+ fixed = list(T0
Hello,
Excuse me for posting two questions in one day, but I figured it would be
better to ask my questions in separate emails. I will again give the caveat
that I'm not a statistician by training, but have a fairly decent
understanding of probability and likelihood.
As before, I'm trying t
Hello,
I will start out with the caveat that I'm not a statistician by training, but
have a fairly decent understanding of probability and likelihood.
Nevertheless, I'm trying to fit a nonlinear model to a dataset which has two
main factors using nlme. Within the dataset there are two Type ca
008 at 1:39 PM, Michael A. Gilchrist <[EMAIL PROTECTED]> wrote:
Wow, that's elegant and simple. It's also faster than my approach.
NB, you don't need to use close(), read.delim() closes the pipe when its
done reading.
If read.delim() close the connection in this case, it
k.edu/gilchrist.asp
-
On Fri, 12 Sep 2008, Prof Brian Ripley wrote:
Why not use
con <- pipe(COMMAND)
foo <- read.delim(con, colClasses="numeric")
close(con)
? See the 'R Data Input/Output Manual'.
On Fri, 12 Sep 2008, Michael A. Gilchrist wrote:
Hello,
I am currentl
Hello,
I am currently using R to run an external program and then read the results
the external program sends to the stdout which are tsv data.
When R reads the results in it converts it to to a list of strings which I
then have to maniuplate with a whole slew of commands (which, figuring out
ECTED]
[mailto:[EMAIL PROTECTED] On Behalf Of Michael A.
Gilchrist
Sent: Tuesday, August 26, 2008 6:03 PM
To: r-help@r-project.org
Subject: [R] switch() usage
Hello,
I am wondering if it is possible to run separate lines of
code within switch()
as one can do in most languages.
For example, I'd
Hello,
I am wondering if it is possible to run separate lines of code within switch()
as one can do in most languages.
For example, I'd like to have switch execute a few lines of code based on the
value passed to it.
For example,
index = 1
switch(index,
command 1; command 2; comm
14 matches
Mail list logo
