Convert X to coefficient, and then use lm (Y ~ X + A). Or you can use dummyvars from caret package -
dummy_train<-dummyVars(" ~ .",data=<insert_data_name>) dummy_train<-data.frame(predict(dummy_train,newdata=<insert_the_same_data_name>))
You can run a regression on this.
source share