similar to: CROSSTABULATION

I have this assignment to do and after ten hours of constant trying my eyes ache and i give up.. all i'm able to get is this plot please help me these are the commands i have used till now read.table(file.choose(), sep=";", header=T) read.table(file.choose(), sep=";", header=T)->areas melt(areas,id=c("Year","State"),m=c("Rice"))->

Hello list members, I'm looking for a way to force the results of a crosstabulation to be square - that is, to include 0 values. For example: table(letters[1:4],letters[c(1:3,3)]) yields: a b c a 1 0 0 b 0 1 0 c 0 0 1 d 0 0 1 I would like to return: a b c d a 1 0 0 0 b 0 1 0 0 c 0 0 1 0 d 0 0 1 0 Any suggestions? Thanks! -- Manuel A. Morales

Hi everyone, I am a new one for that lists....actually i have final year project on VOIP & IMS ...so i want to install asterisk on my pc ...IS it possble that ...we can call on small LAN network without buying any card...i will clear my point as that...suppose i have a linux machine on which i want to install asterisk and HOW it will install...and second point is that ..i have 2

hi , i have downloaded the asterisk@home software from the web ..but i have a little confusion about that ...either i wrote in blank cd as it is or some bootable media is required for it...as it is in zip format...BUT it is a .ISO file ...tell me ...what should i do...it will run automatically when i reboot system and first boot device is CDROm...thankssss -- Muhammad Sohail Arham U.E.T. Lahore

Hello R-users, My toy example: aa<-c(1:5) bb<-c(NA,2,NA,4,5) cc<-c(1,2,NA,4,NA) dd<-c("A","B","B","A","C") df<-data.frame(aa,bb,cc,dd=as.factor(dd)) table(unlist(df[,1:3])) Can anyone point me to what function let's me do a crosstabulation between table(unlist(df[,1:3])) and df$dd? I want to find out when dd==A (or B,

HI people , I am not sure whether this is right mailing list for this question , but i tried subscribing to redhat rpm list , failed to do so . I'm facing a issue while creating RPM for python punjab http://code.stanziq.com/punjab/ I've created a spec file which contains the list of files to packaged http://pastie.org/1236027 And when i try to run rpmbuild -ba command on this spec

how about this: tapply(vector,list(factor1,factor2),function)? -----Original Message----- From: owner-r-help at stat.math.ethz.ch [mailto:owner-r-help at stat.math.ethz.ch]On Behalf Of Murray Jorgensen Sent: 24. september 2000 02:46 To: R-help Subject: [R] Crosstabulation I can't seem to find a function in R similar to Splus crosstabs() for creating a multi-way table from factors and a

I'm havinf a problem with a simple file i have the following data State 1960 1970 1980 1990 1 All India 35988.70 37346.00 39707.30 42321.00 2 Andhra Pradesh 3431.03 3163.27 3687.23 3695.63 3 Assam 1902.93 2001.60 2278.47 2525.33 4 Bihar 5277.07 5133.80 5138.70 4662.57 5 Gujarat 538.13 456.10 484.23 590.47 6

Hello, I have a dataframe (var_1, var_2, ..., var_n) and I would like to export summary statistics to Latex in the form of a table. I want specific summary statistics by crossing numerous variables 2x2 AT ONCE. In each cell I would like sometimes to have the median (Q1 - Q3), or frequency and proportion, etc. CrossTable, xtab, etc... do not allow for multiple 2 by 2 crosstabulation. The table

Hi, I find difficult to read crosstabulated data without percentages, so I wrote this dirty function that may be useful to others. It only works up to 3 dimensions since I am not very good at programming in R. I would appreciate if someone else has a better similar function or can improve this one. As an example of use: > x<-rbinom(100,1,.3) > y<-rbinom(100,1,.3) >

Hello All, Is there a way of producing an ANOVA table split into contrasts, thus showing the contrasts sums of squares and associated p-values? Thanks, Martin. Martin Hoyle, School of Life and Environmental Sciences, University of Nottingham, University Park, Nottingham, NG7 2RD, UK Webpage: http://myprofile.cos.com/martinhoyle

Hi all, Has anyone tested for FA in R? I need to seperate out the variance due to measurement error from variation between individuals (following Palmer & Strobeck 1986). Andy Higginson Animal Behaviour and Ecology Research Group School of Biology University of Nottingham NG7 2RD U.K. This message has been checked for viruses but the contents of an attachment may still contain

Hi all, I'm having some troubles with the Quantreg package. I am using R version 2.10.0, and have downloaded the most recent version of Quantreg (4.44) and SparseM (0.83 - required package). However, when I try to run an analysis (e.g. fit1<-rq(y~x, tau=0.5)) I get an error message saying that the function "rq" could not be found. I get the same message when I try to search

Dear R Users, I suppose this is a school boy question, but here it is anyway. I'm trying to re-create the residuals for a poisson GLM with simulated data; x<-rpois(1000,5) model<-glm(x~1,poisson) my.resids<-(log(x)- summary(model)$coefficients[1]) plot(my.resids,residuals(model)) This shows that my calculated residuals (my.resids) are not the same as residuals(model). p 65 of

Hello, I am trying to print a crosttabulation of mean,sd,n for a continuous variable crossclassified by anoother/s grouping variables. I came up with: xtab2 <- function(x,g1,g2) { funy <- function(z) list(mean(z,na.rm=T),sd(z,na.rm=T),length(z)) aa <- by(x,list(g1,g2),funy) bb <- matrix(unlist(aa),nrow=3 ,dimnames=list(c("mean","sd","n"),

Hello R Users, I am investigating the basic use of the LME function, using the following example; Response is Weight, covariate is Age, random factor is Genotype model.lme <- lme (Weight~Age, random=~ 1|Genotype) After summary(model.lme), I find that the estimate of Age is 0.098 with p=0.758. I am comparing the above model with the AOV function; model.aov <- aov (Weight~Age + Genotype)

Dear All, I'm rather a beginner on nested ANOVAs, so here goes with my 2 questions; Qu 1: I'm modelling the number of galls on a leaf (the response variable) as a function of; the tree on which I find the leaf, the branch on which I find the leaf. Then, the tree and the branch are both random factors, and I'm quite happy that I should write; aov(galls~tree/branch +

Hi list, Is there a standards based reason why the DIFlags enum is set to uint32_t[1]? I am sure my DWARF-std-reading-fu is not up to snuff and so I cannot seem to find it. The reason I ask is that we are running out of space for our own DIFlags and would like to nail this down before deciding on an approach. Thanks! Sohail [1] The code in question:

Dear All, I wonder if anyone can advise me as to whether there is a consensus as to how the effect size should be calculated from GLIM models in R for any specified significant main effect or interaction. In investigating the causes of variation in infection in wild animals, we have fitted 4-way GLIM models in R with negative binomial errors. These are then simplified using the STEP procedure,

Hello! I have dataset where one variable is 'multiresponse', like this: [1] "1 2" "1 2 3" "4" "1 4" "4 3" etc. 'responses' separated by space. observations in different 'rows' of data.frame. I can do strsplit(data$var,' ') and make a list, where multiple responses are elements of character vectors, like this: $