similar to: Inconsistencies (PR#550)

>>>>> "BDR" == Prof Brian D Ripley <ripley@stats.ox.ac.uk> writes: BDR> On Mon, 22 May 2000 berwin@maths.uwa.edu.au wrote: >> Dear all, >> >> I was playing around with some C-code that I dynamically linked to R >> and noticed the following inconsistency: > > dat <- matrix(1,3,3) > >

Dear all, I've tried a multiple regression, and now I want to try a cross-validation. I obtain this error (it must be sth related to df) that I don't understand, any help would be appreciated. cv.lm(df= dat, lm2.52f, m=3) Error en `[.data.frame`(df, , ynam) : undefined columns selected lm2.52f is my lm object, dat is a dataframe where the variables involved in .lm are I tried CVlm

I am trying to export a dateframe created in R: > Duration_summary V1 2.5 % 97.5 % V4 2.5 % 97.5 % [1,] 1 0.46076018 1.128776 1.000000 0.5280828 0.9576338 [2,] 0 0.00000000 0.000000 1.000000 0.1741793 1.2352705 [3,] 1 0.46566719 1.313711 1.000000 0.7233312 1.4097987 [4,] 1 0.38866371 2.453226 0.976024 0.6377314 1.3493957 [5,] 1 0.08894066 1.036830

.... I have a 64-bit CentOS 5.6 VM running on a 64-bit FC14 host. I notice that it has an automatic logout feature after a period of idleness, see below: [wam at centos-5, CFD, 5:52:45pm] 1183 % ll -tr LIST LIST.Update /usr/local/bin/P* /usr/local/bin/*WAM* ; ps -u wam -lf ; uname -a ; date ls: No match. F S UID PID PPID C PRI NI P SZ:RSS WCHAN STIME TTY TIME CMD 5 S

There's something peculiar that I do not understand here. However, did you realize that the thing you are assigning into parts of `a' is NULL? Check you're my.test.boot.ci.1: It's NULL. Be that as it may, I get: > a <- data.frame(matrix(1:4, nrow=2), X3=NA, X4=NA) > a X1 X2 X3 X4 1 1 3 NA NA 2 2 4 NA NA > a[a$X1 == 1,]$X3 <- NULL > a X1 X2 X3 X4 1 1

Dear all, while looking at some R-code submitted by students in a unit that I teach, I came across constructs that I thought would lead to an error. Much to my surprise, the code is actually executed. A boiled down version of the code is the following: > tt <- function(x, i){ + mean(x[i,2])/mean(x[i,1]) + } > dat <- matrix(rnorm(200), ncol=2) > mean(dat[,2])/mean(dat[,1]) [1]

Hi R-people: When I use the command to fit a model with an intercept, only: glm ( formula=haspdata ~ 1, data=dat, family=binomial, weights= dat$hy.wgt.s, subset=(dat$haspdat0!=3) ) I get the message: Warning message: non-integer #successes in a binomial glm! in: eval(expr, envir, enclos) Does anyone know what this means?? The data for this command is listed below. Thanks, Phil Smith CDC

Dear useR, I have a dataset like this below, > prevRND.dat <- read.table("C:\\workdir\\prevRND.txt", + header=FALSE, # No header. + col.names = c("X","Y","Z"), + sep = ",") > prevRND.dat X Y Z 1 A A 0.950933 2 A B 0.143600 3 A C 0.956133 4 B A

I have contributed a "gam" library to CRAN, which implements "Generalized Additive Models". This implementation follows closely the description in the GAM chapter 7 of the "white" book "Statistical Models in S" (Chambers & Hastie (eds), 1992, Wadsworth), as well as the philosophy in "Generalized Additive Models" (Hastie & Tibshirani 1990,

I have contributed a "gam" library to CRAN, which implements "Generalized Additive Models". This implementation follows closely the description in the GAM chapter 7 of the "white" book "Statistical Models in S" (Chambers & Hastie (eds), 1992, Wadsworth), as well as the philosophy in "Generalized Additive Models" (Hastie & Tibshirani 1990,

Following simple program gives different answers when compiling with GCC (4.8.4) and ToT Clang: $ cat builtin_fmodf_bugpoint.cpp #include <cstdio> int main(int argc, char** argv) { const float a = 1.0f; const float b = 0.1f; printf("%f mod %f = %f\n", a, b, __builtin_fmodf(a, b)); return 0; } $ g++ -o builtin_fmodf_bugpoint_gcc builtin_fmodf_bugpoint.cpp $

Dear UseRs called Hadley, or Paul, I am trying to print an edited ggplot2/grid graphics to a metafile. With the commented line below it works, but when I edit the plot by uncommenting the line, it fails, because it's illegal to have 2 graphics in a metafile. It works with pdf, but even then I get two plots, which is a nuisance. I found a workaround by using windows(); savePlot, but it only

Hi guys, I must be missing something obvious with pessimistic locking. Here is what I do (Rails 2.3, mySQL): In my Foo model: def test Foo.transaction do lock = UpdateLockForFoo.find_by_foo_id(self.id, :lock => true) sleep 30 end end For each Foo model, I have an UpdateLockForFoo. This "lock" needs to be acquired before doing any changes to Foo. I believe it works,

Hi Peter, I got all your patches and applied them, and I wanted to say kudos for a great job done. This is the output of upsc that I'm getting right now. I had the same issue when using any cvs version (from two weeks ago) that made my UPS go into a state where it wasn't functioning normally (checking the kernel log shows that it cannot be connected and I had around a page of DEVFS_USB

Hi! I was trying to implement a probit model on a dichotomous outcome variable and found that the predictions were outside the (0,1) interval that one should get. I later tried it with some simulated data with a similar result. Here is a toy program I wrote and I cant figure why I should be getting such odd predictions. x1<-rnorm(1000) x2<-rnorm(1000) x3<-rnorm(1000)

Dear colleagues Please could someone kindly explain the following inconsistencies I've discovered when performing logical calculations in R: 8.8 - 7.8 > 1 > TRUE 8.3 - 7.3 > 1 > TRUE Thank you, Emma Jane [[alternative HTML version deleted]]

Hello all R-users! In Jun-2005, I find the follow discussion about using of LME function ( in NLME library ) for fitting non-replicate data The thread: ANOVA vs REML approach to variance component estimation http://tolstoy.newcastle.edu.au/R/help/05/06/6498.html Someone expose the follow problem: # non-replicate data y <- c(2.2, -1.4, -0.5, -0.3, -2.1, 1.5, 1.3, -0.3, 0.5, -1.4,

Is anyone working on Vorbis over RTP? It's been silent on this frontier for a while. I have been hearing from representatives of the IETF that they were not impressed by the previous draft, and that they were not coherent with RTP thinking and other RTP transport specs for various content. There exist a textbook on RTP which was released recently, if someone wants to learn it from the ground

Hi, On 20-11-15 17:07, Samuel Pitoiset wrote: > > > On 11/20/2015 11:36 AM, Hans de Goede wrote: >> Hi Samual, et al, > > Hi Hans, > >> >> In >> http://cgit.freedesktop.org/mesa/mesa/commit/src/gallium/drivers/nouveau?id=ff72440b40211326eda118232fabd53965410afd >> >> you write: "This compute support has been tested by >> Pierre

Hi- I need to subset the following data by the column 'dal' for values that equal the regular interval seq(0, 150, by=0.5) exactly.... ....excluding rows with irregular 'dal' values such as c(2.888958, 2.891620), etc. data<-data.frame(id=id, dal=dal, date=date, mu.x=mu.x) $dal [1] 0.000000 0.500000 1.000000 1.500000 2.000000 2.500000 2.888958 2.891620 3.000000 3.245405