Displaying 1 result from an estimated 1 matches for "logs647".
2004 Feb 13
5
predict function
...on procedure.
In 1.5.1 I used the following commands to normalise the variable "logratio",
over a 2d surface (defined by coordinates x and y):
> array <- read.table("121203B_QCnew.txt", header=T, sep="\t")
> array$logs555<-log(array$s555)/log(2)
> array$logs647<-log(array$s647)/log(2)
> array$logratio<-array$logs555-array$logs647
> array$logav<-(array$logs555+array$logs647)/2
> library(modreg)
> loess2d<-loess(logratio~x+y,data=array)
> array$logratio2DLoeNorm <-array$logratio - predict.loess(loess2d, array)
However in 1.8.1...