similar to: why doesn't table() have a data=argument?

In the following, I'm fitting a logistic regression model, and using car:::influencePlot. When I run the latter with output to the screen, it calls identify() that lets me label observations with large CookD. However, if I use postscript() to get .eps output, identify() seems not to be called at all. If instead, I use dev.copy2eps() after getting output to the screen, the point labels

I'm pretty sure that I did everything right, but my R is just not drawing the mosaicplot that I want, and there is also no error messege, looks like this right now: > mosaicplot(arthritis) > mosaicplot(~ sex + treatment + improved, data = arthritis, color = TRUE) > -- View this message in context: http://r.789695.n4.nabble.com/problems-with-mosaic-plot-tp3090128p3090128.html Sent

This is my reproducible example: example<-structure(list(SENSOR = structure(1:6, .Label = c("A", "B", "C", "D", "E", "F"), class = "factor"), VALUE = c(270, 292.5, 0, 45, 247.5, 315), DATE = structure(1:6, .Label = c(" 01/01/2010 1", " 01/01/2010 2", " 01/01/2010 3", " 01/01/2010

In http://tolstoy.newcastle.edu.au/R/e2/help/06/10/3064.html a method was given for converting a frequency table to an expanded data frame representing each observation as a set of factors. A slightly modified version was later included in the NCStats package, only on http://rforge.net/ (and it has too many dependencies to be useful). I've tried to make it more general, allowing an input

Hoping someone can help me with xtabs and ftable. I'm trying to get a pair of ftables (possibly more) next to each other. For example: > dunhill_lights_xtab<-ftable(xtabs(grossedupobs ~ gender+age_group + dunhill_lights, data = ciggs)) > dunhill_lights_xtab dunhill_lights No Yes gender age_group Female

For the following model, library(vcd) arth.polr <- polr(Improved ~ Sex + Treatment + Age, data=Arthritis) summary(arth.polr) where Improved is an ordered, 3-level response I'm looking for a *simple* way to test the validity of the proportional odds assumption, typically done via a score test for equal slopes/effects over the predictors. I do find a po.test= option in the repolr package

Hi all, Please, is there any way of controlling factors in row/columns when using ftable/xtabs? As far as I can see, the last cross-clasifing variable in the formula will appear in columns. The previous ones, in rows. For instance, is it possible to make tension and replicate appear in columns? ftable(xtabs(breaks ~ wool + tension + replicate, data = warpbreaks)) After some years using SAS

Hi, I am trying to reproduce some functionalities of Excel pivot table in R, sadly, I couldn't figure out how to do it. I am wondering if this is even possible in R. Does anyone know? Here is an example: year=rep(2003,16) quarter=rep(1:4,each=4) sales=1:16 company=rep(c("a","b","c","d"),4) df=data.frame(year,quarter,sales,company) #this is the

I am getting assembler errors on clang-llvm-2.9 output for a program with exceptions that I do not get when using the installed g++ Mac OS-X 10.4 powerpc-apple-darwin > as -version Apple Computer, Inc. version cctools-590.23.2.obj~17, GNU assembler version 1.38 the folks at gnu-binutils assure me this is an assembler bug, but also that this isn't a recognizable "gnu"

Dear list, I would like to do some calculation using different grouping variables. My 'df' looks like this: # Some data set.seed(345) id <- seq(200,400, by=10) ids <- sample(substr(id,1,1)) group1 <- rep(1:3, each=7) group2 <- rep(1:2, c(10,11)) group3 <- rep(1:4, c(5,5,5,6)) df <- data.frame(id, ids, group1, group2, group3) df <- rbind(df, df, df) df$time <-

Dear all, I have a dataset. I want to make barplot from this data. Zero1 <- " V1 V2 V3 V4 V5 V6 V7 V8 V9 1 1 0 0 0 1 0 0 0 Positive 2 0 0 1 0 1 0 1 1 Negative 3 0 0 1 0 0 0 1 1 Positive 4 0 1 0 1 1 1 0 1 Negative 5 0 0 1 0 1 1 0 0 Positive 6 0 1 0 0 1 1 1 1 Negative 7 1 0 1 1 1 1 1 1 Negative 8 0 0 0 0 1 0 0 1

Dear listmembers, I have a little question: I have my data organized as follow sic level area a 211 2.4 b 311 2.3 b 322 0.2 b 322 0.5 c 100 3.0 c 100 1.5 c 242 1.5 d 222 0.2 where levels and sics are factors. I'm trying to obtain a matrix like this: level 211 311 322 100 242 222 sic a 2.4 0 0

Hello R-helpers. I've tried to recreate a parallel version of tapply() and table() using a combination of the parallel functions mclapply() and pvec() and papply(), but haven't been successful. In the end, I'm trying to get a cross tab of two vectors. I currently (can) use tapply(..., sum) and table(), and even xtabs() and ftable(), but with tens of millions of words and tens of

I have a list with four fields and 3500 records, from a MySQL table. How can i have the total for each level using table, ftable, xtabs? ex: print(ftable(table[,"field1"]~table[,"field2"])) thanks Danilo Maurizio Veneto Lavoro - Area Osservatorio - -.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.- r-help mailing list -- Read

Hello, I'm new to R and have been following many guides including the two-way anova (http://www.personality-project.org/r/r.anova.html). Using that walkthrough including the supplied data I do str(data.ex2) and receive the appropriate types of data as follows: > str(data.ex2) 'data.frame': 16 obs. of 4 variables: $ Observation: int 1 2 3 4 5 6 7 8 9 10 ... $ Gender :

-- apologies if this is a dup - i got a bounce saying the message was unprocessed. Is there a straightforward way to get a table with percents in the cells rather than counts? I've looked at table, ftable, xtabs, and ctab, which did the conversion but returned the results in a single row without labels. any suggestions are appreciated. thank you.

Problem on flexmix when trying to apply signature developed in one model to a new sample. Dear R Users, R Core Team, I have a problem when trying to know the classification of the tested cases using two variables with the function of flexmix: After importing the database and creating a matrix: BM<-cbind(Data$var1,Data$var2) I see that the best model has 2 groups and use: ex2

I am quite new to R, so please bear over with me if I have problems with the R terminology. I want to (try to) use R for some analyses within vegetation ecology, using the vegan package. I have my data in a postgresql database, and I manage to get them into R as a dataframe with columns for respectively: Name of the analysed m2, Name of the species, coverage of species in the square in %. I

I am using R 2.0.2 on a WinXP I am trying to get the code of the Kruskal-Wallis test but > kruskal.test function (x, ...) UseMethod("kruskal.test") <environment: namespace:stats> > ls(3) [1] "acf" "acf2AR" "add.scope" .............. [181] "kruskal.test" "ks.test"

Can anyone help with scripting the commands below? 1. Add/create a new column to df and name it as ShortPause . 2. Convert ShortPause as a factor variable. 3. ShortPause will have two levels (short and long). Assign short if pause is smaller than 0.1, and assign long otherwise. 4. How many ?syntax==intransitive? observations were realized with a tap (tapped==Tapped)? (Use xtabs() and ftable() to