similar to: svychisq??

Hola Buen día, estoy ajustando una regresión multiple, pero sucede que tengo un factor de expanción, digamos que deseo ajustar el modelo ingreso ~ edad + ocupacion mi problema es que cada observación tienen un factor de expanción, es decir cada variable no representa una observación sino más, que no es un número fijo para cada observación. por lo que a mí se me ocurrio usarlas como peso

Hello I want to know if the summary of the logistic model with survey Pr (> | t |) to test if the coefficient of the model is significant, ie is the p_valor wald test for the model coefficients, for I am interested to know if the three levels of the variable educational level are significant to the model (significance of handling 0.2), I present below the results of my model

Hello, I?m hoping you have a few minutes to help out someone very new to R. I?ve done some searching, but cannot find this particular issue. I have survey data from two different time periods (years). Both years are stratified samples and have the same variables (and variable names), but are different people in the community answering in the different years. Everything loads into the survey

Dear All I am trying to use svychisq with a two-dimensional table 4 x 5. The command I am using is summary(svytable(~dietperception+dietstatus,dudleyls1rake,na.rm=TRUE),"C hisq") It is throwing up an error message as follows: Error in NCOL(y) : only 0's may be mixed with negative subscripts In addition: Warning messages: 1: is.na() applied to non-(list or vector) in:

Dear Forum I a running svychisq from the survey package and get errors with the number of dimensions, errors that I do not understand and do not know how to fix. I ask you kindly to help me out. The eror message follows with some information below. I hope there are enough information to help me to fix the problem if not please let me know what is needed Best Stefan Hrafn Jonsson > QHISQ2

Thank you so much Rui. 1. How do I export this table to excel file? I used this tbl1 <- table(Country, IDNum) tbl2=addmargins(tbl1) write.xlsx(tbl2,"tt1.xlsx"),sheetName="summary", row.names=FALSE) The above did not give me that table. 2. I want select those unique Ids that do have records in all countries. From the above data set, this ID

Thank you so much Rui! On Sun, Feb 4, 2018 at 12:20 AM, Rui Barradas <ruipbarradas at sapo.pt> wrote: > Hello, > > Please always cc the list. > > As for the question, I believe the following does it. > > a <- strsplit(mydata$ID, "[[:alpha:]]+") > b <- strsplit(mydata$ID, "[[:digit:]]+") > > a <- sapply(a, `[`, 1) > c <-

Hi R-help, I posted about this late yesterday but got no response. I may have put TMI in the original request. Not to mention I couldn't cut and paste yesterday because I was working R off a non-network computer while asking for help on a network computer. Essentially, I have this user-defined function: test <- function(X){ chisq <- svychisq(~X

Hello, Thanks for the reproducible example. See if the following does what you want. IDNum <- sub("^(\\d+).*", "\\1", mydata$ID) Country <- sub("^\\d+(.*)", "\\1", mydata$ID) tbl1 <- table(Country, IDNum) addmargins(tbl1) tbl2 <- xtabs(Y ~ Country + IDNum, mydata) addmargins(tbl2) Hope this helps, Rui Barradas On 2/3/2018 3:00 AM, Val

Hi all, I have a data set need to be summarized by unique ID (count and sum of a variable) A unique individual ID (country name Abbreviation followed by an integer numbers) may have observation in several countries. Then the ID was changed by adding the country code as a prefix and new ID was constructed or recorded like (country code, + the original unique ID Example original ID

Hola Abro una sesión de R y me aparecen unos objetos de una sesión de ayer. Los borro, cierro R y vuelvo a abrir otra sesión y los tengo otra vez. Es la primera vez que me pasa. ¿sabe alguien por qué me pasa esto y cómo me puedo deshacer de ellos, aunque sean míos y les tenga mucho cariño? Saludos Jesús [[alternative HTML version deleted]]

R-help, I'm trying to evaluate the misclasification rate of an estimated discrimated function (hyDA) hyDA <- with(hyB, lda(Maturity ~ Length + Weight)) , so here is what I do (code from S-plus 6.1 Guide to Statistics Vol. 2 p 101) test <- predict(hyDA) tbl <- table(hyB$Maturity, test$class) tbl <- tbl[-2,-2] tbl1 <- cbind(tbl, error = (apply(tbl, 1, sum) - diag(tbl)) /

Dear R-listers, I want to reperfrom a poisson distribution test that presented in a recent-published biological research paper (Plant Physiology 2008, vol 148, pp. 1189-1200). That test is about the occurrence number of a kind of gene in separate chromosomes. For instance: The observed gene number in chromosome A is 36. The expected gene number in chromosome A is 30. Then, the authors got a

Hi Guys I am back with another thing that's puzzling me. I am creating contingency tables but then I want to filter out certain columns and also find if any entry in the table is 0. Example: gts labels A1 B2 G3 1 21 127 120 2 23 112 0 Here I want to remove B2 column from this table and also if any entry is 0 in this case G3 second row. Missing out on how to do

Dear Md Kamruzzaman, To answer your second question first, you could just use the svychisq() function. The difference-of-proportion test is equivalent to a chisquare test for the 2-by-2 table. You don't say how you computed the confidence intervals for the two separate proportions, but if you have their standard errors (and if not, you should be able to infer them from the confidence

== What is Scratch? Scratch is a little toy I put together in Ruby on Rails. Think of it as a spin-off of my larger Rails app, Elite Journal. If nothing else, you can think of Scratch as a demonstration of the new Action Web Services component of Rails. From the web page: Scratch is the minimalist''s web log. Scratch gives you nothing more than the meta-weblog API for posting.

Thank you very much for the explanations, Josh. | the seekpoints are not the only places in the stream that can | be seeked to, they are just hints about regular locations in | the stream. you can still seek to any sample. the number of | seekpoints are just one factor in determining how fast the | decoder can find any particular sample. I did a lot of testing | before deciding on -S10s. So

Hola Pedro. El INE cambió los ficheros de microdatos no hace mucho, aquí dejo como se haría ahora, (utilizando MicroDatosEs). Lo que cambia es la función para recodificar. http://rpubs.com/joscani/unemplrate El 15/10/14 a las #4, Carlos Ortega escribió: > Hola Pedro, > > Acabo de recordar que hace poco José Luis Cañadas (participa en esta lista) > publicó un enlace suyo a un

Hello Everyone, I was analysing big survey data using survey packages on RStudio. Survey package allows survey data analysis with the design effect.The survey package included functions for all other statistical analysis except two-proportion z tests. I was trying to calculate the difference in prevalence of Diabetes and Prediabetes between the year 2011 and 2017 (with 95%CI). I was able to

Thank you for the new release, Josh. I've downloaded it, but some of the details in the history file puzzle me, so I've not tried to use it yet. Pardon my denseness here ... First, as long as you don't use the --cuesheet option when you encode, are .flac files encoded by 1.1.0 still readable by earlier versions? Second, if the old default for seekpoints was -S100x and the new one is