similar to: structural equation modeling in sem, error, The model has negative degrees of freedom = -3, and The model is almost surely misspecified...

I've been reading the R mail archives and I've found a lot of messages with this same kind of problem, but I can't understand the answers. Can one of you try to explain this to me? Here's my example. Given a regression model and a variable, I want to use unstack() on the vector of residuals and make some magic with the result. But unstack hates me. PCSE <- function

I'm wondering how to save an object (models like lm, loess, etc) in a DB to retrieve and use it afterwards, an example: wind_ms <- abs(rnorm(24*30)*4+8) air_kgm3 <- rnorm(24*30, 0.1)*0.1 + 1.1 wind_dg <- rnorm(24*30) * 360/7 ms <- c(0:25) kw_mm92 <- c(0,0,0,20,94,205,391,645,979,1375,1795,2000,2040) kw_mm92 <- c(kw_mm92, rep(2050, length(ms)-length(kw_mm92))) modelspline

Hi It clear that to do a classification with svm under 10-fold cross validation one uses svm(Xm, newlabs, type = "C-classification", kernel = "linear",cross = 10) What corresponds to the nnet? nnet(.....,cross=10)? Regards

Dear R Community, I am currently facing this seemingly obscure Problem with Panel Corrected Standard Errors (PCSE) following Beck & Katz (1995). As the authors suggest, I regressed a linear model (tmodel) with lm() with option "na.action=na.exclude" (I have also tried other options here). My dataframe is organized in pooled times series fashion, but with various missing values

Hi there, My question is about the 'theta' parameter in specification of a NB GAMM. I have fit a GAM with an optimum structure of: SB.gam4<-gam(count~offset(vol_offset)+ s(Depth_m, by=StnF, bs="cs")+StageF*RegionF, family=negbin(1, link=log), data=Zoop_2011[Zoop_2011$SpeciesF=='SB',]) However, this GAM shows heterogeneity in the

I wonder if someone could explain why, when I perform confirmatory factor-analysis model using polychoric correlations why I do not get an estimated confidence interval for the RMSEA. My experience with these type models is that I would obtain a confidence interval estimate. I did not get any warning messages with the output. RESULTS: Model Chisquare = 1374 Df = 185 Pr(>Chisq) = 0

After spending a few hours reading dozens of blogs and tutorials I am still stumped. Numerous sites mention ruby 1.8.7 and rails 2.3 supporting TLS out of the box. However, I can not find any documentation on setting this up to send emails correctly. Furthermore, after attempting to send emails using Gmail and TLS i get the following error: 530 5.7.0 Must issue a STARTTLS command first.

Dear all, I'd like to discuss about a possible bug in function StructTS of stats package. It seems that the function returns wrong value of the log-likelihood, as the added constant to the relevant part of the log-likelihood is misspecified. Here is an simple example: > data(Nile) > fit <- StructTS(Nile, type = "level") > fit$loglik [1] -367.5194 When computing the

Matt: that’s fair. We’re generally apprehensive of option 2 as well. Eli: Yes, currently we believe that aarch64.neon.addp is the only arm64 one affected, but we don’t know how prevalent this is on other targets. Splitting it is certainly possible combined with the autoupgrader. If disambiguating the intrinsics is the preferred solution, then I think we should also have the langref also specify

On Sun, Nov 24, 2013 at 10:19 PM, Mikael Lyngvig <mikael at lyngvig.org> wrote: > Hi, > > The language reference uses "::" many places where it seems it should be > using ".. code-block:: llvm". The result is that Sphinx does not colorize > these samples. > > Should I fix this? Or is it purposeful? > Could you give an example of what you are

Hi, The language reference uses "::" many places where it seems it should be using ".. code-block:: llvm". The result is that Sphinx does not colorize these samples. Should I fix this? Or is it purposeful? Also, anybody but me that has noticed that Sphinx does not always colorize things even if you put the proper ".. code-block:: llvm" in front? If it is a known

Hi, I am trying to fit the following model: sr.reg.s4.nore <- survreg(Surv(age_sym4,sym4), as.factor(lifedxm), data=bip.surv) Where age_sym4 is the age that a subject develops clinical thought problems; sym4 is whether they develop clinical thoughts problems (0 or 1); and lifedxm is mother's diagnosis: BIPOLAR, MAJOR DEPRESSION, or CONTROL. I am interested in whether or not

Hi! I want to construct a path analysis model that can account for measurement error in totally aggregated parcels, which refer to parcels where all of the items in a scale are summed or averaged. If I am not mistaken, Bollen (1989) advocates the following formula for computing the error variance of each parcel: (1−α(parcel))×variance(parcel), such that α refers to Cronbach's alpha, which

Hi all, I am trying to set up a simple path analysis in the SEM package, but I am having some trouble. I keep getting the following error message or something similar with my model, and I'm not sure what I'm doing wrong: Error in solve.default(C) : system is computationally singular: reciprocal condition number = 2.2449e-20 In addition: Warning message: In sem.default(ram = ram, S = S,

Just upgraded to 5.4 yesterday morning. (Kudos to those who did that!) That's probably irrelevant to this problem: "yum install mplayer" fails. Here's what happens: ================================================== # yum install mplayer Loaded plugins: fastestmirror, filter-data, kernel-module, keys, kmod, list- : data, merge-conf, security, upgrade-helper,

Dear All: Do you have any experience on calling winbugs in R enviorment? I uploaded R2winbgus, pktools, BRugs, ect. in order to call winbugs. In the end, I try to use bugs(). The error message still displays that bugs() can not find. Do you have any similar issue? If so, could you give your advice? Many Thanks! Xin [[alternative HTML version deleted]]

Full_Name: Barnet Wagman Version: 63.2 OS: Linux Submission from: (NULL) (207.229.172.168) After a number of interations, glm(family=binomial) throws an error 'Error: names attribute must be the same length as the vector' The same model works in lm() with the same data frame, so presumably it is not misspecified. In the iterations prior to the error, glm() generates the warnings

Greetings, Fellow R-ians: I'm working with a few different quasi-likelihood formulations for some data I'm analyzing. I'd like to implement the Information Matrix Test (see, e.g., White, 1982, or Lancaster, 1984) for each of them to determine which of the models is more likely. Since the null distribution of the test statistic is chi-squared, I envision calculating the ordinate of

ggplot2 =================================== ggplot2 is a plotting system for R, based on the grammar of graphics, which tries to take the good parts of base and lattice graphics and avoid bad parts. It takes care of many of the fiddly details that make plotting a hassle (like drawing legends) as well as providing a powerful model of graphics that makes it easy to produce complex multi-layered

ggplot2 =================================== ggplot2 is a plotting system for R, based on the grammar of graphics, which tries to take the good parts of base and lattice graphics and avoid bad parts. It takes care of many of the fiddly details that make plotting a hassle (like drawing legends) as well as providing a powerful model of graphics that makes it easy to produce complex multi-layered