similar to: what is wrong with this dataset?

Dear all, I'm trying to use the "optim" function to replicate the results from the "glm" using an example from the help page of "glm", but I could not get the "optim" function to work. Would you please point out where I did wrong? Thanks a lot. The following is the code: # Step 1: fit the glm clotting <- data.frame( u =

Hello, silly question I suppose, but somehow I can't manage to extract the probabilities from a glm.predict() result: > str(res) Named num [1:9] 0.00814 0.01877 0.025 0.02941 0.03563 ... - attr(*, "names")= chr [1:9] "1" "2" "3" "4" ... I got from: # A Gamma example, from McCullagh & Nelder (1989, pp. 300-2) clotting <-

Hello, I am trying to figure out how to plot the profile likelihood curve of a GLM parameter with 95% pCI's on the same plot. The example I have been trying with is below. The plots I am getting are not the likelihood curves that I was expecting. The y-axis of the plots is tau and I would like that axis to be the likelihood so that I have a curve that maxes at the parameter estimate. I am

Hi All, I get the following error when i perform a t.test. studentt<-apply(tlr, 1, function(x) t.test(x[2:6], x[7:11],var.equal=TRUE)$p.value) # tlr is a table of 11 columns and 22500 rows. I am not able to post tlr due to its size. Error in if (stderr < 10 * .Machine$double.eps * max(abs(mx), abs(my))) stop("data are essentially constant") : missing value where

Hello, There's a long-standing bug on nouveau (this is a sample bug, but the issue has been around for a while: https://bugs.freedesktop.org/show_bug.cgi?id=85465) whereby we attempt to use the second PCOPY engine on GF116, and it is sometimes does nothing, despite mmio register 22500 saying that it's not disabled (0x22500 == 0 for this user). In the bug you can see a dump from

I have several corrupted ogg files; I believe they were created several years ago with grip, and that the corruption is not disk corruption but just something wrong with the encoder at the time. I would rather not have to re-rip these discs. The error I see with ogginfo is Warning: Hole in data (18000 bytes) found at approximate offset 702086827724505088 bytes. Corrupted ogg. Warning: Hole

Hi, I'm stuck on creating a plot with x tick labels for every Nth tick mark - how is that done? I don't see a simple solution to this in help(plot) or help(par) and what I've tried is not working, eg, the following does not work, although it seems intuitive to me that it should work: x <- seq(-100,1000,25) y <- x * x % find all the x values that are multiples of 100 tmp <-

On 10 September 2010 14:37, Marcus Daniels <mdaniels at lanl.gov> wrote: >  Hello, > > With GCC, it is possible to compile GIMPLE from an object file back to > assembly.  With Dragonegg, it seems to be a NOP.  (See below for a > comparison.) > > This appears it could be a nice way to mix the strengths of GCC and > LLVM.  If it worked..  Should it? I can't

Hi R users, I want to use dglm Package. I run the examples and it give me an error: Error en dglm(lot1 ~ log(u), ~1, data = clotting, family = Gamma) : no se pudo encontrar la funci?n "as.family" dglm can't find "as.family" function why ? Thank you for your help

On Fri, Nov 21, 2014 at 1:16 AM, Andy Ritger <aritger at nvidia.com> wrote: > Hi Ilia, > > Actually 0x90b8 is different than copy engine. I'm not very familiar > with it, but 0x90b8 is an engine for performing LZO decompression as > part of performing the copy. It has a variety of limitations (e.g., > cannot handle blocklinear format), and was only in a few Fermi

Hi, I have a data set (weight) that does not follow the Gaussian (Normal) distribution. However, I have to transform the data before applying the Gaussian distribution. I used this syntax and used log(weight) as: posJy.model<-glm(log(weight) ~ factor(pos), family=gaussian(link='identity'), subset=Soil=="Jy"). This syntax COULD NOT transform the data. But if I transform the

Hello, With GCC, it is possible to compile GIMPLE from an object file back to assembly. With Dragonegg, it seems to be a NOP. (See below for a comparison.) This appears it could be a nice way to mix the strengths of GCC and LLVM. If it worked.. Should it? Thanks, Marcus [mdaniels at dn002 dragonegg]$ cat hello.c #include <stdio.h> int hello () { printf ("Hello

Hi All, This question may not belong here, but I asked this on the R-SIG-FINANCE list and so far have not got any reply, so was hoping someone here may help. I have a basic query about TSRV and was hoping you all can shed some light on the issue. I have 22500 records for each day. So if I take the slow frequency as 5 mins I can have 75 subgrids. An estimate of the 5-minute realized volatility

Hi Ilia, Actually 0x90b8 is different than copy engine. I'm not very familiar with it, but 0x90b8 is an engine for performing LZO decompression as part of performing the copy. It has a variety of limitations (e.g., cannot handle blocklinear format), and was only in a few Fermi chips, as I understand it. It is probably easiest to just ignore it. You can distinguish this decompress engine

If anyone is interested, I smushed the less than 0 quality rate setting into oggdrop. The diff is below. Alan <p><p>Index: encthread.c =================================================================== RCS file: /usr/local/cvsroot/win32-tools/oggdrop/encthread.c,v retrieving revision 1.9 diff -r1.9 encthread.c 117,118c117,118 < if (quality < 0) < quality = 1; ---

On Fri, Nov 21, 2014 at 01:39:55AM -0500, Ilia Mirkin wrote: > On Fri, Nov 21, 2014 at 1:16 AM, Andy Ritger <aritger at nvidia.com> wrote: > > Hi Ilia, > > > > Actually 0x90b8 is different than copy engine. I'm not very familiar > > with it, but 0x90b8 is an engine for performing LZO decompression as > > part of performing the copy. It has a variety of

I am using CALDERA OPENLINUX 3.1.1 with latest patches and Samba 2.2.5 on it. But recently I start to receive linux error "kernel: VFS MAX FILES REACHED" and Samba error "Too many open files". After it I tuning linux increase open files descriptors in /etc/security/limits.conf * soft nproc 4094 * hard nproc 16384 * soft

r-help at lists.R-project.org ? Hello ? First I want to implement exponential regression in R, with out constant for the following formula. Y = exp(a*X) ?a? is coefficient I wanted to determine. That I could do also in SPSS but my question is rather to estimate the ?standard error of ??i ?at each Xi. This is called in SPSS ?satndard error of mean prediction? or generally known for non-linear

I compiled linux-2.6 from git. I think my last pull/merge? is this one (git show) commit d03ab2d78b6ab62e94f9958da50b4419c27e0f60 Author: Marcin Ko?cielnicki <koriakin at 0x04.net> Date: Mon Mar 1 00:18:39 2010 +0000 drm/nv50: Improve PGRAPH interrupt handling. I also compiled nouveau-drm_2.4.18-1_amd64.deb and nouveau_drv.so that day and installed them to a Debian system. I'm

Hello... Im trying to setup HTB to allow me to shape traffic from two upstreams that meets on single lan0 interface. I prefer to use HTB qdisc within HTB root qdisc for cleaner rules design. Seems that it doesnt work at all. tc -s class show doesnt show any traffic on other classes attached to HTB qdisc. Linux 2.6.20.7 iproute-2.6.20-070313 Weird thing is that tc -s class show that 1: and 2: