similar to: R Code for X-Tab with Row/Col Proportions, Expected Vals and Tests

i have a data frame(dat) which has many variables.and i use the following script to get the crosstable. >danx2<-c("x1.1","x1.2","x1.3","x1.4","x1.5","x2","x4","x5","x6","x7","x8.1","x8.2","x8.3","x8.4","x11",

Hi, In Splus if I call the function crosstabs() the output is a contigency table; in each cell of the table is printed: N, N/RowTotal, N/ColTotal, N/Total. N is the number of observations in each cell. The same call to xtabs() in R will produce the contigency table but the only entry in each cell is N. How can I get the same relative frequencies that crosstabs() gives? Thanks, mike --

Hello R-devel, I am currently attempting to implement an API similar to data.table wherein single bracket subsetting can accept an unquoted expression to be evaluated in the context of my object. A simple example from the data.table package looks like this: DT <- data.table(col1 = c('a', 'b', 'c'), col2 = c('x', 'y', 'z')) DT[col1 ==

Hi, I need to compare two matrices with each other. If you can get one of them out of the other one by resorting the rows and/or the columns, then both of them are equal, otherwise they're not. A matrix could look like this: [,1] [,2] [,3] [,4] [,5] [,6] [,7] [,8] [1,] 0 1 1 1 0 1 1 0 [2,] 1 1 0 0 0 1 0 1 [3,] 1 0 1 0 0

Hi All, I have been able to modify the x-axis to start at zero by adding xlow and xhigh parameters; that was pretty simple. I have been unable to find the location of the code that would turn off the information weighting of the box size (I have smaller randomized trials getting less weight than a much larger non-randomized trial). The function is forestplot() from rmeta. Thanks for any

hi, is there a way of calculating of measuring dependence between two categorical variables. i tried using the chi square test to test for independence but i got error saying that the lengths of the two vectors don't match. Suppose X and Y are two factors. X has 5 levels and Y has 7 levels. This is what i tried doing >temp<-chisq.test(x,y) but got error "the lengths of the two

Library gmodels include a function CrossTable that is useful for crosstabulation. In the help, it is indicated that one can call this function as CrossTable(data), were data is a matrix. However, when I try to use this option, it doesn't help. Any idea? Is there a bug? Thanks for your help. Prof. Albert Sorribas Grup de Biomatem?tica i Bioestad?stica Departament de Ci?ncies M?diques B?siques

Has anyone used an AJAX data grid with Rails? I am rebuilding a copy of the venture source web site as a first ruby/ rails exercise. If anyone has used a javascript data grid with Rails, I would appreciate sample code, or notes on your experience. The controller action would presumably communicate with the data grid component running on the browser by extracting data from a relational

I have written a wrapper for heatmap.2 called heatmap.w.row.and.col.clust which auto-generates breaks using breaks<-round((c(seq(from=(-20 * stddev), to=(20 * stddev))))/20, digits = 2) #(stddev in this case = 2.5) This has always worked well in the past but now I am getting an error that non-finite breaks are being generated. Drilling down, it seems that my wrapper is generating finite

Hi all, After a long time having no time, I have finally done some fresh work on my software fax machine. I have replaced the original carrier tracking with something more robust. I have also added 4800, and 2400 bits per second modes, and cleaned up a few bugs in areas like superfine mode operation. I apologise for this update taking so long. At ftp://ftp.opencall.org/pub/spandsp you will

I re-wrote the compile script to build klibc (cleaning up how I do it, and how to find the source tarballs). It seems to be compiling fine. I tried several of the commands in usr/utils/static and most one. However, the sleep command fails in a strange way (EINVAL on execve): ================================================================ faraday/phil /home/phil/klibc-project-0 581> strace

Hi all, I have a 3x4 contingency table with row totals all being 100. I want to generate 3 x 4 tables from the null distribution. Which R function can do this? -- Thanks, Jim. [[alternative HTML version deleted]]

Rxfax answers, makes handshake, and crashes once the page starts to send. It receives a .tif file of 8 bytes. Asterisk dumps core - gdb shows : #0 0x281fd86c in t4_rx_putbit () from /usr/local/lib/libspandsp.so.0 #1 0x281fea3c in fast_putbit () from /usr/local/lib/libspandsp.so.0 #2 0x28208324 in decode_baud () from /usr/local/lib/libspandsp.so.0 #3 0x2820893f in process_baud () from

Sorry to bore you more with the clock issue, but have you check /proc/zaptel/<span> to make sure it's not missing interrupts? There's also an option to record the audio for the fax, you could listen to that vs a recorded file that will receive correctly on a fax machine and see whether there is an obvious difference? (Good luck, that'll be really scraping the barrel!!) Does it

Hi, I just sent this to Steve Underwood, but then found a bunch of posts on the mailing list about similar issues.. does anyone have the fix? I'm running asterisk CVS-HEAD-06/28/04-18:13:13, spandsp 0.0.1k, libtif 3.5.7 one thing i just noticed is that calls come in with format '72' which is G711A-law or LinearPCM.. it uses PCM for the call, i assume this is ok the results of RxFAX

I'm suppressing the labeling of my rows and columns in heatmap.2 using the commands: labRow = " ", labCol = " " But I'd like to annotate them again using the axis command: axis(1, at=seq(500, 1000, 500)) mtext("Group 1", "Group 2") For some reason however it appears that the axis command is having no effect. My complete

Hello, I'm trying to learn how to use lattice and levelplot in particular. There are three elements of customizing the plots I'm stuck with: a) Is there a way to put borders around each "cell" within a level-plot. I'm trying to do something like the colsep/rowsep/sepcolor/sepwidth parameters of heatmap.2 in gplots b) Can I alter the line-width of dendrogram added

Full_Name: John Hendrickx Version: 1.6.0 OS: Windows 98 Submission from: (NULL) (137.224.174.216) CrossTable in the "gregmisc" package fails when the fisher.exact test produces an error (I suspect this is because the number of cases is too large). This can be fixed using "FTt <- try(fisher.test(t, alternative = "two.sided"))" or by making the test optional.

What I am doing : create a simple correlation matrix on 41 variables, then plot an heatmap with this program : library(gplots) mat<-cor(temp.alln,use="pairwise.complete.obs") hm<-heatmap.2(mat,symm=T) HM<-format(round(mat[hm[[1]],hm[[2]]],2)) library(RColorBrewer) brewer.pal(10,"Spectral")->mp heatmap.2(mat, symm = TRUE, col = mp,

Dear R community, I'm currently using R to analyze functional Magnetic Resonance Imaging data. Each analysis involves running ~120,000 repeated-measures ANOVAs. I would like to know if there is any automatic way to access the F- and P-value data that are associated with each of these 120,000 ANOVAs. For example, if the summary output (for the 1st ANOVA of 120,000) shows the following value