Data Structure
Networking
RDBMS
Operating System
Java
MS Excel
iOS
HTML
CSS
Android
Python
C Programming
C++
C#
MongoDB
MySQL
Javascript
PHP
Physics
Chemistry
Biology
Mathematics
English
Economics
Psychology
Social Studies
Fashion Studies
Legal Studies
- Selected Reading
- UPSC IAS Exams Notes
- Developer's Best Practices
- Questions and Answers
- Effective Resume Writing
- HR Interview Questions
- Computer Glossary
- Who is Who
How to convert MANOVA data frame for two-dependent variables into a count table in R?
MANOVA refers to multivariate analysis of variance, in this method we have more than one dependent variable and multiple independent variables. We want to compare each level of the independent variable combination for each of the dependent variables. To convert MANOVA data frame for two-dependent variables into a count table, we can use cast function of reshape package but we need to melt the data frame first so that the casting can be done appropriately.
Example
Consider the below data frame −
Gender<−sample(c("Male","Female"),20,replace=TRUE)
Class<−sample(c("I","II","III"),20,replace=TRUE)
Score<−sample(1:100,20)
Rating<−sample(1:10,20,replace=TRUE)
df1<−data.frame(Gender,Class,Score,Rating)
df1Output
Gender Class Score Rating 1 Male II 96 9 2 Male I 38 3 3 Female III 32 5 4 Male I 77 2 5 Male I 62 2 6 Female II 81 9 7 Male II 90 2 8 Female III 79 8 9 Male III 34 8 10 Male II 36 9 11 Male I 57 5 12 Male I 29 1 13 Female III 100 7 14 Female II 94 5 15 Male I 35 9 16 Female III 78 4 17 Female I 18 3 18 Female I 47 9 19 Female III 61 1 20 Male III 60 3
Loading reshape package −
library(reshape)
Melting df1 −
df1_melt<−melt(df1)
Using Gender, Class as id variables
Finding the counts based on Gender and Class −
cast(df1_melt,Gender~Class+variable)
Aggregation requires fun.aggregate: length used as default
Gender I_Score I_Rating II_Score II_Rating III_Score III_Rating 1 Female 2 2 2 2 5 5 2 Male 6 6 3 3 2 2
Let’s have a look at another example −
ID<<sample(c("Y1","Y2","Y3","Y4"),20,replace=TRUE)
Grade<<sample(LETTERS[1:3],20,replace=TRUE)
Sal<<sample(20000:50000,20)
Count<<sample(200:210,20,replace=TRUE)
df2<<data.frame(ID,Grade,Sal,Count)
df2Output
ID Grade Sal Count 1 Y3 B 28528 204 2 Y3 C 40854 207 3 Y3 A 31199 207 4 Y4 B 25338 207 5 Y3 B 30180 209 6 Y2 B 29921 209 7 Y4 C 46134 210 8 Y4 B 46829 205 9 Y3 B 42607 205 10 Y1 A 38174 202 11 Y2 A 41451 207 12 Y1 C 23912 200 13 Y4 B 44047 209 14 Y2 B 32236 200 15 Y2 A 24851 203 16 Y2 B 36341 207 17 Y3 B 37003 208 18 Y2 C 37285 207 19 Y3 B 45113 207 20 Y3 A 40034 203 df2_melt<−melt(df2)
Using ID, Grade as id variables
Finding counts based on ID and Grade −
cast(df2_melt,ID~Grade+variable)
Aggregation requires fun.aggregate: length used as default
ID A_Sal A_Count B_Sal B_Count C_Sal C_Count 1 Y1 1 1 0 0 1 1 2 Y2 2 2 3 3 1 1 3 Y3 2 2 5 5 1 1 4 Y4 0 0 3 3 1 1
144 Views
