On Fri, 2005-07-01 at 19:40 +0800, ronggui wrote:
> I write a function to get the frequency and prop of a variable.
>
> freq<-function(x,digits=3)
> {naa<-is.na(x)
> nas<-sum(naa)
> if (any(naa))
> x<-x[!naa]
> n<-length(x)
> ta<-table(x)
> prop<-prop.table(ta)*100
> res<-rbind(ta,prop)
> rownames(res)<-c("Freq","Prop")
> cat("Missing value(s) are",nas,".\n")
> cat("Valid case(s) are",n,".\n")
> cat("Total case(s) are",(n+nas),".\n\n")
> print(res,digits=(digits+2))
> cat("\n")
> }
>
> > freq(sample(letters[1:3],48,T),2)
> Missing value(s) are 0 .
> Valid case(s) are 48 .
> Total case(s) are 48 .
>
> a b c
> Freq 11.00 20.00 17.00
> Prop 22.92 41.67 35.42
>
> and i want the result to be like
> a b c
> Freq 11.00 20.00 17.00
> Prop 22.92% 41.67% 35.42%
>
> how should i change my function to get what i want?
Here is a modification of the function that I think should work. Note
that part of the output formatting process has to take into account the
a priori unknowns involving your 'digits' argument, the lengths of the
dimnames resulting from the table and the lengths of the frequency
counts in the table. Thus, a fair amount of the code is establishing the
'width' argument, which is then used in formatC() so that the output can
be column aligned properly.
Note that by default, table() will exclude "NA", so you do not need to
subset 'x' before using table().
Also, note that I change "Prop" to "Pct".
freq <- function(x, digits = 3)
{
n <- length(x)
missing <- sum(is.na(x))
ta <- table(x)
pct <- prop.table(ta) * 100
width <- max(nchar(unlist(dimnames(ta))) + 1,
nchar(ta) + digits + 1,
5 + digits)
Vals <- paste(formatC(unlist(dimnames(ta)), format = "s",
width = width),
collapse = " ")
Freq <- paste(formatC(ta, format = "f", digits = digits,
width = width),
collapse = " ")
Pct <- paste(formatC(pct, format = "f", digits = digits,
width = width),
"%", sep = "", collapse = " ")
cat("Missing value(s) are", missing, ".\n")
cat("Valid case(s) are", n - missing,".\n")
cat("Total case(s) are", n, ".\n\n")
cat("", Vals, "\n")
cat("Freq", Freq, "\n")
cat("Pct ", Pct, "\n")
cat("\n")
}
Thus:
> freq(sample(letters[1:3], 48, TRUE), 2)
Missing value(s) are 0 .
Valid case(s) are 48 .
Total case(s) are 48 .
abc
Freq 28.00 8.0012.00
Pct58.33% 16.67% 25.00%
> freq(sample(c(letters[1:3], NA), 1000, TRUE), 2)
Missing value(s) are 257 .
Valid case(s) are 743 .
Total case(s) are 1000 .
abc
Freq 250.00 218.00 275.00
Pct33.65% 29.34% 37.01%
> freq(iris$Species)
Missing value(s) are 0 .
Valid case(s) are 150 .
Total case(s) are 150 .
setosa versicolorvirginica
Freq 50.000 50.000 50.000
Pct 33.333% 33.333% 33.333%
> freq(iris$Species, 0)
Missing value(s) are 0 .
Valid case(s) are 150 .
Total case(s) are 150 .
setosa versicolorvirginica
Freq 50 50 50
Pct 33% 33% 33%
HTH,
Marc Schwartz
__
R-help@stat.math.ethz.ch mailing list
https://stat.ethz.ch/mailman/listinfo/r-help
PLEASE do read the posting guide! http://www.R-project.org/posting-guide.html