ludvigla / semla

Other
47 stars 6 forks source link

neighborhood mapping error #14

Closed chrkuo closed 10 months ago

chrkuo commented 11 months ago

Hi,

@ludvigla I transitioned to semla after hitting a wall with proceeding with neighborhood analysis in STUTility but the same error is popping up when I am trying to map neighborhoods. It consistently keeps saying the rows don't match although I see my meta-data with a lot of NAs and those that were annotated to be neighbors. I have posted this issue with dplyr and they think its something wrong inherently with the package of semla and not dplyr. Can I please get some assistance.

Thank you

ludvigla commented 11 months ago

Hi,

I'll try to help you out but I will need more information to go on. If you have NA values in the meta data that could be the issue.

Would you be able to share the object? That would make it easier for me to figure out what's wrong.

/Ludvig

chrkuo commented 11 months ago

@ludvigla absolutely what is your email. I will send it right now

I also tried it with the example (Mouse Brain) given in the vignette and it's the same output so I'm wondering if there is a bug in the regionneighbor code.

I thought the region-neighbor code is supposed to generate NA for any spots that are not neighbors. I had a bunch of NAs along with few nbs_X.

An example of the output datatable:

[1] NA NA NA NA NA NA NA NA NA NA NA NA NA NA [15] NA NA NA NA NA NA NA NA NA NA NA NA NA NA [29] NA NA NA NA NA NA NA NA NA NA NA NA NA NA [43] NA NA NA NA NA NA NA NA NA NA NA NA NA NA [57] NA NA NA NA NA NA NA NA NA NA NA NA NA NA [71] NA NA NA NA NA NA NA NA NA NA NA NA NA NA [85] NA NA NA NA NA NA NA NA NA NA NA NA NA NA [99] NA NA NA NA NA NA NA NA NA NA NA NA NA NA [113] NA NA NA NA NA NA NA NA NA NA NA NA NA NA [127] NA NA NA NA NA NA NA NA NA NA NA NA NA NA [141] NA NA NA NA NA NA NA NA NA NA NA NA NA NA [155] NA NA NA NA NA NA NA NA NA NA NA NA NA NA [169] NA NA NA NA NA NA NA NA NA NA NA NA NA NA [183] NA NA NA NA NA NA NA NA NA NA NA NA NA NA [197] NA NA NA NA NA NA NA NA NA NA NA NA NA NA [211] NA NA NA NA NA NA NA NA NA NA NA NA NA NA [225] NA NA NA NA NA NA NA NA NA NA NA NA NA NA [239] NA NA NA NA NA NA NA NA NA NA NA NA NA NA [253] NA NA NA NA NA NA NA NA NA NA NA NA NA NA [267] NA NA NA NA NA NA NA NA NA NA NA NA NA NA [281] NA NA NA NA NA NA NA NA NA NA NA NA NA NA [295] NA NA NA NA NA NA NA NA NA NA NA NA NA NA [309] NA NA NA NA NA NA NA NA NA NA NA NA NA NA [323] NA NA NA NA NA NA NA NA NA NA NA NA NA NA [337] NA NA NA NA NA NA NA NA NA NA NA NA NA NA [351] NA NA NA NA NA NA NA NA NA NA NA NA NA NA [365] NA NA NA NA NA NA NA NA NA NA NA NA NA NA [379] NA NA NA NA NA NA NA NA NA NA NA NA NA NA [393] NA NA NA NA NA NA NA NA NA NA NA NA NA NA [407] NA NA NA NA NA NA NA NA NA NA NA NA NA NA [421] NA NA NA NA NA NA NA NA NA NA NA NA NA NA [435] NA NA NA NA NA NA NA NA NA NA NA NA NA NA [449] NA NA NA NA NA NA NA NA NA NA NA NA NA NA [463] NA NA NA NA NA NA NA NA NA NA NA NA NA NA [477] NA NA NA NA NA NA NA NA NA NA NA NA NA NA [491] NA NA NA NA NA NA NA NA NA NA NA NA NA NA [505] NA NA NA NA NA NA NA NA NA NA NA NA NA NA [519] NA NA NA NA NA NA NA NA NA NA NA NA NA NA [533] NA NA NA NA NA NA NA NA NA NA NA NA NA NA [547] NA NA NA NA NA NA NA NA NA NA NA NA NA NA [561] NA NA NA NA NA NA NA NA NA NA NA NA NA NA [575] NA NA NA NA NA NA NA NA NA NA NA NA NA NA [589] NA NA NA NA NA NA NA NA NA NA NA NA NA NA [603] NA NA NA NA NA NA NA NA NA NA NA NA NA NA [617] NA NA NA NA NA NA NA NA NA NA NA NA NA NA [631] NA NA NA NA NA NA NA NA NA NA NA NA NA NA [645] NA NA NA NA NA NA NA NA NA NA NA NA NA NA [659] NA NA NA NA NA NA NA NA NA NA NA NA NA NA [673] NA NA NA NA NA NA NA NA NA NA NA NA NA NA [687] NA NA NA NA NA NA NA NA NA NA NA NA NA NA [701] NA NA NA NA NA NA NA NA NA NA NA NA NA NA [715] NA NA NA NA NA NA "nbs_6" NA "nbs_6" NA "nbs_6" NA NA NA [729] NA NA NA NA NA NA NA NA NA NA NA NA NA NA [743] NA NA NA NA NA NA NA NA NA NA NA NA NA NA [757] NA NA NA NA NA NA NA NA NA NA NA NA NA NA [771] NA NA NA NA NA NA NA NA NA NA NA NA NA "nbs_6" [785] "nbs_6" "6" "6" "6" "6" "6" "nbs_6" "6" NA NA NA NA NA NA [799] NA NA NA NA NA NA NA NA NA NA NA NA NA NA [813] NA NA NA NA NA NA NA NA NA NA NA NA NA NA [827] NA NA NA NA NA NA NA NA NA NA NA NA NA NA [841] NA NA NA NA "nbs_6" "nbs_6" "nbs_6" "nbs_6" "6" "6" NA NA NA NA [855] NA NA "6" "6" "nbs_6" NA NA NA NA NA NA NA NA NA [869] NA NA NA NA NA NA NA NA NA NA NA NA NA NA [883] NA NA NA NA NA NA NA NA NA NA NA NA NA NA [897] NA NA NA NA NA NA NA NA NA NA NA NA NA NA [911] "nbs_6" "nbs_6" "6" "6" "6" NA NA NA NA NA NA NA NA NA [925] "6" "6" "nbs_6" NA NA NA NA NA NA NA NA NA NA NA [939] NA NA NA NA NA NA NA NA NA NA NA NA NA NA [953] NA NA NA NA NA NA NA NA NA NA NA NA NA NA [967] NA NA NA NA NA NA NA NA NA NA NA NA NA NA [981] NA NA "nbs_6" "nbs_6" "6" "6" "nbs_6" "6" "6" NA NA NA NA NA [995] NA NA NA NA NA NA

chrkuo commented 11 months ago

@ludvigla any updates. is the regionneighbour code supposed to generate NA values?

ludvigla commented 11 months ago

Here's my email: ludvig.larsson@scilifelab.se.

I think I will need to take a closer look on the data in order to figure out what's going on. Could you aslo share the code that you ran?

/Ludvig

chrkuo commented 10 months ago

@ludvigla any updates.

Thank you so much

ludvigla commented 10 months ago

For anyone who might be interested.

While we couldn't precisely identify which specific R package was responsible for the conflict, the code executed smoothly after restarting R.

Closing this issue as resolved.