search for: dichotomous

Hello R-User! I appologise in advance if this should also go into statistics but I am presently puzzled. I have a data.frame (about 300 rows and about 80 variables) and my variables are dichotomous factors, continuous (numerical) and ordered factors. I would like to calculate the linear correlation between every pair of my variables, because I would like to perform a logistic regression (glm()) without the correlation between variables. I thought I could use for the continous (numerical) an...

I am trying to run a GLMM on some binomial data. My fixed factors include 2 dichotomous variables, day, and distance. When I run the model: modelA<-glmmPQL(Leaving~Trial*Day*Dist,random=~1|Indiv,family="binomial") I get the error: iteration 1 Error in MEEM(object, conLin, control$niterEM) : Singularity in backsolve at level 0, block 1 >From looking at previ...

Hi everyone, Have sample of items for each one, a set of 20 dichotomous (absent-present) variables are expressed. I'm trying to understand how to explore the co-occurence of each variable. Read some papers concerning smallest space analysis, but it does not seems implemented in any R package (and my protamming skills are =0). Non metric MDS gives error messages, pr...

Newsgroup members, Does anyone have a clever way to simulate a correlation matrix such that each column contains dichotomous variables (0,1) and where each column has different prevalence rates. For instance, I would like to simulate the following correlation matrix: > CORMAT[1:4,1:4] PUREPT PTCUT2 PHQCUT2T ALCCUTT2 PUREPT 1.0000000 0.5141552 0.1913139 0.1917923 PTCUT2 0.5141552 1.0000000 0.2913...

Dear community I'm pretty new to R and I'm trying to do a Point-baserial correlation for a nominal dichotomous variable with a interval scaled variable. It works fine, but the output just shows me the correlation and nothing else (p-Value would be important). I tried it with the following codes: - biseral.cor() - cor.biseral() - I also tried a polyserial() I've found on this question page: http://s...

Hello, I'd like to know if there is a way to fit a Structural equation model with dichotomous indicators (ex: problem with a phone solved/ or not) having effects on a ordinal variable. How I do that using R? Do you have an example with the code in R that you can send to me? Thanks a lot! Renata Estrella UFRJ, Brasil, Rio de Janeiro Renata Leite Estrella Assistente de Pesqui...

hi, How can I compute a reliability score of a scale consisting only of dichotomous items? thanks for any help!

Newsgroup members, I appreciate the help on this topic. David Duffy provided a solution (below) that was quite helpful, and came close to what I needed. It did a great job creating two vectors of dichotomous variables with a known correlation (what I referred to as a phi-coefficient). My situation is a bit more complicated and I'm not sure it is easily solved. The problem is that I must assume one of the vectors is constant and generate one or more vectors that covary with the constant vector. I...

I'd like to estimate an ordinal logistic regression with a random effect for a grouping variable. I do not find a pre-packaged algorithm for this. I've found methods glmmML (package: glmmML) and lmer (package: lme4) both work fine with dichotomous dependent variables. I'd like a model similar to polr (package: MASS) or lrm (package: Design) that allows random effects. I was thinking there might be a trick that might allow me to use a program written for a dichotomous dependent variable with a mixed effect to estimate such a model. The...

df$LAP <- with(df, ifelse(G=='male', (WC-65)*TG, (WC-58)*TG)) That will do both calculations and merge the two vectors appropriately. It will use extra memory, but it should be much faster than a 'for' loop. Regards, Jorgen Harmse. ------------------------------ Message: 8 Date: Fri, 3 Nov 2023 11:10:49 +1030 From: "Md. Kamruzzaman" <mkzaman.m at gmail.com>

Hello, I am running polychoric correlations on a dataset composed of 12 ordinal and binary variables (N =384), using the polycor package. One of the association (between 2 dichotomous variables) is very high using the 2-step estimate (0.933 when polychoric run only between the two variables; but 0.801 when polychoric run on the 12 variables). The same correlation run with ML estimate returns a singularity message. First, I would like to know why the estimations between only the...

Just a minor point in the suggested solution: df$LAP <- with(df, ifelse(G=='male', (WC-65)*TG, (WC-58)*TG)) since WC and TG are not conditional, would this be a slight improvement? df$LAP <- with(df, TG*(WC - ifelse(G=='male', 65, 58))) -----Original Message----- From: R-help <r-help-bounces at r-project.org> On Behalf Of Jorgen Harmse via R-help Sent: Friday,

...am planning. All I want to do is to fit a predictive model to this data and to graphically represent the 'best fit'. Any help will be greatly appreciated. Thanks in advance, Andreas -- View this message in context: http://r.789695.n4.nabble.com/predictive-logistic-model-cell-biology-non-dichotomous-data-tp3595812p3595812.html Sent from the R help mailing list archive at Nabble.com.

Greetings- I'm working with data that are multilevel in nature and have a dichotomous outcome variable (presence or absence of an attribute). As far as I can tell from reading archives of the R and S lists, as well as Pinheiro and Bates and Venables and Ripley, - nlme does not have the facility to do what amounts to a mixed-effects logistic regression. - The canonical alternative i...

Yes, that will halve the number of multiplications. If you?re looking for such optimisations then you can also consider ifelse(G=='male', 65L, 58L). That will definitely use less time & memory if WC is integer, but the trade-offs are more complicated if WC is floating point. Regards, Jorgen Harmse. From: avi.e.gross at gmail.com <avi.e.gross at gmail.com> Date: Friday,

..... I am trying to plot boxplots with side-by-side option.. I tried some of the posted suggestions and could not make it work due to unequal sizes of categories... e.g. weekly measured water depth values are categorized into 5 levels based on their values such measurement is again categorized into dichotomous levels - based on the result of a test I would like generate boxplot of water depth, with side-by-side display at dichotomous levels (0,1), for each category (1-5). However, dichotomous levels are not always available in each category.....e.g. category 1 may have all categorized as 0s, category 5...

Is it possible to fit a structural equation model with link functions in R? I am trying to build a logistic-regression-like model in sem, because incorporating the dichotomous variables linearly seems inappropriate. Mplus can do something similar by specifying a 'link' parameter, but I would like to be able to do it in R, ofcourse. I have explored the 'sem' package from John Fox, but it does not seem to be able to fit non-linear relations. Is there some...

Dear R users: I have a continuous outcome variable and four predictors, two continuous and two dichotomous. i would like to use the lattice plot to create scatter plots for the continuous predictors and boxplots for the dichotomous predictors. with 4 continuous variables, this is what i have been doing: trial = rbind ( cbind ( cimt$ant.mean, cimt$age, 1 ), cbind ( cimt$ant.mean, cimt$sbp, 2 ), cbind (...

...o this. > > I'm wanting to create a table in gtsummary that displays entries for "Yes" and "No" as fractions rather than a relative to the total (e.g. {n} / {N} in the documentation). > > I realize I can split them into rows with something like: > type = all_dichotomous() ~ "categorical" > > But I would honestly much rather display the variable as something like: > > Condition? (Y/N) count "Yes"/ count "No" > > > Is this possible to do in gtsummary? > > Thank you in advance. > The i...

Avi: Thank you for checking. I think the optimization is limited. If test is all TRUE or all FALSE then at most one vector is evaluated. Anything beyond that would be very complicated. (Inspect the two expressions and verify that both specify elementwise computations. Then use indexing to shrink the input properly. Take into account all recycling rules for binary operations.) > ifelse(0:1,