setwd("d:/nevs docs/rwork/test/")
library (adegenet)
library (RODBC)
channel <- odbcConnectAccess("Linkage_fe.mdb")
alldata <- sqlFetch(channel, "Toadegenet")
mydata<-subset(alldata,Site_Name!="ANW")
n<-ncol(mydata)
geno=mydata[,c(7:n)]
Site=mydata[,c(1)]
BirdID=mydata[,c(2)]
XY=mydata[,c(5:6)]
mygenind<-df2genind(geno, sep="",  ind.names=BirdID,ncode=3,pop=Site,loc.names=NULL,   missing=NA, ploidy=2, type=c("codom"))
mygenind@other$xy<-XY
myspca<-spca(mygenind)