similar to: Degree of freedom for contrast t-tests in lme

Dear R users I would like to ask a question regarding to icc (intraclass correlation) or many biologists refer it to as repeatability. It is very useful to get icc for many reasons and it is easy to do so from linear mixed-effects models and many packages like psy, psychometric, aod and irr have functions to calculate icc. icc = between-group variance/(between-group variance + residual

Dear R users I would like to simulate underdispersed Poisson and binomial distributions somehow. I know you can do this for overdispersed counterparts - using rnbinom() for Poisson and rbetabinom() for binomial. Could anyone share functions to do this? Or please share some tips for modifying existing functions to achieve this. Thank you very much for your help and time Shinichi

Hi, Does anybody have a good tip of how to treat bimodal data to perform statistical analyses? My data set ranges from -1 to 1 (any values are posssible in between) and most data are either close to -1 or close to 1. They are the results of a two choice experiment where individuals could choose more than once in either direction and scores were calculated. Simone Simone Immler

I have been trying to make use of the odesolve library on my university's Linux grid - currently R version 2.0.1 is installed and the system runs 64-bit Scientific Linux based on Redhat. I cannot seem to get lsoda working when I define the model as a shared C library. For example, the following snippet uses the mymod.c example bundled with the package: ### START rm(list=ls())

Dear all - I have followed the thread the reply to which was lead by Thomas Lumley about using pweibull to generate fitted survival curves for survreg models. http://tolstoy.newcastle.edu.au/R/help/04/11/7766.html Using the lung data set, data(lung) lung.wbs <- survreg( Surv(time, status)~ 1, data=lung, dist='weibull') curve(pweibull(x, scale=exp(coef(lung.wbs)),

Hi all, I have a questions regarding mixed effects models: I'm trying to plot e.g. fitted vs residuals for each level of the random effects, and i'm getting the same error. I guess this might be a problem of the graphic capabilities of R. Is there any way to obtain those plots? library(nlme) attach(ergoStool) names(ergoStool) [1] "effort" "Type"

Let's try this again. I have a Quintum AX Tenor gateway sending calls to Asterisk from BT analogue lines connected to FXO. The agents hear an echo on their side but incoming callers hear the conversation fine. I can't seem to find the problem. Anyone seen this issue before? <p style="margin: 0; padding: 0; border-collapse: collapse; font-family: Tahoma, Arial, Sans-Serif;

I tried rsync 2.6.1cvs on my FreeBSD 4.X boxes. It dumps core when preserve_hard_links and backup. And SUPPORT_HARD_LINKS is true. rsync -avH --delete --backup --backup-dir=/some/backup/dir \ srchost::dir/ dir/ The core says... ------------------------------------------------------------------------ Program terminated with signal 11, Segmentation fault. #0 0x80536c0 in keep_backup

Hi i am having a problem with the 'integrate' function the function i want to integrate has the form sum(vector^x) i have defined the function with a for loop first - integrandtotest <- function(x) {a<-rep(0,len=2) for (i in 1:2) {a[i]<-t[i]^x} sum(a) } the results gives errors ########### Error in integrate(integrandtotest, lower = 0.1, upper =

Dear all, I have a large data set and would like to fit a logistic regression model using the bigglm function. I need to include an offset in the model but when I do this the bigglm function seems to ignore it. For example, running the two models below produces the same model and the offset is ignored bigglm(y~x,offset=z,data=Test,family=binomial(link = "logit"))

double-colon mode doesn't work in CVS version. -vv says port number is 0. % rsync -vv host.domain::any opening tcp connection to host.domain port 0 rsync: failed to connect to host.domain: Can't assign requested address (49) rsync error: error in socket IO (code 10) at clientserver.c(94) And --port=873 cannot be used in client mode now. I made a small sample patch. (Sorry, I

Dear all, I am trying to calculate the proximity matrix for a data set with 16 variables and 6804 observations using random forests. I have a Pentium 4, 3.00GHz processor with 1 GB of RAM. When I use the command randomForest(data.scale,proximity=T) I get the warning message Error: cannot allocate vector of size 361675 kb Is this because I have reached the limit of what my computer is

I think this is one for Gabor. I don't seem to be able to find my way to an answer despite numerous rereadings of factor and table. Here is a toy example: ### Some data EthnicCode <- c("European/Other", NA, "European/Other", "European/Other", "Pacific", "European/Other", "European/Other",

Dear All, we've been having some intermittent problems with students who's home directories are on a Sun E450 running Samba 2.0.6. I'm not totally convinced that the problem is caused by the server, but may be the network. However, I've had a look at the log.smb file and I'm seeing some error messages that I've not come across before: [2001/02/23 18:20:40, 0]

Dear R experts I need to write a function that incorporates double summation, the problem being that the upper limit of the second summation is the index of the first summation, i.e: sum_{j=0}^{x} sum_{i=0}^{j} choose(i+j, i) where x variable or constant, doesn't matter. The following code obviously doesn't work: f=function(x) {j=0:x; i=0:j; sum( choose(i+j,i) ) } Can you help? Thanks

Now that I am at least mapping the right repository to my Debian version (doh!) I have another question. In the installation page at https://cran.r-project.org/bin/linux/debian/#debian-stretch-stable there is the line: "Please note that R packages from the Debian stretch distribution are not compatible with R 3.5.x, as it provides r-api-3.5, while the stretch packages depend on

# Your mailer is set to "none" (default on Windows), # hence we cannot send the bug report directly from R. # Please copy the bug report (after finishing it) to # your favorite email program and send it to # # r-bugs@r-project.org # ###################################################### Dear R developers, this morning I started experiencing a bizarre problem with the

Hello, I've been running some HLMs using the lme function quite happily; it does what I want and I'm pretty sure I understand it. The issue is that I'm currently trying to estimate a model with a 14-level "nusiance" factor as an independent variable...which makes the output quite ugly. All I'm really interested in is the question of whether these factor as a whole

Wow, this has amplified my respect and gratitude for the R Debian repositories. Doing all the install.packages() and following the clues about missing dependencies has taken hours! However, I have almost all the packages I want installed on this machine and Rstudio and the shiny server are installed and running fine. Bizarrely, to my mind, I am totally failing on rJava and hence a number of

Dear All, We are seeing the following warning messages in event logs on PCs that are accessing a Samba server (Solaris 2.7, samba 2.0.6):- 03/03/04/2001,14:44:04,Rdr,Warning,None,3025,N/A,CIS-C033-04,A write-behind operation has failed to the remote server ilex. The data contains the amount requested to write and the amount actually written. What does this message actually mean? Does it