[R-gregmisc-users] SF.net SVN: r-gregmisc: [943] trunk/PathwayModeling/thesispaper/R
Brought to you by:
warnes
From: <r_b...@us...> - 2006-03-20 11:49:54
|
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. |