Skip to content
New issue

Have a question about this project? Sign up for a free GitHub account to open an issue and contact its maintainers and the community.

By clicking “Sign up for GitHub”, you agree to our terms of service and privacy statement. We’ll occasionally send you account related emails.

Already on GitHub? Sign in to your account

Error when no compensation control for a channel #35

Open
wudustan opened this issue Aug 4, 2021 · 2 comments
Open

Error when no compensation control for a channel #35

wudustan opened this issue Aug 4, 2021 · 2 comments

Comments

@wudustan
Copy link

wudustan commented Aug 4, 2021

I am trying to use spillover() but not all the channels in the experiment have an appropriate stain - one of our lines expresses Katushka protein which we use in the flow cytometry experiment, but we don't have an equivalent stain for beads.

I unwound the package function for spillover and just manually ran the steps. At the bottom I also added one more row with 0s for all other channel and 1 for Katushka to just try and get across the problem but accepting it will affect the resulting data. I just wanted to double check how this will affect the rest of my channels and wanted to get the take of people who have worked on the package and their data

Here I include the matrix heatmap of the computed spillover matrix. Stained control samples are rows with Katushka added manually, columns are the $P channels detected by flowCore/flowStats.

Screenshot 2021-08-04 at 10 35 50

@gfinak
Copy link
Member

gfinak commented Aug 6, 2021

What's the question?
Are you computing the spillover matrix from controls or are you just extracting it from the file?
One problem is your off-diagonal elements are 1.0, and your diagonal elements are < 1.
What version of flowCore?
Any chance you can share a reproducible example that generates this matrix?

@wudustan
Copy link
Author

Hi,

Thanks for the response. Unfortunately, I can't provide my data.

I have two questions:

  1. What is the best course of action when one of the $P channels doesn't have an associated stained control?
    -- in my example, I simply added an extra row to the matrix and set every value to 0 except the column associated with that channel. But I'm conscious that's a horrible hack.

  2. For cases where as you correctly point out off-diagonal elements are 1.0, and diagonal elements are < 1 despite knowing that the stain is the correct stain for that channel, do I manually amend this?
    -- I have made sure to double check this isn't a channel/stain mix up. For some reason, some of the unstained controls have a higher signal in channels out of their normal target.

flowCore: 2.4.0
R: 4.1.0

Sign up for free to join this conversation on GitHub. Already have an account? Sign in to comment
Labels
None yet
Projects
None yet
Development

No branches or pull requests

2 participants