similar to: plotrix 2.8-3

Hi, I've been looking at ways to make pyramid plots in R. I like the pyramid.plot method in plotrix as it seems the simplest to use and building them in ggplot looks a bit more code intensive than I'd like, being as I'm new to R. This package does pretty much what I need it to do, however I can't seem to scale the x axis labels. The other labels scale fine with labelcex. I

Dear R Users, Greetings! I am confused with an error message related to "plotrix" library in the newer versions of R. I used to run an R script without fail in the earlier versions (R 2.8.1) of R; but the same script is now throwing up an error message in the newer versions (Now I have R 2.13.0 & R 2.14.0). Herewith I am furnishing the same code for your perusal. It would

I am using pyramid.plot() from the plotrix package. I have something like this ############################################ xy.pop<-dados$masfr xx.pop<-dados$femfr #agelabels<-dados$femlab xycol<-color.gradient(c(0,0,0.5,1),c(0,0,0.5,1),c(1,1,0.5,1),11) xxcol<-color.gradient(c(1,1,0.5,1),c(0.5,0.5,0.5,1),c(0.5,0.5,0.5,1),11) xylab<-dados$maslab xxlab<-dados$femlab

Hi everybody, An update to plotrix (2.5-3) is now on CRAN. Apart from the usual polishing, this includes an interesting new function (hierobarp) that was suggested on this very list a few weeks ago. This function will plot nested breakdowns of numeric values by one or more factors. It seems quite useful for visualizing the deviations in successively finer divisions of the data. Comments are

Hi all, After a lot of rewriting, I am happy to say that plotrix_2.6-4 is on CRAN. The reason I am announcing this is that the newish hierobarp function has had a complete rewrite. This was the result of my learning that the old version was horrendously slow on large data sets. Let me know if there are any other problems in the package. Next up - a function to center plot titles to the figure

Hi all, After a lot of rewriting, I am happy to say that plotrix_2.6-4 is on CRAN. The reason I am announcing this is that the newish hierobarp function has had a complete rewrite. This was the result of my learning that the old version was horrendously slow on large data sets. Let me know if there are any other problems in the package. Next up - a function to center plot titles to the figure

Hi all Another question related to my earlier -- how might I get transparency to work with pyramid.plot in library(plotrix)? I can't figure it out, and the help page says it is there but doesn't give an example. Thanks!

Dear List: I am working to understand some differences between the results of the svymean() function in the survey package and from code I have written myself. The results from svymean() also agree with results I get from SAS proc surveymeans, so, this suggests I am misunderstanding something. I am never comfortable with "I did what the software" does mentality, so I am working to

Dear WizaRds, Working through sampling theory, I tried to comprehend the concept of stratification and apply it with Survey to a small example. My question is more of theoretic nature, so I apologize if this does not fully fit this board's intention, but I have come to a complete stop in my efforts and need an expert to help me along. Please help: age<-matrix(c(rep(1,5), rep(2,3),

R-help, Sorry for this question (I guess it has been addressed before but I could not find it in the archives) but how can I see a function code when the following comes up: > svymean function (x, design, na.rm = FALSE, ...) { .svycheck(design) UseMethod("svymean", design) } <environment: namespace:survey> Thanks in advance

Hi, I have been experimenting with the new Survey package. Specifically, I was trying to use some of the functions on the public-use survey data from NHIS (2000 Sample Adult file). Error 1): The first error I get is when I try to specify the complex survey design. nhis.design<-svydesign(ids=~psu, probs=~probs, strata=~strata, data=nhis.df, check.strata=TRUE) Error in svydesign(ids =

Hello, Although my R code for the svymean () and svyquantile () functions works fine, I am stuck with the svyby () and make.formula () functions. I got the following error messages. - Error: object of type 'closure' is not subsettable # svyby () - Error in xx[[1]] : subscript out of bounds # make.formula () A reproducible example is appended below. I would appreciate if

I have the following code which produces the output below it clus1 <- svydesign(ids = ~schid, data = lower_dat) items <- as.formula(paste(" ~ ", paste(lset, collapse= "+"))) rr1 <- svymean(items, clus1, deff='replace', na.rm=TRUE) > rr1 mean SE DEff W525209 0.719748 0.015606 2.4932 W525223 0.508228 0.027570 6.2802 W525035 0.827202

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)

Dear WizaRds, dear Thomas, First of all, I want to tell you how grateful I am for all your support. I wish I will be able to help others along one day the same way you do. Thank you so much. I am struggling with a multistage sampling design: library(survey) multi3 <- data.frame(cluster=c(1,1,1,1 ,2,2,2, 3,3), id=c(1,2,3,4, 1,2,3, 1,2), nl=c(4,4,4,4, 3,3,3, 2,2), Nl=c(100,100,100,100,

Dear list, I got some surprising results when using the svytotal routine from the survey package with data containing missing values. Some example code demonstrating the behaviour is included below. I have a stratified sampling design where I want to estimate the total income. In some strata some of the incomes are missing. I want to ignore these missing incomes. I would have expected that

Hi R users, I am trying to use the aggregate function with a survey design object and survey functions, but get the following error. I think I am incorrectly using the syntax somehow, and it may not be possible to access variables directly by name in a survey-design object. Am I right? How do I fix this problem? I have used aggregate with "mean" and "weighted.mean", and

How do I retrieve the replicates estimates from a crosstab done using svyby? Here is an example from the help page for svyby in the package: > data(api) > dclus1<-svydesign(id=~dnum, weights=~pw, data=apiclus1, fpc=~fpc) > rclus1<-as.svrepdesign(dclus1) > > a <- svyby(~api99, ~stype, rclus1, svymean, return.replicates=TRUE) > a$replicates NULL But, compare to > b

Hi all, When I try the following with pkg Survey it returns the error below: ftable(svyby(~INCOME, ~AGECL+RACECL, svymean, design=q50), rownames=list(AGECL=c("<35", "35-44", "45-54", "55-64", "65-74", ">=75"), RACECL=c("white non hispanic", "non white or

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