data(stjd)
genes = seq(100)
geneLocations = abs(stjd$geneLocation[genes])
position = c(1000) # location for which you want the kernelweights
kernelweights = rbf(as.matrix(geneLocations), as.matrix(position),
list(gamma=1/10e13))
hist(kernelweights)
Run the code above in your browser using DataLab