Calculate population allele frequencies from a vcf file including multiple populations

I have a vcf file with about 800 individuals (diploids) and millions of SNPs. The individuals can be divided in 15 to 25 populations.
I would like to calculate the allele frequencies for each SNP on each population.
Has someone got a R script doing this?
Thank you