similar to: body bug?

When fitting model 'y~1', the 'lm' response is ok, but it is not possible to get predictions using 'predict'. May be this is a bug, at least i think it is somehow inconsistent. > rm(list=ls()) > d <- data.frame(x=runif(50),y=rnorm(50)) > plot(d) > d.lm <- lm(y~1,data=d) > d.lm Call: lm(formula = y ~ 1, data = d) Coefficients: (Intercept)

On Tue, 25 Sep 2001 12:27:27 GMT, Jorge Luis Ojeda Cabrera <jojeda@posta.unizar.es> wrote in message <20010925.12272700@macdesp8.unizar.es>: >When fitting model 'y~1', the 'lm' response is ok, but it is not possible >to get predictions using 'predict'. May be this is a bug, at least i >think it is somehow >inconsistent. It still appears in 1.3.1,

There's more to it than that (which was already fixed in R-devel). The internal code expects there to be columns in the model frame. I have already fixed this in R-devel (by faking an unused column). On Tue, 25 Sep 2001 murdoch@stats.uwo.ca wrote: > On Tue, 25 Sep 2001 12:27:27 GMT, Jorge Luis Ojeda Cabrera > <jojeda@posta.unizar.es> wrote in message >

I 've been working with R since last year and I am quite "happy" with R capabilities. Nevertheless, I need to present my job to other people (S-users) and I am facing a problem with R & S compatibility: my R-scripts make intensive use of functions local scope, and as far as I know(I am a new to S) S does not have this capability. I use to do things like: create.f <-

I would like you consider that the function ecdf could be extended in the following way to handle weights when computing Empirical distribution Functions. There exist particular cases that supports this kind of extension, see for example: Rao, C. R., 1997. Statistic and True. Putting chance to work. World Scientific Publishing. Cox, D. R., 1969. Some Sampling Problems in Technology. New

Hi, I am using fanny and I have estrange results. I am wondering if someone out there can help me understand why this happens. First of all in most of my tries, it gives me a result in which each object has equal membership in all clusters. I have read that that means "the clustering is entirely fuzzy". Looking at the graphics it is really difficult to understand how objects with so

Hello, I'm suffering a estrange problem in a WinNT-Samba3 environment. I have two servers: WinNT4 (PDC of domain A-DOMAIN) and Samba3 (PDC of B-DOMAIN). A-DOMAIN and B-DOMAIN trust each other (I had followed the procedures described in HOWTO Chapter 16 successfully). The problem arises when I assign permission in WinNT server's folders (A-DOMAIN) for users in the Samba domain (B-DOMAIN).

Hello, I am running a i7 MacBook Pro 2011. If I write: @g = global double 0.000000e+00 define i32 @main() { entry: %0 = load double* @g %1 = fmul double 1.000000e+06, %0 store double %1, double* @g ret i32 0 } in test.ll and I run > llc test.ll > gcc test.s I get: test.s:12:no such instruction: `vmovsd _g(%rip), %xmm0' test.s:13:no such instruction: `vmulsd LCPI0_0(%rip),

I have mysql running on a different machine. I''ve editted database.yml, here''s the entry development: adapter: mysql database: cookbook username: root password: mypassword host: posta.rmasonfamily.info test: adapter: mysql database: cookbook username: root password: mypassword host: posta.rmasonfamily.info production: adapter: mysql database: cookbook

Hola, Dado que hay varios matemáticos en el foro, ?alguien puede recomendar algún buen libro en español o ingles que ayude a entender el álgebra lineal de primero de matemáticas ?. Gracias y un saludo JM Arbones

Dear list users, let me ask you this trivial question. I worked on that for a long time, by now. Suppose to have a data frame with NAs and to sum some columns with rowSums: df <- data.frame(A = runif(10), B = runif(10), C = rnorm(10)) df[1, ] <- NA rowSums(df[ , which(names(df) %in% c("A","B"))], na.rm=T) If all the elements of the selected columns are NA, rowSums

Hello, This is a little off-topic. But I am writing a compiler to llvm ir for a language that admits recursive nested functions and am stuck as to how to translate them. Concretely, I'm trying to lift them all to the topmost level and pass all their free variables explicitly as arguments. To do this, I have to determine all their free variables in their bodies. In particular when I come

Desde el terminal de Ubuntu, he conseguido descargar los datos usando esto: DESCARGAR LISTA DE ESTACIONES wget --method GET --header 'cache-control: no-cache' --no-check-certificate --output-document - https://opendata.aemet.es/opendata/api/valores/climatologicos/inventarioestaciones/todasestaciones/?api_key=*vuestra api key* DESCARGAR DATOS HISTÓRICOS Diciembre 2017, ALICANTE wget

Dear R-users, suppose I have three dataframes like these df1: mydate min_temp 31032013 12 01042013 8 02042013 -999 df2: mydate min_temp 31032013 10 01042013 11 02042013 14 df3: mydate min_temp 31032013 4 01042013 3 02042013 5 where -999 means that the temperature data is not available (at the moment I cannot change it to NA because I am not the db administrator); suppose also that oggi is

Hi, Compiling a simple function with the 3.0 SVN branch llc produces .s output containing cfi_startproc, cfi_endproc, cfi_def_cfa_offset that gcc refuses to compile (OS X 10.7 on X86-64) with the following error: test.s:6:Unknown pseudo-op: .cfi_startproc test.s:10:Unknown pseudo-op: .cfi_def_cfa_offset test.s:10:Rest of line ignored. 1st junk character valued 51 (3). test.s:38:Unknown

Dear list, I want to add SJLJ exception handling to my frontend. Unfortunately, there doesn't seem to be any examples in the documentation as to how to use the intrinsics @llvm.eh.sjlj.setjmp @llvm.eh.sjlj.longjmp @llvm.eh.sjlj.lsda @llvm.eh.sjlj.callsite Is there a way to force Clang to use SJLJ exception handling for C++? That way I would be able to look at its output to learn how to use

Hi, When calling a function, does the llvm code generator support passing structured types (arrays, structs, etc.) by _value_? I wrote some small examples, and it seemed to work, but I was wondering if anything can go wrong if the structured types are very large... Thanks, N

Hi, Is this correct behaviour? test.ll: declare {} @__ex__print_int(i64) define i32 @main() { entry: %0 = call i64 @f.1() %1 = call {} @__ex__print_int(i64 %0) ret i32 0 } define internal fastcc i64 @f.1() { entry: ret i64 7 } > opt -std-compile-opts test.ll -S ; ModuleID = 'test.ll' define i32 @main() noreturn nounwind { entry: tail call void @llvm.trap()

Dear, I need to create a custom device extension in order to dial a local queue. Suppose my queue number is 8888, how can fill the Dial field from the custom extension ??? Because if I put just 8888 or Local/8888, I don't succeed. Thanks a lot

Dear R-list users, I need to use strptime because I have to deal with date with hours and minutes. I read the manual for strptime and I also looked at many examples, but when I try to apply it to my code, I always encounter some problems. I try to change the default format, with no success. Why? How can I change the format? 1. init_day <- as.factor("2015-02-24-00-30")