Revision: 943
Author: r_burrows
Date: 2006-03-20 03:49:41 -0800 (Mon, 20 Mar 2006)
ViewCVS: http://svn.sourceforge.net/r-gregmisc/?rev=943&view=rev
Log Message:
-----------
minor update
Modified Paths:
--------------
trunk/PathwayModeling/thesispaper/R/plotDensities.R
trunk/PathwayModeling/thesispaper/R/plotVi.R
Modified: trunk/PathwayModeling/thesispaper/R/plotDensities.R
===================================================================
--- trunk/PathwayModeling/thesispaper/R/plotDensities.R 2006-03-20 11:41:09 UTC (rev 942)
+++ trunk/PathwayModeling/thesispaper/R/plotDensities.R 2006-03-20 11:49:41 UTC (rev 943)
@@ -6,5 +6,6 @@
hist(output[,i],breaks=30,xlim=c(low,high),xlab=paste("d",i,sep=""),main="")
par(new=T)
plot(seq(low,high,step),dnorm(seq(low,high,step),modes[i],SD[i]),type='l', xlim=c(low,high),xlab="",yaxt='n',ylab="")
+ abline(v=mcmcML[i],col=2)
}
}
Modified: trunk/PathwayModeling/thesispaper/R/plotVi.R
===================================================================
--- trunk/PathwayModeling/thesispaper/R/plotVi.R 2006-03-20 11:41:09 UTC (rev 942)
+++ trunk/PathwayModeling/thesispaper/R/plotVi.R 2006-03-20 11:49:41 UTC (rev 943)
@@ -3,7 +3,8 @@
# attach(as.data.frame(input))
times <- input[,1]
CI <- get95CI(2,input,1,Dsamp)
- temp <- fitPaper(paramsML$par,input)[[1]]
+ temp <- fitPaper(optimML,input)
+ temp2 <- fitPaper(mcmcML,input)
low <-min(CI,input[,3])
high <- max(CI,input[,3])
plot(times,input[,3],ylim=c(low,high),xlab="time",ylab="v2")
@@ -12,12 +13,13 @@
par(new=T)
plot(times,CI[,2],ylim=c(low,high),type="l",lty=2,xlab="",ylab="",col=2)
par(new=T)
- plot(times,temp,ylim=c(low,high),type="l",xlab="",ylab="",col=2)
+ plot(times,temp[[1]],ylim=c(low,high),type="l",xlab="",ylab="",col=3)
+ par(new=T)
+ plot(times,temp2[[1]],ylim=c(low,high),type="l",xlab="",ylab="",col=2)
abline(h=0)
legend(32,2000,c("data","fitted curve","95% CI"),lty=c(0,1,2),pch=c(1,NA,NA),col=c(1,2,2))
for (i in 3:5) {
CI <- get95CI(i,input,1,Dsamp)
- temp <- fitPaper(paramsML$par,input)[[i-1]]
low <-min(CI,-400)
high <- max(CI,800)
plot(times,input[,i+1],ylim=c(low,high),
@@ -27,7 +29,9 @@
par(new=T)
plot(times,CI[,2],ylim=c(low,high),type="l",lty=2,xlab="",ylab="",col=2)
par(new=T)
- plot(times,temp,ylim=c(low,high),type="l",xlab="",ylab="",col=2)
+ plot(times,temp[[i-1]],ylim=c(low,high),type="l",xlab="",ylab="",col=3)
+ par(new=T)
+ plot(times,temp2[[i-1]],ylim=c(low,high),type="l",xlab="",ylab="",col=2)
abline(h=0)
}
# detach(as.data.frame(input))
This was sent by the SourceForge.net collaborative development platform, the world's largest Open Source development site.
|
