similar to: Translating NLMIXED in nlme

minimal reproducible example plot(1,1) identify(1,1) # or locator() now, trying to close the window by clicking on the cross of the upper right corner causes Rgui (and Rterm) to crash. I see the same behaviour on 2 different Windows PC (one with Win 8.1 and one with Win 10). I did not see the problem in linux (see below) WINDOWS ************** > sessionInfo() R version 3.3.0 beta

Thanks, I'll track this down. Duncan Murdoch On 05/04/2016 9:35 AM, Simone Giannerini wrote: > minimal reproducible example > > plot(1,1) > identify(1,1) # or locator() > > now, trying to close the window by clicking on the cross of the upper > right corner causes Rgui (and Rterm) to crash. > > I see the same behaviour on 2 different Windows PC (one with Win 8.1

On 05/04/2016 11:56 AM, Henrik Bengtsson wrote: > If of any help, > > I can reproduce this (on Windows 7) back to at least R 3.0.3 but it's > not there in R 3.0.0. (I have *not* checked with R 3.0.1 and 3.0.2 > which I don't have installed). That doesn't necessarily mean that 3.0.0 was fine. It's a segfault (I'd guess some memory being accessed after being

If of any help, I can reproduce this (on Windows 7) back to at least R 3.0.3 but it's not there in R 3.0.0. (I have *not* checked with R 3.0.1 and 3.0.2 which I don't have installed). /Henrik On Tue, Apr 5, 2016 at 8:23 AM, Duncan Murdoch <murdoch.duncan at gmail.com> wrote: > Thanks, I'll track this down. > > Duncan Murdoch > > > On 05/04/2016 9:35 AM,

dear all, I report a problem very simple, that I does non know how to handle. look at the following code: > a = rep(16.256, 5) > sum(a[1:5]^2) - (sum(a[1:5])^2/5) [1] 2.273737e-13 as you can see i retrieve a non 0 value, when i am expected to. what can I do? > sessionInfo() R version 2.6.2 (2008-02-08) i386-pc-mingw32 locale:

Dear list, I do not know what would be the 'correct' answer to the following but I think that they should return the same value to avoid potential problems and hard to debug errors. Regards, Simone --------------------------------------- > NCOL(NULL) [1] 1 > NROW(NULL) [1] 0 > sessionInfo() R version 4.3.1 RC (2023-06-08 r84523 ucrt) Platform: x86_64-w64-mingw32/x64 (64-bit)

Hello, I just installed R-2.15.0 on windows XP and cannot load package tcltk2 (which I just downloaded from CRAN as tcltk2_1.2-1.zip; package install reported no problems): > library(tcltk2) Carico il pacchetto richiesto: tcltk Loading Tcl/Tk interface ... done Error : .onLoad failed in loadNamespace() for 'tcltk2', details: call: system("cat /etc/issue", intern = TRUE,

According to ?parCapply: parRapply and parCapply always return a vector. This appears not to be the case in the following minimal reproducible example: > library(parallel) > nslaves <- 2 > cl <- makeCluster(nslaves) > X <- matrix(2,nrow=3,ncol=4) > X <- rbind(c(1,1,0,1),X) > tv <- parCapply(cl,X,FUN=function(x){ +

Dear all, I think I have found a buglet in plot.ts plot.ts(x=1,type="n") # correct: does not show the plot plot.ts(x=1,y=1,type="n") # not correct: does show the plot I did not investigate the problem in depth but it could be related to the switch xy.labels, in fact plot.ts(x=1,y=1,type="n",xy.labels=TRUE) # does show the plot

Dear all, it looks like that tabulate() does not check for the bounds of the input. Reproducible example: > b <- 1:2 > tabulate(b[1:100]) [1] 1 1 > R.version _ platform i386-pc-mingw32 arch i386 os mingw32 system i386, mingw32 status Patched major 2 minor 11.1 year 2010 month 09 day

Dear all, from the help page of quantile: "x ??? numeric vectors whose sample quantiles are wanted. Missing values are ignored." from the help page of IQR: "x ??? a numeric vector." as a matter of facts it seems that both quantile() and IQR() do not check for the presence of a numeric input. See the following: set.seed(11) x <- rbinom(n=11,size=2,prob=.5) x <-

R-Listers: I am doing a quasi-maximum likelihood estimation and I get a "subscript out of bound" error message, Typically I would think this means that a subscript used in the function is literally out of bounds however I don't think this is the case. All I change in the code is a constant, that is hard-wired in (not data dependent and not parameter dependent), furthermore,

Hi, I was very much hoping someone could help me with the following. I am trying to convert some SAS NLMIXED code to NLME in R (v.2.1), but I get an error message. Does anyone have any suggestions? I think my error is with the random effect "u" which seems to be parametrized differently in the SAS code. In case it's helpful, what I am essentially trying to do is estimate parameters

To account for likely differences between families in naturalization rates, we fitted a generalized linear mixed model, using PROC NLMIXED in SAS10, with the naturalization rate per genus (that is, the number of naturalized species in a genus as a proportion of the total number of introduced species in a genus) as the response variable, a variable coding genera as containing at least one native

Hello! Can anyone help me. I am using the posterior.mode from the result of glmmML. It apears to be different from the BLUe estimate of the RANDOM statement in PROC NLMIXED in SAS. Why is that? Thank you Ronen [[alternative HTML version deleted]]

Dear R friends A while a go I sent an email to the epi-list and later to the help-list and no answer could fully illuminate my question. So Im trying again with a more specific matter. Im trying to work on a script (function) to analyse data from a diagnostic test meta-analysis with random effects. This was first described by an author using SAS witn PROC NLMIXED. Im not an expert in R and much

I'm trying to run a zero-inflated regression for a continuous response variable in R. I'm aware of a gamlss implementation, but I'd really like to try out this algorithm by Dale McLerran that is conceptually a bit more straightforward. Unfortunately, the code is in SAS and I'm not sure how to re-write it for something like nlme (if at all possible - with conditions etc). Does

On Windows XP SP3 > format(pippo, format="%I:%M:%S %p") [1] "11:43:12 PM" In any case, if there is a problem it is in Windows, not in R, since this is an Windows service that is being called. So there is no evidence at all of any bug in R here. On Wed, 10 Sep 2008, pego.enrico at tiscali.it wrote: > Full_Name: Enrico Pegoraro > Version: 2.7.2 > OS: Windows

Hi, I am trying to convert the following NLMIXED code to NLME, but am running into problems concerning 'Singularity in backsolve'. As I am new to R/S-Plus, I thought I may be missing something in the NLME code. NLMIXED *********** proc nlmixed data=kidney.kidney; parms delta=0.03 gamma=1.1 b1=-0.003 b2=-1.2 b3=0.09 b4=0.35 b5=-1.43 varu=0.5; eta=b1*age+b2*sex+b3*gn+b4*an+b5*pkn+u;

Dear All I am trying to estimate a mixture of regression and get the following error using the mixreg package: Error in y - yhat : non-conformable arrays The instruction I used were: x <- as.matrix(LRHUN) y <- as.matrix(LRINTER) TS <- list(list(beta=c(3.0,1.0),sigsq=1,lambda=0.4), list(beta=c(0.0,1.0),sigsq=1,lambda=0.6)) prova <- mixreg(x,y, ncomp=2, theta.start=TS)