similar to: Forcing coefficients in lm object

I'm using the 'pnt' C function from Rmath library in some C-code. How can I catch the warning message: "full precision was not achieved in 'pnt'" in R. I call the function using the .C(). (options(warn=-1) didn't work) Thanks in advance -- Maarten van Iterson Center for Human and Clinical Genetics Leiden University Medical Center (LUMC) Research Building,

Dear list, If have some repeated measurement data which looks something like: time <- rep(1:5 , each=2*4) groups <- rep(c("Case", "Control"), each=4) subjects <- factor(rep(1:(2*4), 5)) responses <- time + rnorm(5*2*4) + as.integer(factor(groups)) data <- data.frame(responses, time, groups, subjects) Now I want to plot each subject in a separate panel:

Dear R users, Suppose I have the following data.frame: myID myType myNum1 myNum2 myNum3 a Single 10 11 12 b Single 15 25 35 c Double 22 33 44 d Double 4 6 8 and I want to have new records: myID myType myNum1 myNum2 myNum3 e Single 12.5 18

Hi, Can someone of you tell me why these two simple examples works as "expected", plot(c(1,1,2,2)~c(1,2,1,2),col=c(1,2,3,4),pch=c(1,2,3,4)) coplot(c(1,1,2,2)~c(1,2,1,2)|c("A","A","B","B"),col=c(1,2,3,4),pch=c(1,2,3,4 )) but when I try it with xyplot in the lattice package it fails.

Hi all, I should have given a better explanation of my problem. Here it is. I extracted from my code the bit that gives the error. Place this in a file called test.c #include <math.h> #include <R.h> #include <Rmath.h> #include <float.h> #include <R_ext/Print.h> int main(){ double prob[3] = {0.0, 0.0, 0.0}; double prob_tot = 0.; prob[0] = 0.3*dnorm(2, 0,

Hello all, Currently we are looking at the possibility of moving from Courier-IMAP to Dovecot (the one that comes with Debian Lenny: version 1.0.15). I'm currently testing it on a non-production server and have a few questions: 1) Currently all our clients use Thunderbird as their mail client to read mail from our IMAP server. Thunderbird is configured with "INBOX." as the

Full_Name: Jelle Goeman Version: 2.0.0 OS: windows Submission from: (NULL) (145.88.209.33) I get a strange error: > is.vector(1:10) Error: recursive default argument reference What's recursive about is.vector? Kind regards, Jelle Goeman My R: platform i386-pc-mingw32 arch i386 os mingw32 system i386, mingw32 status major 2

Dear useRs! I would like to generate a list of all possible (unique) partitions of N units into k classes. For example, all possible partitions of 4 units into 2 classes are (I hope I have not missed anyone): 1,1,1,2 (this can be read as {1,2,3},{4}) 1,1,2,1 1,2,1,1 2,1,1,1 1,1,2,2 1,2,1,2 1,2,2,1 The partitions 1,1,2,2 and 2,2,1,1 are the same and are therefore not two unique

The virtio_ring header file uses the struct device without a forward declaration. Signed-off-by: Shunsuke Mie <mie at igel.co.jp> --- include/linux/virtio_ring.h | 1 + 1 file changed, 1 insertion(+) diff --git a/include/linux/virtio_ring.h b/include/linux/virtio_ring.h index 8b95b69ef694..77a9c2f52919 100644 --- a/include/linux/virtio_ring.h +++ b/include/linux/virtio_ring.h @@ -58,6

I use model comparison with glms without mixed effects with anova(modelA,modelB), with mixed effects glm (glmmPQL), this doesn't work. Is there a way to compare model fits with glmmPQL's? Paula M. den Hartog Behavioural Biology Institute of Biology Leiden Leiden University [[alternative HTML version deleted]]

On Mon, Apr 10, 2023 at 08:00:33AM -0400, Michael S. Tsirkin wrote: > On Mon, Apr 10, 2023 at 08:28:45PM +0900, Shunsuke Mie wrote: > > Fix the build dependency for virtio_test. The virtio_ring that is used from > > the test requires container_of_const(). Change to use container_of.h kernel > > header directly and adapt related codes. > > > > Signed-off-by:

Hi, I use samba to share some files between some clients. I would like than these files are always view with -rw-rw-rw- rights by all clients event if someone change file access permissions. Can I configure smb.conf to do that ? Thanks for your help. Stephane

Hi all :-) I''m experiencing strange document loss with Ferret 0.11.4 / Ruby 1.8.4. The index handles around 3000 documents, and sometimes during a query_update operation, the delete is done but not the add. That''s very random and happens after a lot of queries (many K). I''ve tried to look into Ferret''s C code, everything looks like going well, the add *LOOKS

Hello, It's difficult to help without a sample of the format. Can you provide a short sample like 10 lines and a few columns.? Best regards, Jeremie On Wed, Jan 10, 2018 at 6:12 PM, John Kane via R-help <r-help at r-project.org> wrote: > Have a look at http://stackoverflow.com/questions/5963269/how-to-make- > a-great-r-reproducible-example > and >

The virtio_ring header file uses the struct device without a forward declaration. Signed-off-by: Shunsuke Mie <mie at igel.co.jp> --- Changes from v2: https://lore.kernel.org/virtualization/20230410074929-mutt-send-email-mst at kernel.org/ - Fix a typo of commit title include/linux/virtio_ring.h | 1 + 1 file changed, 1 insertion(+) diff --git a/include/linux/virtio_ring.h

Full_Name: Jelle Goeman Version: 1.9.0 OS: mingw32, windows 2000 Submission from: (NULL) (145.88.209.33) Hi Everyone, I get very strange results using the function hasArg with the ... function argument. In my own function: > gt <- globaltest(X,Y) > sampling(gt) works fine, but > sampling(globaltest(X,Y)) results in: Error in eval(expr, envir, enclos) : "missing"

Hello, I didn't get any response on this before, leading me to believe I've missed something fundamental. Can anybody guide me in the correct direction for more help on this? Paul ================================================= I want to be able to create a design matrix with two factors. For instance, if I have: > t1 <- factor(c(1,1,2,2)); > t2 <-

Dear All, We recently reinstalled our computing cluster. We were using CentOS 5.3 (32 bits). It is now CentOS 6.3 (64 bits), installed from the CentOS 6.2 x64 CD, then upgraded to 6.3. We have some issues with the memory needs of our running jobs. They require much more than before, it may be due to the switch from 32 to 64 bits, but to me this cannot explain the whole difference. Here are

Dear all, I have a questions regarding using the c function rmultinom.c. I got the following error message "rbinom: probability sum should be 1, but is 0.999264" Which is thrown by: if(fabs((double)(p_tot - 1.)) > 1e-7) MATHLIB_ERROR(_("rbinom: probability sum should be 1, but is %g"), (double) p_tot); I understand my probabilities do not sum to one close enough. I

Hello, I have a problem with the point'n'print driver installation. I followed the instructions on https://www.samba.org/samba/docs/man/Samba-HOWTO-Collection/classicalprinting.html#id2626650 Now if I try to print on the Windows client I get an error 0x00000006. Here is my configuration. Printer: Kyocera FS-C2026MFP+ Server: Debian Squeeze, Samba version 3.5.6, CUPS cersion 1.4.4 Client: