setwd("d:/nevs docs/rwork/test/") library (adegenet) library (RODBC) #channel <- odbcConnectAccess("Linkage_fe.mdb") #mydata <- sqlFetch(channel, "ToadegenetNoANW") #odbcClose(channel) mydata<-read.csv("mydata.csv") n<-ncol(mydata) geno=mydata[,c(7:n)] #geno <- as.matrix(dat) # conver to matrix #y <- which(is.na(dat)==TRUE) # get index of NA values #geno[y] <- "000/000" # replace all na values with "000/000" #geno Site=factor(mydata[,c(3)]) 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<-df2genind(geno, sep="/", ind.names=BirdID,ncode=3,loc.names=NULL, missing="NA", ploidy=2, type=c("codom")) #change line above to include populationon data mygenind@other$xy <-XY myspca<-spca(mygenind)