Skip to content

Add markers and fixed intializations, fix an error#15

Draft
almdudleer wants to merge 4 commits intoalmdudleer/formattingfrom
almdudleer/markers_init
Draft

Add markers and fixed intializations, fix an error#15
almdudleer wants to merge 4 commits intoalmdudleer/formattingfrom
almdudleer/markers_init

Conversation

@almdudleer
Copy link
Collaborator

No description provided.

Error: "no match for ‘operator=’ (operand types are ‘arma::subview_row<double>’ and ‘<brace-enclosed initializer list>’)"
@almdudleer almdudleer requested a review from denklewer June 20, 2023 17:10
marker_means <- lapply(markers_list, function(x) {
sel <- x[x %in% rownames(self$new_points)]
if (sum(sel) == 0) stop("One or more cell types has zero markers present")
colMeans(self$new_points[sel, ])
Copy link
Collaborator Author

@almdudleer almdudleer Jun 21, 2023

Choose a reason for hiding this comment

The reason will be displayed to describe this comment to others. Learn more.

Essentially, markers initialization is this line only. For the set of marker genes for each cell type, we take the mean of the corresponding self$new_points column, and that is the X. Judging by the plots, it works fine, and proportions, generated from that seem to be good.

self$setInitParamsByX()
},

setInitParamsByX = function() {
Copy link
Collaborator Author

Choose a reason for hiding this comment

The reason will be displayed to describe this comment to others. Learn more.

However the basis generated from this function doesn't seem to be good in 10d:
image

Here's the math for it:
image

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

Successfully merging this pull request may close these issues.

1 participant