Error: Aesthetics must be either length 1 or the same as the data (1020)

[Marjory]

Hey, I'm fairly new to R and came across an error I cannot seem to solve. My data has a fair amount of missing values, so when constructing a new column of the residuals using:

linear <- lm(LN_ratio_CD3_Blast ~ ID, data = CMI_LNratioCD3)

CMI_LNratioCD3$residual <- residuals(linear)

I get the following error: Error in $<-.data.frame(*tmp*, residual, value = c(317 = -1.27035435466668, : replacement has 244 rows, data has 1020

I've tried adding na.exclude, na.omit, na.pass, na.fill etc. but I seem to be doing something wrong...

In order to surpass this problem I've tried to not create the residual column, and just plot this in the qq plot directly:

residual <- residuals(linear)

g5 <-ggplot(data = CMI_LNratioCD3,
aes(sample = residual))+
geom_qq() +
geom_qq_line(colour = "red") +
labs(title = "Quantile plot of residuals")

However, this will result in the following error: Error: Aesthetics must be either length 1 or the same as the data (1020): sample

Can someone please tell me what might be going on? And how I can solve this?

Kind regards,

Marjory

[mara]

Could you please turn this into a self-contained reprex (short for reproducible example)? It will help us help you if we can be sure we're all working with/looking at the same stuff.

install.packages("reprex")

If you've never heard of a reprex before, you might want to start by reading the tidyverse.org help page. The reprex dos and don'ts are also useful.

There's also a nice FAQ on how to do a minimal reprex for beginners, below:

FAQ: How to do a minimal reproducible example ( reprex ) for beginners

A minimal reproducible example consists of the following items: A minimal dataset, necessary to reproduce the issue The minimal runnable code necessary to reproduce the issue, which can be run on the given dataset, and including the necessary information on the used packages. Let's quickly go over each one of these with examples: Minimal Dataset (Sample Data) You need to provide a data frame that is small enough to be (reasonably) pasted on a post, but big enough to reproduce your issue. Let's say, as an example, that you are working with the iris data frame head(iris) #> Sepal.Length Sepal.Width Petal.Length Petal.Width Species #> 1 5.1 3.5 1.4 0.…

What to do if you run into clipboard problems

If you run into problems with access to your clipboard, you can specify an outfile for the reprex, and then copy and paste the contents into the forum.

reprex::reprex(input = "fruits_stringdist.R", outfile = "fruits_stringdist.md")

For pointers specific to the community site, check out the reprex FAQ.

[Marjory]

After running repress, I got a more detailed error and figured it out by myself. In short:
I've should have included the 'elimination' of the missingdata(NA) in the lm function when calculating the residuals with na.action=na.exclude

linear <- lm(LN_ratio_CD3_Blast ~ ID, data = CMI_LNratioCD3,na.action=na.exclude)
> CMI_LNratioCD3$residual <- residuals(linear)