search for: dipa

...ut.3932.100")) #print (snp); # pi: desired proportion of variation due to QTN pi = 0.05; print (paste("pi:", pi)); MAF = 0.05; print (paste("MAF:", MAF)); # S: number of segregating sites S = length(snp[1,]); # N: number of samples N = length(snp[,1]); Dips = sample(1:N,N) DipA = Dips[1:50] DipB = Dips[51:100] disnp = snp[DipA,]+snp[DipB,] snp = as.data.frame(disnp, row.names=NULL); N = length(snp[,1]); # get allele freq for all SNPs allele_f <- mean(snp[,1:S])/2; print (allele_f); sites = rfind(allele_f); print(sites); # collapse sites that have perfect correlation...

Hello Sir, I am using mgcv package for my data. My model is y~x1+f(x2),I want to find out the function f(x2) . Following is the code.   sm1=gam(y~x1+s(x2),family=binomial, f) summary(sm1) plot(sm1,residuals=TRUE, xlab="AGE",pch=20)   In this plot I am getting S(x2,1.93) on y axixs  How should I get the function for x2 from this plot.or Is there anyother procedure in R  to get this

I am using semiparametric Model  library(mgcv) sm1=gam(y~x1+s(x2),family=binomial, f) How should I  find out standard error for ed50 for the above model ED50 =( -sm1$coef[1]-f(x2)) / sm1$coef [2]   f(x2) is estimated value for non parametric term.   Thanks [[alternative HTML version deleted]]