similar to: mixed-effects model using lmer

The following fit does not make sense to me, please, correct me if I have a logical error. > moddowsn Call: coxph(formula = Surv(start, stop, resp) ~ sn + matfac2, data = coxsn1, method = "efron") coef exp(coef) se(coef) z p sn2 0.0497 1.051 0.02030 2.450 1.4e-02 sn3 -0.0532 0.948 0.02038 -2.610 9.0e-03 sn4 -0.0410 0.960

Dear R members, I am new in using frailty models in survival analyses and am getting some contrasting results when I compare the Wald and likelihood ratio tests provided by the r output. I am testing the survivorship of different sunflower interspecific crosses using cytoplasm (Cyt), Pollen and the interaction Cyt*Pollen as fixed effects, and sub-block as a random effect. I stratified

Hello, I acquired a Diva Server 4BRI for an AsteriskPBX system and installed a 2.4.20 kernel modified with the latest eicon capi2.0 support from www.melware.net. Capi support seems to work fine: Asterisk+chan_capi answers the calls. However if I try to use my old isdnutils apps by loading a i4l compatibility layer then I have the following problems: - if I use the diva2i4l module and start the

Hello, I'm using the "sem" package to do a confirmatory factor analysis on data collected with a questionnaire. In the model, there is a unique factor G and 23 items. I would like to calculate the standardized residual variance of the observed variables. "Sem" only gives the residual variance with the "summary" function, or the standardized loadings with the

Hi everybody, I'm trying to fit a Garch(1,1) process to the DAX returns. My data consists of about 2300 10day-logreturns in chronologically descending order (see attachment). But if I use the garch function I get a very high alpha_1 and a quite low beta, which doesn't make that much sense. I think I am missing something, but have no idea what it might be. I'd appreciate it a lot

Dear experts, I'm attempting to solve a constrained optimization problem using the Rdonlp2 package. I created a Lagrange function (L=f(x)-lambda(g(x)-c)), where x is a vector of 16 parameters. This is what I'm using as objective function in the code below. In addition, I set bounds on these parameters (par.u and par.l). When I run the code, I get the error message shown below. Any idea

Hi, I'm trying to create a plot showing the density distribution of some shipping data. I like the look of violin plots, but my data is not continuous but rather binned and I want to make sure its binned nature (not smooth) is apparent in the final plot. So for example, I have the number of individuals per vessel, but rather than having the actual number of individuals I have data in the

Full_Name: Gang Liang Version: 1.4.1 OS: mandrake-8.2 Submission from: (NULL) (128.32.81.135) tseries version: 0.9-1 quadprog version: 1.4-4 mva version: 1.4.1 version: platform i686-pc-linux-gnu arch i686 os linux-gnu system i686, linux-gnu status major 1 minor 4.1 year 2002 month 01

Greetings, I have a data set look like this: > print(MyResult) 2 5 9 10 11 13 mean 6.108394 3.272211e+01 3.951703e+02 2.086217e+03 79.838867 2.592272e+04 median 1.288889 2.743511e+00 4.061491e+00 4.943710e+01 -2.821667 -3.040625e+00 var 2885.350595 2.343144e+05 1.022680e+06 1.200589e+07 365179.508831

On 20/11/2017 6:38 AM, Maria Lathouri via R-help wrote: > Dear all > I know that it is a very simple question but it seems that I cannot change the colour in the bars. > I have the following dataframe: > A ? ? ? ? ? ? ? ?? B ? ? ? ? ? ? ? ? C ? ? ? ? ? D ? ? ? ?? E ? ? ? ? ? ? ? ? ?? F ? ? ? ? ? ? ? ? ?? G ? ? ? ? ? ? ? ?0.0.24 ? ? ? ?? 152460 ? ? ? ? 474 ? ? ? 5.5 ? ? ?? 612000 ? ?

I need a function like signif(), but returns the rounded values as character strings, formatted with trailing zeros where appropriate. If anyone has one, I would sure appreciate a copy. Thanks -Don Details: signif() rounds a number to a specified number of significant digits, for example: > x <- c(2.503,2.477,0.1204) > signif(x[1],3) [1] 2.5 > signif(x[2],3) [1] 2.48

On 20/11/2017 7:09 AM, Maria Lathouri wrote: > Dear all, > > I am really sorry for this. I have attached the script and a .csv file > with an example. The problem is that you have dat as a matrix with only one row. barplot() applies the colours to the rows of dat; since you have only one row, everything ends up blue. You can fix this by creating dat as a vector using

Dear all,? I am really sorry for this. I have attached the script and a .csv file with an example.? Hope this will help. Many thanks,Maria ???? 11:53 ?.?. ???????, 20 ????????? 2017, ?/? Duncan Murdoch <murdoch.duncan at gmail.com> ??????: On 20/11/2017 6:38 AM, Maria Lathouri via R-help wrote: >? Dear all > I know that it is a very simple question but it seems that I cannot

I have a model with quant vars only and the error message does not make sense: (mod1 <- coxph(Surv(time=strt,time2=stp,event=(resp==1))~ +incpost+I(amt/1e5)+rate+strata(termfac), subset=dt<"2010-08-30", data=inc,method="efron")) Call: coxph(formula = Surv(time = strt, time2 = stp, event = (resp == 1)) ~ +incpost + I(amt/1e+05) + rate + strata(termfac),

Dear WizaRds, A run of nls produces the following concise summary: > summary(cs.wt) Formula: 0 ~ wt.MM(conc, time, A1, a1, A2, a2) Parameters: Estimate Std. Error t value Pr(>|t|) A1 4.814e+02 2.240e+01 21.495 0.0296 * a1 7.401e-01 7.435e-02 9.956 0.0637 . A2 1.613e+02 1.738e+01 9.280 0.0683 . a2 1.770e-02 7.324e-03 2.417 0.2497

I see a package "x246-0.120-5.20120303.el6 (i686)" available in my Add/Remove Software tool. I can't find said package in any of the repos I've searched through. Can anyone point me to a link where I can download that package to my local repo? -- _ ?v? /(_)\ ^ ^ Mark LaPierre Registerd Linux user No #267004 https://linuxcounter.net/ ****

Dear all I know that it is a very simple question but it seems that I cannot change the colour in the bars. I have the following dataframe:? A ? ? ? ? ? ? ? ?? B ? ? ? ? ? ? ? ? C ? ? ? ? ? D ? ? ? ?? E ? ? ? ? ? ? ? ? ?? F ? ? ? ? ? ? ? ? ?? G ? ? ? ? ? ? ? ?0.0.24 ? ? ? ?? 152460 ? ? ? ? 474 ? ? ? 5.5 ? ? ?? 612000 ? ? ? ? ? 59061000 ? ? ? 1540313 and here is the script:

Dear list members, I am trying to estimate parameters of the AR(1)-GARCH(1,1) model. I have one additional dummy variable for the AR(1) part. First I wanted to do it using garchFit function (everything would be then estimated in one step) however in the fGarch library I didn't find a way to include an additional variable. That would be the formula but, as said, I think it is impossible to add

I currently work on a draft of an aquatic bioassessment. The conditions tested are the following: ER river water T dechlorinated water control 0.5 + 0.5mg / L of malate T + 1 dechlorinated water control + 1g / L of malate T ED dechlorinated water control SED + ER + river water sediment SED ED + sediment + water dechlorinated. It is the result of AChE in muscle (fillet of fish). The production of

Hello, I am hoping for some advice regarding warning/error messages I received when running a Cox regression # message 1 - obtained while creating a plot of residuals > plot (NV.zph, main = "groupNUSM - UNFIT", var= 'groupNUSM') Warning messages: 1: In approx(xx, xtime, seq(min(xx), max(xx), length.out = 17)[2 * : collapsing to unique 'x' values 2: In