pfx <- file.path(system.file("extdata", package="plinkFile"), "m20.rel")
rel <- readBSM(pfx) # relatedness kernel matrix
re2 <- rel^2 # 2nd order polynomial kernel
tmp <- tempdir()
dir.create(tmp, FALSE)
out <- file.path(tmp, 'm20.re2')
saveBSM(out, re2) # save the polynomial kernel
dir(tmp) # show new files, then clean up
unlink(tmp, recursive=TRUE)
Run the code above in your browser using DataLab