similar to: general location model

Dear all, I am implementing a redundant variables F-test. For that I need to compute 2 models, restricted and unrestricted, then extracting the residuals to calculate the test statistic. I borrowed this elegant solution from the LMtest package to rebuild the first of the matrices involved (the unrestricted one) on the basis of model spec. and data >red.var.test<-function(formula,

Hello all, When I ran my R program in R version 2.1.1 in windows XP and 2000, and I downloaded a package call "MIX", I got an application error: RGui:Rgui.exe - Application Error: The instruction at "0x1009d8a1" referenced memory at "0xbde48f58". The memory could not be "read". Click on ok to terminate the program click on cancel to debug the program

Can anyone recommend a package that can be used to randomize subjects? I am looking for a generalized package, or several packages that can accomplish unrestricted randomization (i.e. simple random assignment) restricted randomization including stratified randomization, blocked randomization, and adaptive randomization. Thanks, John John David Sorkin M.D., Ph.D. Chief, Biostatistics and

Dear list members, I'm trying to use DCC-GARCH model to estimate the correlation. I have downloeaded ccgarch packeage but can't understand some argument in the formula. dcc.estimation(inia, iniA, iniB, ini.dcc, dvar, model, method="BFGS", gradient=1, message=1) which is on R.Help I understand others except "ini.dcc" which is described as "a vector of initial

Hi All, I a C function code, which has been successed in integrating with Splus on Unix. Now I want to move it on windows. So I want to integrate this C function with R on windows XP. I download tools.zip and saved it in c:\bin directory, download MinGW and saved it in c:\MinGw directory, download Perl and saved it in c:\Perl directory. My path is set as:

Daniel, When Group is entered as a factor, and the factor has two levels, the ANOVA table gives a p value for each level of the factor. What I am looking for is the omnibus p value for the factor, i.e. the test that the factor (with all its levels) improves the prediction of the outcome. You are correct that normally one could rely on the fact that the model

Dear Professor Schaefer I am experiencing a technical difficulty with your mix package. I would appreciate it if you could help me with this problem. When I run the following code, R 2.5.1 and R 2.6.0 crashes. It's been tested on at least 2 windows machine and it is consistent. Execution code it's self was coped from the help file of imp.mix. Only thing I supplied was a fake dataset.

All - I am evaluating an arima(2,1,3) and arima(3,1,3) and notice the log-likelihood of the restricted model is higher than the log-likelihood of the unrestricted. Since these are nested models, I thought the unrestricted model would have a log-likelihood at least as large as that of the restricted model. Am I interpreting the "loglik" output incorrectly? Regards, Stephen

Hello, all. I'm investigating the rate at which skeletal joint surfaces pass through a series of ordered stages (changes in morphology). Current statistical methods in this type of research use various logit or probit regression techniques (e.g., proportional odds logit/probit, forward/backward continuation ratio, or restricted/unrestricted cumulative probit). Data typically include the

Hello, I was studying several packages related to time series analysis (urca, vars, tseries). I understand that we can estimate a VECM and also test restrictions on alphas and betas. However, I couldn't find a function that allows me to specify the five cases of VECM (restricted constant, unrestricted constant, restricted and unrestricted trend and no constant). Is there any function that

Dear All, My question is: how can I estimate VECM system with "unrestricted trend" (aka "case 5") option as a deterministic term? As far as I know, ca.jo in urca package allows for "restricted trend" only [vecm <- ca.jo(data, type = "trace"/"eigen", ecdet = "trend", K = n, spec = "transitory"/"longrun")].

Dear R-users, Half a year ago we put out the R package "glmmADMB" for fitting overdispersed count data. http://otter-rsch.com/admbre/examples/glmmadmb/glmmADMB.html Several people who used this package have requested additional features. We now have a new version ready. The major new feature is that glmmADMB allows Bernoulli responses with logistic and probit links. In addition there

A copy of this question can be found on Cross Validated: https://stats.stackexchange.com/questions/645362 I am estimating a system of seemingly unrelated regressions (SUR) in R. Each of the equations has one unique regressor and one common regressor. I am using `gmm::sysGmm` and am experimenting with different weighting matrices. I get the same results (point estimates, standard errors and

Hi I constructed a mixed-effects model from longitudinal repeated measurements of lab values in 22 patients seperated into two groups with the groups as fixed effect using lme. I thought about using the jackknife procedure, i. e., removing any one subject and calculating the fixed effect, to assess the stability of the fixed effect and thereby validate the model. I suppose this has been done in

Generally speaking, this sort of detailed statistical question about a speccial package in R does not get a reply on this general R programming help list. Instead, I suggest you either email the maintainer (found by ?maintainer) or ask a question on a relevant R task view, such as https://cran.r-project.org/web/views/Econometrics.html . (or any other that you judge to be more appropriate).

https://bugzilla.mindrot.org/show_bug.cgi?id=2713 Bug ID: 2713 Summary: Please provide a StrictModes-like setting (command line parameter) for ssh (client) Product: Portable OpenSSH Version: 7.5p1 Hardware: Other OS: Other Status: NEW Severity: enhancement Priority: P5

Keywords T1 Q.931 isdn disconnect cause codes itu standard libpri Dont know if anyone wondered what q.931 cause codes are but i wishwe could get these back into the dial plan as a var Standard Q931 Codes Decimal Value Hexadecimal Value Definition 1 01 Unallocated (unassigned) number. This number is not in the routing table or it has no path across the ISDN cloud (network). 1. Check routing

I know the call parking feature changed in 1.4.23.1 to fix some serious issues. I'm seeing a major change though which I find disturbing. A person parks a call by transferring it to the parking position (700). When the timeout value is reached, the call is NOT returned to that device, but rather the 's' extension of the phone's registered context (in this case [unrestricted]).

Dear R user, I am a newcomer and need help concerning 'draw a random number for a restricted area of a prareto distribution'. (1) For estimation of pareto distribution: >http://stats.stackexchange.com/questions/27426/how-do-i-fit-a-set-of-data-to-a-pareto-distribution-in-r< We calculate the pareto distribution (parameter estimation) as follows: pareto.MLE <- function(X) { n

Hallo List, I'm trying to implemement a restricted permutation scheme in permutest(). More precisely I have dependence in my data that should be allowed for in the permutation - I simulated the problem in the example of the vegan documentation p.24: library(vegan) data(varespec) ## Bray-Curtis distances between samples dis <- vegdist(varespec) ## First 16 sites grazed, remaining 8 sites