Видео

Thank you very much Muhammad for such a nice feedback! Sure, in the beginning we'll all had difficulties to switch to R. I came from Matlab and NCSS to R. And also needed to box myself through the error messages. The good news is - the error messages are finite. The are only a few (20-50) error messages, you quickly learn how to deal with. After it error message will become a help. Levels are easy, you can determine the order yourself: library(dplyr) library(forcats) # install the packages, if they don't load df <- data.frame( category = c("B", "A", "C") ) # Reorder levels df <- df %>% mutate(category = fct_relevel(category, "A", "C", "B")) # Print the reordered factor levels levels(df$category)

@@yuzaR-Data-Science Thank you very much. Will also be looking forward to more video.

thanks

Glad it helped! I don't know about the methods switch after update. But I know that the author of the package is responsive and you can ask him via github or twitter or any other way. He'll give much better answer.

@yuzaR-Data-Science i will try that, thank you

@@thomaswiggersmller8983 👍

no, you can sort them out via p-values, e.g. <0.2, but the variable importance should only be asked from the multivariable models. folks in the medicine have little idea about stats, thus, take their methods with a grain of salt and consult a statistician ;)

@@yuzaR-Data-Science thanks for replying! Just to clarify would you put all of the available predictors in a multivariate model and then based on p-values <0.2 adjust the model accordingly? (like backwards selection process?)

hey, unfortunately the robustbase is not supported yet. but you can make a request on the github of the author. most other models work well

Thank you! ❣

@@FabianMueller-n6g you are welcome! :)

Glad you enjoyed it, Elias! I am working on roc curve and optimal cutpoint video right now. Hope it will deliver the things you are interested in. Stay tuned. Kind regards from holidays in Australia

ok

Hi Anthony, thanks, and sure, feel free to use my content and direct others to my videos. I would actually appreciate that! :) The janitor video is useful, but kind of old, so, you might have a look at the newer ones, I think they are even more useful. Let me know when your content is ready, I would also love to see it! Cheers

You are very welcome 🙏 thank you for watching

Glad it was helpful! Sure, there are several functions for ROC curves in R. Several packages provide good results, but I like two of them more then the rest: Epi::ROC(form = survived ~ predicted_glm, data = d, plot = "ROC", grid = F, MX = T, MI = F, lwd = 3) cutpointr() - I am workind on a whole video about this one, it's just amazing

Unfortunately not. Only two older theoretical videos on survival, but they low quality and no programming. Plan to do the similar one in the future. So, please, stay tuned.

@ looking forward to that. Thanks for this great vid nonetheless!

welcome!

Thanks 🙏

Hi Nike, thanks for the positive feedback. And I am interested in your 'chisquare' package. Unfortunately I did not find much info online on it. I have actually already made one video on chi-squared test. If you have seen this one, what does your package does better and differently? If you send me the code for what your package can do and explanations why it is useful and why it is better than usuals chi-square function or ggbarstats, I would love to make a video on your package!

@@yuzaR-Data-Science Hello, and thanks for your reply. The package is on CRAN, and it's currently in its version 1.1.1 (it started from vers 0.1 in 2022). In few words, the package is meant to provide a one-stop shop for chi-square analysis of cross-tabs, and provides a number of facilities that are not coherently integrated in existing packages (to the best of my knowledge). For example, it provides (in just one simple line of code), different types of chi-sq residuals (with adjustements for multiple comparisons, and color coded for easy visual interpretation) and a extensive suite of association coefficients (for both 2x2 and larger tables), some of which not currently implemented elsewhere (maximum-corrected version of the phi and Cramer's V coeff, corrected version of Goodman-Kruskal's lambda, both asymmetric and symmetric). Also, it provides different versions of the chi-sq test itself, like the N-1-corrected version, which (again) is not currently provided elsewhere. As for post-hoc-analysis, it provides measures not currently available elsewhere, like the so-called Quetelet index and the IJ association factor. Further, it computes independent odds ratios for tables larger than 2x2, while for 2xK tables it can optionally produce a plot of pair-wise odds ratios (plus confidence intervals). Also, it provides suggestions as to a 'viable' chi-sq test given the input table characteristics. Effect size verbal articulation for relevant association coefficients (both chi-square-based and marginal-free) are also reported. Finally, all the outputs are nicely formatted via the 'gt' table package. I think that should be almost pretty much all. Everything can be obtained by just running: chisquare(mytable). Cheers.

hey, your package is impressive, I found the visualization of odds ratios good. I have two questions: - first, do you have more info, like article or so on post hoc pairwise tests with all the significance, like when we have a table 4x4 or 3x5, so that all categories (percentages) are checked automatically. till now I use a pairwise_fishers_test() function which is cool, but an extra code. It would be amazing when we could just use your function and get all we need - ORs plot with significances and all the pairwise 2x2 tests from bigger contingency table in some form of a table. - second, may be more important: I could not get chisquare() function work with a simple table() function: > chisquare(table(mtcars$cyl, mtcars$am) ) Error in `gt::tab_style()`: ! Failed to style the body of the table. Caused by error in `cells_body()`: ! Can't select columns that don't exist. ✖ Column `0` doesn't exist. Run `rlang::last_trace()` to see where the error occurred. so, when this can be allowed and we could do bigger tables, like this one: chisquare(table(ISLR::Wage$jobclass, ISLR::Wage$education) ), this could be awesome!

@@yuzaR-Data-Science Hello. Thanks for taking the time to check that and for replying. I do not want to hijack your comments section here. If you want to contact me on the email you find in the package documentation, I will more than happy to discuss things further. Looking forward. Cheers.

hey mate, no worries, you don't hijack the comments section! :) I am actually glad to read and answer the comments. the next weeks I'll be on holidays, but we can talk about your package next month. generally, as I said before, I would love to be able to apply your chisquare function to a simple cross table, like that "chisquare(table(mtcars$cyl, mtcars$am) )". do you think it's possible?

You are very welcome 🙏

Thank you! Cheers!

hi, yes, this is similar to subgroup modelling, or stratification, but I don't know what you mean with moderator

You are very welcome!

Thanks for your positive feedback and for the insight with glht function! Emmeans function is also used in most of my favorite packages, like gtsummary for example, where the contrasts can be extracted and p-values can be corrected with any method

@@yuzaR-Data-Science Thank you for one of the best channels on RUclips on the topic of data science using R. By the way, I had a Jewish grandpa from Belarus and Latvia (he was born in Latvia, but moved to Belarus, before arriving at South Africa). Your accent is quite similar to his, but with an extra German sound. Is there any chance that you are from Belarus?

@@itamar.j.rachailovich hey man, well recognized! respect! :) I was born and grew up in Belarus, but live in Germany since > 20 years. Where are you located? And thanks for your generous and nice feedback, Itamar!

@@yuzaR-Data-Science I am located in Israel, as my Belarusian grandfather was Jewish (note please that I have nothing to do with Israeli politics, military and operations, I have never served in their army, and I am totally for one country for Palestinians and Jews with equal rights for all, and right of return for the millions of Palestinians that were expelled) I am originally from South Africa, my grandpa migrated to that country. The other half of my family are South African Boers - descendants of Dutch and Germans ( 23andme l am 45% German from souther Jutland peninsula (near the Dane border). I have degrees in psychology and biology, and currently at the Weizmann institute of science, brain science department. Just started my journey in computing, I knew nothing about it, and I am fascinated by it and by the advanced statistics required to the extent that I would like to change the direction of my career to a more statistics/programming oriented biology / brain sciences

Such a rich history you have, Belarus, South Africa and Germany + Dutch :) cool! Yeah, I am also not political and wish the peace in the world! So, we are in the same boat there ;) I am not Jewish myself, but some Jewish relatives in Israel (Haifa), Belarus and Germany. Wanted to visit them before corona, then corona came and now those conflicts. So, I guess I'll wait a few years before going to Israel. I studied Animal Science and Marine Biology originally, but did my PhD in computational biology. That was the moment where I fall in love with scripting and data science... I don't say programming, because I usually use preprogrammed functions, but don't write many of them myself. Love scripting though. So, I could imagine you will enjoy that too, and for sure programming in R would be super beneficial for your carrier! Especially in bioinformatics. :)

thanks a lot !!!

That’s a brilliant insight! 😀 thanks mate 🙏 I’ll use it to tell my colleagues who take outliers either too seriously or too not-seriously 👍

thanks :) of coarse ... that's why I have few further videos on logistic regression ... one on only categorical predictor, and another - multivariable logistic regression ;) hope you enjoy those too

the best way to reach me is to comment here. king regards!

@@yuzaR-Data-Science sir i want to show you some figures and want to know how can i made them in R

@@nihargupta7358 I have created several videos on ggplot graphs. so, have a look at those. or ask chatgpt how to make this plots in R. you can copy-paste pics into chatgpt

first, xlab & ylab, like in ggplot, second, p < 0.001 not possible. only scientific notation.

Thanks a ton, Marco 🙏 I’ll do my best to keep the content going 😉 hope you like other videos too. Kind regards