myprojects-hg-reborn

Fork

rm(list=ls())
data<-read.csv("data-ven215.csv",header=FALSE)
data<-as.matrix(data)

datapar<-read.csv("data-par.csv",header=FALSE)

datapar<-as.matrix(datapar)

f <- function(x, A1,mu1,myvar1) {
    expr <- expression(
log(10)*A1/sqrt(2*pi)/log(myvar1)*exp(-((log(x/mu1))^2)/2/log(myvar1)/log(myvar1)))
    eval(expr) ### expression I want to use for my NLS
}
print("OK")

N<-data[ ,2]
Dp<-data[ ,1]
pdf("plots-venturi-215.pdf")
plot(Dp, N, bg = "black", pch = 1, cex = 1.5,xlab="Dp [nm]",
ylab=expression("dN/dlogDp ["*cm^-3*"]"),
main=expression("Venturi Mass Flow Rate 21.5 "*m^3/min),log="x"
,ylim=range(c(0,1.8e7)))
N<-data[ ,3]
 lines(Dp,N,"p", bg = "black", pch = 21, cex = 1.)
lines(Dp,N,"p", bg = "black", pch = 1, cex = 1.5)

N<-data[ ,4]
lines(Dp,N,"p", bg = "black", pch = 1, cex = 1.5)


 mypar<-c(datapar[ 1,1],datapar[2,1],datapar[3,1])

newx<-seq(min(Dp),max(Dp),len=1000)


g <- function(x, A1,mu1,myvar1) {
    expr <- expression(
log(10)*A1/sqrt(2*pi)/log(myvar1)*exp(-((log(x/mu1))^2)/2/log(myvar1)/log(myvar1)))
    eval(expr) ### expression I want to use for my NLS
}


M1 <- do.call("g", c(list(x = newx), mypar))

lines(newx, M1, col="blue", lwd=2.5)

mypar<-c(datapar[ 1,2],datapar[2,2],datapar[3,2])

M1 <- do.call("g", c(list(x = newx), mypar))

lines(newx, M1, col="orange", lwd=2.5)

mypar<-c(datapar[ 1,3],datapar[2,3],datapar[3,3])

M1 <- do.call("g", c(list(x = newx), mypar))

lines(newx, M1, col="green", lwd=2.5)
dev.off()