Can you please help confirm if my code is correct and makes sense

I am trying to predict the following (Violent~ year + state + law) to see if gun laws impact violent crime rates

*******P,Q,N, I am worried I am giving wrong variables ****** for hypothesis testing

#degrees of freedomp 

Hide

#F-stat and its p-valueMSR 

[1] 280.5283

************all code below use CTRL+F "hide" to find where my code starts************

rm(list=ls())library(lattice)data("Guns")library(AER)library(DAAG)library(dplyr)library(ggplot2)str(Guns)

'data.frame': 1173 obs. of 13 variables: $ year : Factor w/ 23 levels "1977","1978",..: 1 2 3 4 5 6 7 8 9 10 ... $ violent : num 414 419 413 448 470 ... $ murder : num 14.2 13.3 13.2 13.2 11.9 10.6 9.2 9.4 9.8 10.1 ... $ robbery : num 96.8 99.1 109.5 132.1 126.5 ... $ prisoners : int 83 94 144 141 149 183 215 243 256 267 ... $ afam : num 8.38 8.35 8.33 8.41 8.48 ... $ cauc : num 55.1 55.1 55.1 54.9 54.9 ... $ male : num 18.2 18 17.8 17.7 17.7 ... $ population: num 3.78 3.83 3.87 3.9 3.92 ... $ income : num 9563 9932 9877 9541 9548 ... $ density : num 0.0746 0.0756 0.0762 0.0768 0.0772 ... $ state : Factor w/ 51 levels "Alabama","Alaska",..: 1 1 1 1 1 1 1 1 1 1 ... $ law : Factor w/ 2 levels "no","yes": 1 1 1 1 1 1 1 1 1 1 ...

Hide

sum(is.na(Guns))

[1] 0

Hide

pairs(~ violent + state + year + law , data = Guns, main = "Scatterplot Matrix")

\f\f\fO O O 50 100 150 200 250 |Residuals) O O OC O O O O O O O O O O O O O O O 500 1000 1500 2000 2500 Fitted\f\f\f\f\f\f\f\f\f\f\f\f\f\f\f\f\f\f\f\f\f\f\f\f\f\f\f\f\f\f\f\f\f\f\f\f\f\f\f\f\f