Hi folks i tried deconvoluting with PSEA package

first, according to vignette; i defined the markers before this; none of the probesets contain any NA values;

      neuronref<-marker(neuronprobeset,expr=counts,targetMean = 100)
    mgref<-marker(mgprobeset,expr=counts,targetMean = 100)

 astrocyteref<-marker(astrprobeset,expr=counts,targetMean = 100)

Defined the groups (corresponding to my expression matrix as well:

    x<-c(0,1)
groups<-rep(x,each=4)

However when i defined linear model

model.group0 <- lm(counts["Apoe",] ~ neuronref + astrocyteref +
           mgref,subset=which(groups==0))

I realised that mgref is NA and

summary (model.group0)

only shows Estimates of neuronref, astrocyteref. The rest; whether it is p-values or degree of freedom are NA. this problem does not happen if i remove subset=which(groups==0);

Would appreciate some advice, thank you.