Mercurial > repos > mnhn65mo > regionalgam
diff gls.R @ 0:5b126f770671 draft
Uploaded
author | mnhn65mo |
---|---|
date | Fri, 03 Aug 2018 08:28:22 -0400 |
parents | |
children |
line wrap: on
line diff
--- /dev/null Thu Jan 01 00:00:00 1970 +0000 +++ b/gls.R Fri Aug 03 08:28:22 2018 -0400 @@ -0,0 +1,16 @@ +#!/usr/bin/env Rscript +library(nlme) +library(MASS) + +args = commandArgs(trailingOnly=TRUE) +input1 = read.table(args[1], header=TRUE) #input1=col.index =glmmpql-output +input2 = read.table(args[2], header=TRUE,sep=" ") #input2=data.index =abundance_index-output + +input1<-as.matrix(input1) + +year <- unique(input2$YEAR) +mod <- gls(input1~year) +summary<-summary(mod) + +save(mod, file = "mod.rda") +capture.output(summary, file="mod-summary.txt")