Как использовать функцию summary.df в пакете psych в R с тремя интересующими переменными?

Question

Я бы хотел найти показатели центральной тенденции для моего набора данных с помощью функции desribeBy в пакете psych.Кажется, он не работает, и говорит: «аргумент не числовой или логический: возвращение NA».Для некоторой предыстории я смотрю на связь между использованием Facebook (измеряется как «Пользователь» или «Не пользователь») и Политическим интересом (измеряется по шкале 1-7).

setwd("~/Desktop")

df <- read.csv("2016 ANES.csv")

# DV: V162174a (POST: Days in past week discussed politics)  
# IV: V162370 (POST: Facebook account used recently)
# Control: V161270 (PRE: Highest level of Education)

# Cleaning up my Dependent Variable; Political Involvement

table(df$V162174a) 

# I ignored -1, -6, -7, -9 and -8 (non-responses) when recoding my variable; 'Political interest' 

df$Political_interest[df$V162174a=="0. Zero days"] <- 0
df$Political_interest[df$V162174a=="1. One day"] <- 1
df$Political_interest[df$V162174a=="2. Two days"] <- 2
df$Political_interest[df$V162174a=="3. Three days"] <- 3
df$Political_interest[df$V162174a=="4. Four days"] <- 4
df$Political_interest[df$V162174a=="5. Five days"] <- 5
df$Political_interest[df$V162174a=="6. Six days"] <- 6
df$Political_interest[df$V162174a=="7. Seven days"] <- 7

table(df$Political_interest)

# Cleaning up my Independent variable; Facebook Usage

table(df$V162370)

df$Facebook_usage[df$V162370=="1. Yes have a Facebook account I have used in the past month"] <- "User"
df$Facebook_usage[df$V162370=="2. Have a Facebook account but have not used it in past month"] <- "Non_user"
df$Facebook_usage[df$V162370=="3. No, do not have a Facebook account"] <- "Non_user"

table(df$Facebook_usage)

# Comparing facebook usage and Political engagement

as.numeric(df$Facebook_usage)
na.omit(df$Facebook_usage)

table(df$Facebook_usage, df$Political_interest)

# Creating a grouped bar plot 

Political_interest <- table(df$Political_interest)
Political_interest

Facebook_usage <- table(df$User, df$Non_user)
User
Non_user

install.packages("plotly")
library(plotly)

Political_interest <- c("0", "1", "2", "3", "4", "5", "6", "7") 
User <- c(88, 195, 266, 196, 171, 181, 79, 514)
Non_user <- c(72, 169, 190, 156, 119, 123,  51, 308)
df <- data.frame(Political_interest, User, Non_user)

p <- plot_ly(df, x = ~Political_interest, y = ~User, type = 'bar', name = 'User') %>%
  add_trace(y = ~Non_user, name = 'Non User') %>%
  layout(yaxis = list(title = 'Respondents'), barmode = 'group')
p

# Looking at measures of central tendancy 

install.packages("psych")
library(psych)

describeBy(df$Political_interest, df$User, df$Non_user)

mean(df$Political_interest)