search for: depress

...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,.1...

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,.1...

...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...

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 id="db_bfilter1" class="raised"><%= filter_link_helper "This month" %></span> when hit, the link execute a new filtered search (by week or by month according to th...

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 l...

...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(Soi...

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 p...

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.3...

...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

...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. Schizophrenic 28, depressed 26, normal 30. When I am trying to make a contingency table in R, using table(name of variable one,name of variable two), I am getting an error that all arguments must have the same length. I believe that there would be two rows and three columns according to categories of IV and DV, But I guess...

...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 r...

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 ZDIVERS...

...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 fa...

..., 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 list. As a refresher, here are all of the lists, and how to join without a Google account: On Feb 13, 2008, at 11:30 AM, Luke Kanies wrote: > > In the meantime, here are the lists: > > http://group...

...ld standard directory tree (3) clone in standard Rakefile (4) get started speccing. (2) and (3) are my problem. Thanks much! ----- Bill Tozier AIM:vaguery at mac.com ? Twitter:Vaguery http://williamtozier.com/slurry "The brotherhood of man is not a mere poet''s dream; it is a most depressing and humiliating reality." -- Oscar Wilde

...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...

...hat 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 loss depresses me. I started using RGoogleData which works. I have noticed that both packages read data slowly. RGoogleData is much slower than RGoogleDocs used to be. Both seem a lot slower than if one manually downloaded a google spreadsheet as a csv and then used read.csv function - but then I would not be...

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 i...