R/getNe.R

"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
}

Try the optiSel package in your browser

Any scripts or data that you put into this service are public.

optiSel documentation built on Jan. 16, 2021, 5:30 p.m.