relativeweight <- function(n1=1,n2=1,sigsqalpha=1,sigsqerror=1)
{
theweight <- (n1/n2)*(n2*sigsqalpha + sigsqerror)/(n1*sigsqalpha +
                                                      sigsqerror)

return(theweight)
}

### Plot a function over a grid

gridcalc <-function(n1lim=10,n2lim=10,sigsqalpha=1,sigsqerror=1){
x <- c(1:n1lim)
y <- c(1:n2lim)
return(outer(x,y,FUN=relativeweight,sigsqalpha,sigsqerror))
}



showthis <- function(n1lim=10,n2lim=10,sigsqalpha=1,sigsqerror=1){
x <- c(1:n1lim)
y <- c(1:n2lim)
persp(x,y,outer(x,y,FUN=relativeweight,sigsqalpha,sigsqerror))
return()
}