similar to: problems with glm

Hi Stephen, Check the help for predict.glm(). The argument for passing new data is actually 'newdata', as in: > pred = predict(glm.model, newdata=form[150001:200000,-1], > type="response") Cheers Joe -----Original Message----- From: r-help-bounces at r-project.org [mailto:r-help-bounces at r-project.org] On Behalf Of stephenc at ics.mq.edu.au Sent: Tuesday, 16

I am trying to train on part of my data and test on another part: > glm.model = glm(as.factor(h_finished) ~ . , family=binomial, data=form[1:150000,]) > pred = predict(glm.model, data=form[150001:200000,-1], type="response") > t = table(pred, form[150001:200000,1]) Error in table(pred, form[150001:2e+05, 1]) : all arguments must have the same length but try as I might

Hello. My setup: _ one Samba 3.5 domain (XXXXXXXX), with a PDC and a BDC, both running FreeBSD; _ one AD domain (YYYYYYYY) running on two Windows 2003 DCs; _ bidirectional trust between the two domains. Everything used to work until I moved the PDC from Samba 3.5 (EOL'ed) to 3.6; now, users from domain YYYYYYYY cannot access the PDC's shares. I used to have in smb.conf: >

Hi. I introduced "is_partially_uptodate" aops for OCFS2. A page can have multiple buffers and even if a page is not uptodate, some buffers can be uptodate on pagesize != blocksize environment. This aops checks that all buffers which correspond to a part of a file that we want to read are uptodate. If so, we do not have to issue actual read IO to HDD even if a page is not uptodate

dear all, i've looked at the r-intro (chapter 10, writing your own functions) and searched the r-help archives but am still stuck at the following. i have a simple function, something like: myhist<-function(yvar) { y<-subset(myframe,yvar>1 & yvar<=150000,select=yvar) attach(y) hist(yvar) } calling it as follows: myhist(x1) gives the following error: Error in

Hello all, I have setup my gentoo samba server as a domain member server of a windows 2003 Active Directory Domain using winbind. I am able to see all the users and groups with wbinfo -u and wbinfo -g andwbinfo -t does not report any errors. 'Getent passwd' and 'getent group' return local and domain users and groups. My problem is that when I login as a domain user, my prompt

How do I convert a list to a matrix? --8<---------------cut here---------------start------------->8--- list(c(50000, 101), c(1e+05, 46), c(150000, 31), c(2e+05, 17), c(250000, 19), c(3e+05, 11), c(350000, 12), c(4e+05, 25), c(450000, 19), c(5e+05, 16)) as.matrix(a) [,1] [1,] Numeric,2 [2,] Numeric,2 [3,] Numeric,2 [4,] Numeric,2 [5,] Numeric,2 [6,] Numeric,2 [7,]

I am trying to write a wrapper for the last example in help(cv.glm) that deals with leave-one-out-cross-validation (LOOCV) for a logistic model. This wrapper will be used as part of a bigger program. Here is my wrapper funtion : logistic.LOOCV.err <- function( formu=NULL, data=NULL ){ cost.fn <- function(cl, pred) mean( abs(cl-pred) > 0.5 ) glmfit <- glm(

Dear all, I've been having problems running update() to re-fit a cpglm model inside a function (as in the code below). The solution is probably simple, but I'm stuck. If anyone could help, I'd greatly appreciate it. Regards, Rubem ## R code library(cplm) ## Data simulation period<-factor(1:4)                        herd<-factor(1:50)  

Hi, I have been trying to join a Samba Domain member server to the AD and use LDAP for IDMAP storage. I have run into many strange issues and I was hoping someone can please take time to clarify things for me. I have read quite a bit (I own both the Samba books by Terpstra) and done a lot of Google searching. I think part of my problem is the unusual setup I have, as all the examples in the

I'm trying to setup Samba in ADS security mode so I can run winbind for NSS and Kerberos for user authentication, chiefly for shell accounts for developers. These hosts will not provide any file or printer services, at least in the near-term. My hosts are CentOS 3 (a free RHEL3 clone) and my ADS servers are Windows 2000 (not 2003), in hybid mode. I am using stock RPMs for both Kerberos and

I am trying to use winbind on a domain with 150000 users, but I can't make it work. Is that because I have done something wrong, or doesn't winbind scale to that many users? It works a little, I can use wbinfo to translate between sid, uid and gid, and I can also get the groups a user is in. But when I try to get the list of all users with 'wbinfo -u', I get a timeout and winbindd

Hi, 1. How do I construct 95% prediction interval for new x values, for example - x = 30000? 2. How do I construct 95% confidence interval? my dataframe is as follows : >dt structure(list(y = c(26100000, 60500000, 16200000, 30700000, 70100000, 57700000, 46700000, 8600000, 10000000, 61800000, 30200000, 52200000, 71900000, 55000000, 12700000 ), x = c(108000, 136000,

Dear r-users,   I have saved data in word pad.  I would like to extract certain part of data only, for example   19710629 080000(PARTIAL) 39.3 at interval beginning 19701126 010326 19720629 080000(PARTIAL) 33.8 at interval beginning 19720517 144507 19730629 080000(PARTIAL) 32.2 at interval beginning 19720910 135747   and so on...   The original data set are given below:      

Dear R users, I just join this forum a few minutes ago. My friend recommend this forum to me.   I have data in excel csv. My problem is I would like to extract some number, for example as highlighted below;   Year Ending           Maximum Value (Rain mm)  Gap from 19710509 090000 to 19710607 080000 of   28.96 Days 19710629 080000(PARTIAL)       39.3 at interval beginning 19701126 010326  Gap from

Kelly Lesperance wrote: > I did some additional testing - I stopped Kafka on the host, and kicked > off a disk check, and it ran at the expected speed overnight. I started > kafka this morning, and the raid check's speed immediately dropped down to > ~2000K/Sec. > > I then enabled the write-back cache on the drives (hdparm -W1 /dev/sd*). > The raid check is now running

<!DOCTYPE html PUBLIC "-//W3C//DTD HTML 4.01 Transitional//EN"> <html> <head> </head> <body bgcolor="#ffffff" text="#000000"> <font size="-1"><font face="Verdana">Hi all,<br> <br> In last year, i had made some research and benchmarks based on CentOS 4 to know which filesystem is better for

Hi, I have a little problem with my Samba PDC setup. I have a Debian / Etch with Samba 3.0.24 with OpenLDAP Backend and use it as PDC with domain logons and roaming profiles enabled. The logon for normal domain users work without any problems, but when I login with my domain administrator account my Windows XP Professional client I have a big problem. After the logout, when it tries to copy

Dear R helpers, I need to estimate a probit model with box constraints placed on several of the model parameters. I have the following two questions: 1) How are the standard errors calclulated in glm (family=binomial(link="probit")? I ran a typical probit model using the glm probit link and the nlminb function with my own coding of the loglikehood, separately. As nlminb does not

Dear all, I want to make a cross table out of a data set which is 2 columns wide and more than 150000 rows long. When I use the table() function I get an error message This is the code I have used: >Dataset <- read.table("test.txt", header=TRUE, sep=",", na.strings="NA", dec=".", strip.white=TRUE) > .T <-table(Dataset$K1,Dataset$K2) This