R/getNe.R

In optiSel: Optimum Contribution Selection and Population Genetics

"getNe"<-function(Year, Div,df=3.5,I){
	for(j in 2:length(Div)){if(is.na(Div[j]) || Div[j]>Div[j-1])Div[j]<-Div[j-1]}
	DivEst <- c(predict(sm.spline(Year, Div, df=df), Year, 0))
	DivAbl <- c(predict(sm.spline(Year, Div, df=df), Year, 1))
	deltaf <- -I*DivAbl/DivEst
	
	Ne       <- 1/(2*(deltaf))
	names(Ne) <-Year
	Ne[Ne<0]<-NA
	Ne
}