0

I have 2 categorical variables income level and Temporary visa status and the count for each combination.

All i need is a crosstab for creating a proportion crosstab bar chart to get the proportion of different temporary visa categories within an income level 

library(readxl)
Crosstab_Temporary_visas_income <- read_excel("C:/Users/axelp/Documents/RMIT/Semester 2/Data Visualisation/Assignment 3/Crosstab Temporary visas income.xls")

str(Crosstab_Temporary_visas_income)

margin.table(Crosstab_Temporary_visas_income,1) #Row marginals

Error in margin.table(Crosstab_Temporary_visas_income, 1) : 
  'x' is not an array

> str(Crosstab_Temporary_visas_income)
Classes ‘tbl_df’, ‘tbl’ and 'data.frame':   9 obs. of  6 variables:
 $ Income                  : chr  "Negative / nil income" "$1– 299" "$30 - 649" "$650– 999" ...
 $ Temporary Work (Skilled): num  405 2364 6496 19248 41595 ...
 $ Student                 : num  2169 33846 104569 27140 6737 ...
 $ New Zealand Citizen     : num  2446 16045 51337 104133 98986 ...
 $ Working Holiday Maker   : num  515 3670 18119 24476 7869 ...
 $ Other Temporary visa    : num  887 5325 24234 31975 16269 ...

structure(list(...1 = c("0", "$1– 299", "$30 - 649", "$650– 999"
), `Temporary Work (Skilled)` = c(405, 2364, 6496, 19248), Student = c(2169, 
33846, 104569, 27140), `New Zealand Citizen` = c(2446, 16045, 
51337, 104133), `Working Holiday Maker` = c(515, 3670, 18119, 
24476), `Other Temporary visa` = c(887, 5325, 24234, 31975)), row.names = c(NA, 
-4L), class = c("tbl_df", "tbl", "data.frame"))

I used the table function on the imported csv to create a crosstab but all i get is more than 6000 matrix slices

axel_p
  • 51
  • 9
  • 2
    When asking for help, you should include a simple [reproducible example](https://stackoverflow.com/questions/5963269/how-to-make-a-great-r-reproducible-example) with sample input and desired output that can be used to test and verify possible solutions. – MrFlick Nov 05 '19 at 02:20
  • Hi Mr Flick i have added some data to my OP. Kindly let me know if you need anything else . thanks – axel_p Nov 05 '19 at 02:39

1 Answers1

0

There was a problem when you read your data in because the row names were read as a column labeled "...1". Usually R will recognize a row names if there is one fewer column name than the number of columns. Nothing will work until you fix that.

library(tidyverse)
CTVI <- structure(list(...1 = c("0", "$1– 299", "$30 - 649", "$650– 999"), 
`Temporary Work (Skilled)` = c(405, 2364, 6496, 19248), Student = c(2169, 
33846, 104569, 27140), `New Zealand Citizen` = c(2446, 16045, 
51337, 104133), `Working Holiday Maker` = c(515, 3670, 18119, 
24476), `Other Temporary visa` = c(887, 5325, 24234, 31975)),  
row.names = c(NA, -4L), class = c("tbl_df", "tbl", "data.frame"))

Now we need to delete the first column, use it for the row names, and convert the tibble to a matrix since some table functions such as addmargins and margin.table do not accept tibbles:

CTVI.mat <- as.matrix(CTVI[, -1])
rownames(CTVI.mat) <- unlist(CTVI[, 1])
CTVI.mat <- CTVI.mat[, -1]
names(dimnames(CTVI.mat)) <- c("Income", "Visa")

Now we can compute margins or proportions:

margin.table(CTVI.mat, 1) 
addmargins(as.matrix(CTVI.mat))
round(prop.table(as.matrix(CTVI.mat), 1), 3)
dcarlson
  • 10,936
  • 2
  • 15
  • 18
  • Thanks for looking in @dcarlson. I read in CTVI as advised in the first step then got an error "Setting row names on a tibble is deprecated.Error in `.rowNamesDF<-`(x, value = value) : invalid 'row.names' length" so changed the order to CTVI.tbl <- CTVI[, -1] CTVI.mat <- as.matrix(CTVI.tbl) rownames(CTVI.mat) <- CTVI[, 1] names(dimnames(CTVI.mat)) <- c("Income", "Temporary Visa Class") now im getting an error in the 3rd line Error in dimnames(x) <- dn : length of 'dimnames' [1] not equal to array extent – axel_p Nov 05 '19 at 07:28
  • Try this. You didn't indicate that you had loaded the `tidyverse` package which was throwing the errors, but I should have guessed that you were using it since you were using `readxl`. Now it should work. – dcarlson Nov 05 '19 at 15:13