Displaying 2 results from an estimated 2 matches for "yjg45".
2004 Oct 23
1
Legend/Substitute/Plotmath problem
Hello,
I seem unable to construct a legend which contains a substitution as
well as math symbols. I'm trying to do the following:
strain2 <- "YJG48"
legend.txt <- c(
substitute(
strain *
%==% *
"YJG45, rpn10" *
%Delta%,
list(strain=strain2)
),
"Verhulst/Logistic",
"Malthus"
)
legend(
100,2.5,
legend.txt,
cex=0.75,
bty="n",
pch=c(20,NA,NA),
lty=c(NA,1,2)
I derived this from the following "text" command which works,
text(
160,2.2,
cex=...
2004 Aug 25
3
Beginners Question: Make nlm work
...or conceptual problem which can't be solved by me staring at it ...
Since I'm repeating something I have done with gnuplot I know that 3 and
4e-3 as starting values for the fit are appropriate ...
Thanks for any hint,
Joh
setwd("~/Biology/R_versuch")
mydata<-read.table("YJG45-7_Growth.dat")
x<-mydata$V1
y<-mydata$V2
VH <- function(p) y ~ p[1]/(1+((p[1]-0.008)/0.008)*exp(-(p[2]*x)))
plot(x, y, xlab="Time (h)",ylab=expression(OD[600][~nm]),las=1)
out <- nlm(VH, p = c(3, 4e-3), hessian = TRUE)