justin bem
2006-Nov-14 07:40 UTC
[R] Variance of a complex estimator using survey package ...
Un texte encapsul? et encod? dans un jeu de caract?res inconnu a ?t? nettoy?... Nom : non disponible Url : https://stat.ethz.ch/pipermail/r-help/attachments/20061114/6c1fa862/attachment.pl
justin bem
2006-Nov-14 14:02 UTC
[R] Variance of a complex estimator using survey package ...
Un texte encapsul? et encod? dans un jeu de caract?res inconnu a ?t? nettoy?... Nom : non disponible Url : https://stat.ethz.ch/pipermail/r-help/attachments/20061114/e4439e11/attachment.pl
Thomas Lumley
2006-Nov-14 15:25 UTC
[R] Variance of a complex estimator using survey package ...
On Tue, 14 Nov 2006, justin bem wrote:> Hi, > > I want to compute the variance of two complex statistics. The first > statistic is the a ratio > R1=Q(a1)/Q(a2) > where Q denote de quantile at a1 and a2. The second is also a ratio but > not a classic one this ratio is > R2=sum(x_{i}|x_{i}>Q(a1))/sum(x_{i} }|<Q(a2)) > > Linearisation for those statistic is very complex. I want to use > bootstrap and JKn How can I procede since I dont have a function like > svrepanyfunction where you just specify the function to compute the > statistic.You can specify any function of data and weights as an argument to withReplicates(). JKn weights will not work for the first statistic (they don't work for quantiles), but bootstrap weights should work. -thomas