similar to: combining data from 2 independent surveys

Good afternoon! My question is more of a "is there a package for doing........?" an inference regarding the mean, median, regression estimates, etc. by using the information from dual frame surveys? (this methodology is based on the work of BANKIER 1986, SKINNER 1991, LOHR and RAO 2000...) If one has 2 independent studies which both independently measure a trait of a

Good morning! I have 4 questions which trouble me: 1. I want to test the hypothesis that the 2 proportions (the mean of a binomial) which come from 2 different samples are equal. I want to use the following function z= (p1-p2)/ sqrt((p1(1-p1)/n1)+(p2(1-p2)/n2)) which is one of the standard formulas for this case. Is there such a function in R? p1=the proportion from the first sample n1=the

Hi R users! I am trying to learn some multilevel analysis, but unfortunately i am now very confused. The reason: http://www.ats.ucla.edu/stat/hlm/seminars/hlm_mlm/mlm_hlm_seminar.htm http://www.ats.ucla.edu/stat/sas/seminars/sas_mlm/mlm_sas_seminar.htm and MlmSoftRev. pdf from mlmRev package. >From what i see, the first two links seem to declare the level one variable as a random part (i

Hi R users I want to use the lme package for a multilevel analysis on the following example: > math<-c(2, 3,2, 5, 6 ,7 , 7) > sex<-c(1, 2, 1, 2, 2, 2, 1) > school_A<-c(1,1,1,2,2,2,2) > school_B<-c(10,10,10,20,20,20,20) > mydata<-data.frame(math, sex, school_A, school_B) > mydata School_A and school_B are two different school characteristics, math is an

Hello R-users, Checking the archives, I recently came across this topic: "export tables to Excel files" (http://r.789695.n4.nabble.com/export-tables-to-Excel-files-td1565679.html#a1565679), and the following interesting references have been proposed: http://rwiki.sciviews.org/doku.php?id=tips:data-io:ms_windows

I never heard of Asterisk before today, but from what i'm looking at on the website and hearing, it sounds pretty incredibly. If I understand correctly with a 1,500.00 Wildcard TE410p T1 card, a good BSD or Linux Server, and a couple IP phones or Netmeeting on a few workstations, and of course, Asterisk which is free; I call have a small call center. This can't be? I was looking at

Ein eingebundener Text mit undefiniertem Zeichensatz wurde abgetrennt. Name: nicht verf?gbar URL: <https://stat.ethz.ch/pipermail/r-help/attachments/20130611/09dcb017/attachment.pl>

Hi list, Using the lmList function from ?nlme package? I get the following error message: ?Error in !unlist(lapply(sum.lst, is.null)) : invalid argument type? # this is the syntax used # fm_cream <-lmList(MULTDV~TIME|cod, data=zz, na.action=na.omit) # summary(fm_cream) When I generate a simple example the function works, but on the real data I get that message. I imagine that it is a

Hi List, I have the following problem: I am using the multilevel package and make.univ function for available in the package and then xyplot from lattice and I want to know how could I be able to use the ?coefficient? for the straight line that passes the data ? Example from help: library(multilevel) data(univbct) #a data set already in univariate or stacked form for job satisfaction

Dear all, I wanted to let you know about four phylogenetic comparative methods (PCM) packages that have become available on (3 on CRAN and 1 on GitHub) recently that hopefully will be interesting to somebody. Three of them go significantly beyond the Brownian motion (BM) and Ornstein-Uhlenbeck (OU) processes. 1) There is a new version of mvSLOUCH available. The most important change is that the

Dear all, I wanted to let you know about four phylogenetic comparative methods (PCM) packages that have become available on (3 on CRAN and 1 on GitHub) recently that hopefully will be interesting to somebody. Three of them go significantly beyond the Brownian motion (BM) and Ornstein-Uhlenbeck (OU) processes. 1) There is a new version of mvSLOUCH available. The most important change is that the

Hi R-users!! I have some trouble with the survey pakage and i would be very glad if you can give me an advice. I have a sample from a survey where household were interviewed. The sample has 4 criteria on which the stratification was based: REGION, SIZE OF HOUSEHOLD, SIZE OF LOCALITY, AGE OF HEAD OF HOUSEHOLD. Since i don't have the whole information in each cell of the cross

Hi R-users! Can anyone tell me what this error message refers to exactly > a<-c('m','f','m') > b<-c(1:3) > st<-c(1:3) > data.frame(a,b,st) a b st 1 m 1 1 2 f 2 2 3 m 3 3 > data.frame(a,b,st)->l > p<-svydesign(ids=~1, data=l) > pop<-c('(intercept)'=100, st1=10,st2=20) > calibrate(p,~st, pop) Error in

Hi List, Can anyone tell me how could i put the "BACK" button in the following code, just under the "AAA" menu? I want this button to go back to the previous page, and since it has nothing to do with the "1" and "2" buttons, i want it somehow separated from these two buttons, but i don't know how. I searched the web for some examples but my results

This is partly an R and partly a general statistics question. I'm trying to get confidence intervals of proportions (sometimes for subgroups) estimated from complex survey data. Because a function like prop.test() does not exist for the "survey" package I tried the following: 1) Define a survey object (PSU of clustered sample, population weights); 2) Use svyglm() of the package

Hello, I still apologise if the question i ask is naive. I would like to use R instead of usual softs (SAS, Stata) or at list to do things that the others do not. As a consequence I am trying to get familiar with R, and i am astonished because it does not seem to do things that are easy to do on other softs. I have got data with missing values (I am working on a questionnaire filled by students)

Hi, I have some problems to configure the xen I''ve installed (3.3.1). The computer is a Intel Core 2 Duo, I''m using Ubuntu 8.10 and have linux in my dom0 and winxp pro in my domU. I have two cores and I''d like to set the affinity of dom0 to cpu0 and domU to cpu1 but I haven''t find the way of making this permanent. I''ve set cpus=1 in the domU config

Good afternoon! I'm trying to use the Survey package for a stratified sample which has 4 criteria on which the stratification is based. I would like to get the corrected weights and for every element i get a weight of 1 E.g: tipping design <- svydesign (id=~1, strata= ~regiune + size_loc + age_rec_hhh + size_hh, data= tabel) and then weights(design) gives

Hi R users! I have the following problem: how appropriate is my aov model under the violation of anova assumptions? Example: a<-c(1,1,1,1,1,1,1,1,1,1,2,2,2,3,3,3,3,3,3,3) b<-c(101,1010,200,300,400, 202, 121, 234, 55,555,66,76,88,34,239, 30, 40, 50,50,60) z<-data.frame(a, b) fligner.test(z$b, factor(z$a)) aov(z$b~factor(z$a))->ll TukeyHSD(ll) Now from the aov i found that my model

Greetings, I have the following? Problem: Given k (=10) discrete independent random variables X_i with n_i (= 5 to 20) values each,compute quantiles of the distribution of the sum X = X_1+...+X_k. Here X has n=n_1 x n_2 ... n_k distinct values which is too large to list them all together with their probabilities. I tried several approaches: (A) Convolution: each X_j is approximated with