search for: depression

...on, just the warning message: "Could not compute QR decomposition of Hessian. Optimization probably did not converge. in: sem.default(ram = ram, S = S, N = N, param.names = pars, var.names = vars, " What does it mean? Maybe I am doing something wrong? I have 3 latent factors (Anxiety, Depression, and Fear) - each of them has 3 observed indicators (a1, a2, a3; d1, d2, d3, and f1, f2, f3) Below is my script in R: ! ANALYSIS OF ANXIETY, DEPRESSION, AND FEAR - LISREL P.31 ! Creating the ANXIETY, DEPRESSION, AND FEAR intercorrelation matrix (KM): KM<-matrix( c(1,.8,.83,.2,.21,.19,.18,.18,....

Hello, I finally run my measurement model in sem - successfully. Now, I am trying to print out the path diagram that is based on the results - but for some reason it's not working. Below is my script - but the problem is probably in my very last line: # ANALYSIS OF ANXIETY, DEPRESSION, AND FEAR - LISREL P.31 library(sem) # Creating the ANXIETY, DEPRESSION, AND FEAR intercorrelation matrix (KM) KM<-matrix( c(1,.8,.83,.2,.21,.19,.18,.18,.18, 0,1,.81,.22,.24,.18,.19,.19,.21, 0,0,1,.22,.19,.2,.2,.2,.22, 0,0,0,1,.84,.82,.22,.22,.21, 0,0,0,0,1,.84,.19,.18,.19, 0,0,0,0,0,1,.18,.18,....

...yone see why dffits might be failing? Is there a problem with the data? Consider: # Load data dep <- read.table("http://www.sci.usq.edu.au/staff/dunn/Datasets/Books/Hand/Hand-R/factor1-R.dat", header=TRUE) attach(dep) dep # Fit Poisson glm dep.glm2 <- glm( Counts ~ factor(Depression) + factor(SLE) + factor(Children) + factor(Depression):factor(SLE), family=poisson( link=log) ) # Compute dffits dffits( dep.glm2 ) This produces the output: 1 2 3 4 5 6 1.4207746 -0.1513808 NaN 0.9079142 -0.1032664 -1.08...

How should I define the :onsuccess option in my - link_to_remote - tag, to execute a Js function I wrote : :success => "toggleButton(this, /db_bfilter/i);" ''this'' is the current element defined in my view : <span id="db_bfilter0" class="depressed"><%= filter_link_helper "This week" %></span> <span

Dear all, I noticed the following bug in points.formula > library(DAAG) > data(roller) > fm <- lm(depression ~ weight, data=roller) > plot( depression ~ weight, data=roller, type="n") > abline(fm) > attach(roller) > points( depression~weight, subset=1:7) > points( depression~weight, subset=8:10, col="blue") Error in if (length(x) == l) x[s] else x : argument is of leng...

...variables > resp<- model.response(model.frame(soils.mod)) > # 1-factor means: term="Contour" > m1<-aggregate(resp, list(Soils$Contour), mean) > rownames(m1) <- m1[,1] > ( m1 <- m1[,-1] ) pH N Dens P Ca Mg K Na Conduc Depression 4.692 0.08731 1.343 188.2 7.101 8.986 0.3781 5.823 6.946 Slope 4.746 0.10594 1.333 159.4 8.109 8.320 0.4156 6.103 6.964 Top 4.570 0.11256 1.272 150.9 8.877 8.088 0.6050 4.872 5.856 > > # 2-factor means: term="Contour:Depth" > m2<-aggregate(resp, list(Soils$...

Hi, I recently built my new graphics workstation with the express idea of using Linux, Wine, and an AutoCAD clone called IntelliCAD. http://intellicadms.com/ Running Ubuntu 9.04 (Jaunty) on my old Intel system with Wine works just fine. IntelliCAD functions beautifully. However, when I installed the above on my all-new AMD Phenom II X4 955 system, IntelliCAD simply does not launch. I use: -

*Dear experts,* Please excuse me for disturbing... Right now I am struggling with GLM a bit... Would you be so kind to provide me a solution on using nuisance variables. The problem is that I have data on Depression (volumetric measurements of different brain regions) and I want to include age, gender and education as nuisance parameters in the model. In the other words I would like to model the model depression scores in respect to volumes, considering the mentioned variables as nuisance ones. Could you plea...

Hi, Here is (I hope) all the relevant output from R. > mean(s1$ZDEPRESSION, na.rm=T) [1] -1.041546e-16 > mean(s1$ZDIVERSITY_PA, na.rm=T) [1] -9.660583e-16 > mean(s1$ZMEAN_PA, na.rm=T) [1] -5.430282e-15 > lm.ridge(ZDEPRESSION ~ ZMEAN_PA * ZDIVERSITY_PA, data=s1)$coef ZMEAN_PA ZDIVERSITY_PA ZMEAN_PA:ZDIVERSITY_PA -0.3962254 -0.3636...

...sign is a 2*2 factorial with 10 replicates per treatment. Looking at the code for summary.manova, the error involves a problem with qr(). Does anyone have a suggestion as to how to deal with this error? The analysis works fine when there are only 32 variables. Thanks, Simon. Simon Blomberg, PhD Depression & Anxiety Consumer Research Unit Centre for Mental Health Research Australian National University http://www.anu.edu.au/cmhr/ Simon.Blomberg at anu.edu.au +61 (2) 6125 3379

Hi all, I want to make a contingency table in R. I want to tabulate two variables, one as the independent and second as the dependent variable. The IV has two categories, namely, birth complications, and no birth complications. The frequency of birth complication category is fifty, and the frequency of no birth complication category is 34. The categories and frequencies of DV follows.

...ope and intercept at the second level of a 3-level model. Specifically, I am testing a 3-level model in which time (WEEK) is nested in participants (PARTICIP) and participants are nested in dyads (DYADID). The goal is to examine how an interpersonal style (CORUMTO) one week predicts changes in depression the following week (BDIAFTER) controlling for levels of depression (BDI) from the previous week. I have approximately 80 dyads and approximately 160 participants. When modeling a random slope and intercept at the second level (the participant or person level), I receive output, and I also rece...

Thanks, I was getting to try this, but got side tracked by actual work... Your analysis reproduces the SPSS unscaled estimates. It still remains to figure out how Nick got > coefficients(lm(ZDEPRESSION ~ ZMEAN_PA * ZDIVERSITY_PA, data=s1)) (Intercept) ZMEAN_PA ZDIVERSITY_PA ZMEAN_PA:ZDIVERSITY_PA 0.07342198 -0.39650356 -0.36569488 -0.09435788 which does not match your output. I suspect that ZMEAN_PA and ZDIVERSITY...

...lm() gives different results to lm.ridge() and SPSS Dear Nick, On 2017-05-05, 9:40 AM, "R-devel on behalf of Nick Brown" <r-devel-bounces at r-project.org on behalf of nick.brown at free.fr> wrote: >>I conjecture that something in the vicinity of >> res <- lm(DEPRESSION ~ scale(ZMEAN_PA) + scale(ZDIVERSITY_PA) + >>scale(ZMEAN_PA * ZDIVERSITY_PA), data=dat) >>summary(res) >> would reproduce the SPSS Beta values. > >Yes, that works. Thanks! That you have to work hard in R to match the SPSS results isn?t such a bad thing when you facto...

Hi all, Is everyone switched over to the Google lists? I''m basically ready for the switch, I think, and doing so would be one more closed loop for me. I notice we''ve only got a bit over 300 subscribers to the Puppet Users list, while my madstop list has nearly 600. This is a bit depressing, but it probably represents a better idea of who''s actually using the

I''m trying to set up cucumber so we can try it in some Ruby automation scripts I''m working on. These will be straight Ruby, for use in BBEdit and TextMate (on the Mac). NOT Rails. The big early-stage stumbling block I''m having is setting up the initial Rakefile and directory structure. I think I''m rtfm''ing, but clearly I''m missing

...mm correlation structure is modelled, rather than known. I have also searched the R archives, but no luck. I think it should be of the form gls(response~A*B*C, data=dat, correlation=corSymm(...?)) but I don't understand the arguments to corSymm. Thanks in advance, Simon. Simon Blomberg, PhD Depression & Anxiety Consumer Research Unit Centre for Mental Health Research Australian National University http://www.anu.edu.au/cmhr/ Simon.Blomberg at anu.edu.au +61 (2) 6125 3379

...9 responses (pH, N, ... Conduc). [An R-readable version of the data is appended at the bottom.] > soils <- read.delim("soils.dat") > str(soils) 'data.frame': 48 obs. of 14 variables: $ Group : int 1 1 1 1 2 2 2 2 3 3 ... $ Contour: Factor w/ 3 levels "Depression","Slope",..: 3 3 3 3 3 3 3 3 3 3 ... $ Depth : Factor w/ 4 levels "0-10","10-30",..: 1 1 1 1 2 2 2 2 3 3 ... $ Gp : Factor w/ 12 levels "D0","D1","D3",..: 9 9 9 9 10 10 10 10 11 11 ... $ Block : int 1 2 3 4 1 2 3 4 1 2...

Both of these applications fulfill a great need of mine: to read data directly from google spreadsheets that are private to myself and one or two collaborators. Thanks to the authors. I had been using RGoogleDocs for the about 6 months (maybe more) but have had to stop using it in the past month since for some reason that I do not understand it no longer reads google spreadsheets. I loved it. Its

Dear Nick, On 2017-05-05, 9:40 AM, "R-devel on behalf of Nick Brown" <r-devel-bounces at r-project.org on behalf of nick.brown at free.fr> wrote: >>I conjecture that something in the vicinity of >> res <- lm(DEPRESSION ~ scale(ZMEAN_PA) + scale(ZDIVERSITY_PA) + >>scale(ZMEAN_PA * ZDIVERSITY_PA), data=dat) >>summary(res) >> would reproduce the SPSS Beta values. > >Yes, that works. Thanks! That you have to work hard in R to match the SPSS results isn?t such a bad thing when you factor in t...