Displaying 6 results from an estimated 6 matches for "wtloss".
Did you mean:
loss
2008 Jul 31
1
nls weights warning message
...w what is causing it and how I can fix it.
As an example, the following, from Venables and Ripley, which does not
produce the warning message and then on the data which does procude the
warning message.
~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
##Venables and rippley
library(MASS)
data(wtloss)
wtloss.st = c(b0=90, b1=95, th=120)
n1<-length(wtloss)
Weight=c(1,rep(0.1,n1-1))
wtloss.fm = nls(Weight ~ b0 + b1*2^(-Days/th), data=wtloss, start=wtloss.st,
trace=T, weights=Weight )
wtloss.fm
##My data
y<-c(73,73,70,74,75,115,105,107,124,107,116,125,102,144,178,149,177,124,157,128,
169,1...
2002 Mar 19
0
MASS Book Exercise
...th)
model.func <- b0 + b1 * temp
# Differentiate temp
derive.negexp <- (x * temp) / th^2
Z <- cbind(1, temp, derive.negexp)
dimnames(Z) <- list(NULL, c("b0","b1","th"))
attr(model.func, "gradient") <- Z
model.func
}
wtloss.gr <- nls(Weight ~ expn(b0, b1, th, Days),
data = wtloss, start = wtloss.st, trace = F)
expn2 <- deriv(~b0 + b1*((w0 - b0)/b1)^(x/d0),
c("b0","b1","d0"), function(b0, b1, d0, x, w0) {})
wtloss.init <- function(obj, w0) {
p &...
2001 May 22
2
MASS data sets
I'm running R 1.2.2 under windows 98 on a Pentium 133 laptop.
I can't seem to retrieve the package MASS data sets:
> library(MASS)
> data(wtloss)
Warning message:
Data set `wtloss' not found in: data(wtloss)
> data(abbey)
Warning message:
Data set `abbey' not found in: data(abbey)
And yet all the .rda files for the MASS datasets are in
D:\Program Files\rw1022\library\MASS\data
Help, anyone?
Murray Jorgensen, Department...
2001 Sep 30
1
MASS data sets
All the usual MASS data sets are nicely sitting in
D:\Program Files\rw1022\library\MASS\data
on my laptop
and yet:
> library(MASS)
> data(wtloss)
Warning message:
Data set `wtloss' not found in: data(wtloss)
what could be wrong? I'm running 1.2.2 on Win98.
Murray Jorgensen
-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-
r-help mailing list -- Read http://www.ci.tuwien.ac.at/~hornik/R/R-FAQ.htm...
2003 Nov 07
1
y label after axis (4)
Hi,
I am trying to figure out how to lable the second y-axis after the following
codes:
plot(x, y,
xlab="time",
ylab="pay1"
)
par(new=TRUE)
plot(x,y2,
ann=FALSE,
xaxt="n",
yaxt="n",
pch=7
)
axis(4)
Then, I want to label the second y axis "pay2". I tried "title(ylab="pay2")",
but it put this lable on the first
2003 Apr 09
1
'Apparently' trouble with name spaces and Sweave...
...s, nottem, npk, npr1, oats, OME, painters, petrol, phones,
Pima.te, Pima.tr, Pima.tr2, quine, Rabbit, road, rock, rotifer, Rubber,
ships, shoes, shrimp, shuttle, Sitka, Sitka89, Skye, snails, SP500,
steam, stormer, survey, synth.te, synth.tr, topo, Traffic, UScereal,
UScrime, VA, waders, whiteside, wtloss
Error in library(MASS) : package/namespace load failed
Error in driver$runcode(drobj, chunk, chunkopts) :
Error in library(MASS) : package/namespace load failed
Error in buildVignettes(dir = ".") : Error in driver$runcode(drobj,
chunk, chunkopts) :
Error in library(MASS) :...