similar to: predict glm() with offset

Dear R friends. I´m trying to fit a Logistic Regression using glm( family='binomial'). Here is the model: *model<-glm(f_ocur~altitud+UTM_X+UTM_Y+j_sin+j_cos+temp_res+pp, offset=(log(1/off)), data=mydata, family='binomial')* mydata has 76820 observations. The response variable f_ocur) is a 0-1. This data is a SAMPLE of a bigger dataset, so the idea of setting the offset is to

Dear R friends I´m interested into apply a Jackknife analysis to in order to quantify the uncertainty of my coefficients estimated by the logistic regression. I´m using a glm(family=’binomial’) because my independent variable is in 0 - 1 format. My dataset has 76000 obs, and I´m using 7 independent variables plus an offset. The idea involves to split the data in let’s say 5 random subsets and

Hello Trying to develop variograms and kriged surfaces from a point file. Here is what I've done so far. library(gstat) # also loads library(sp) library(lattice) soilpts$x <- soilpts$UTM_X soilpts$y <- soilpts$UTM_Y soil.dat <- subset(soilpts, select=c(x, y, Area, BulkDensity, LOI, TP, TN, TC, Total_Mg)) dim(soil.dat) [1] 1292 7 coordinates(soil.dat) <- ~ x+y

Hi all, I have an Active Record problem. Having this classes declaration. class Sale < ActiveRecord::Base has_many :sale_line_items, :class_name ''SaleLineItem'' ... end class SaleLineItem < ActiveRecord::Base ... end When I try to delete an item in a Sale class method from sale_line_items (sale_line_items.delete(li)), it extracts ''li'' from

Hello List! I''m glad to announce the update of Rails Plugin Fossilize (0.3.0): Fossilize ========= Allow you distribute your rails application as a Gem. It add a "fossilize" generator and a package task to rake. What it offer ============= Fossilize offer you two ways of using your distributable rails application: * a ''unpack'' method (by default):

Hi, How can the volume serial number of a share can be especified in the smb.conf? I have a machine that each times it connects to the share displays two differents volume serial number, not at the same time.... The server uses samba 3.2.5-4 over debian and the disks are in a software raid 1.

Hola a todos! Estoy intentando crear un vector (alt) a partir de la repetición de valores provenientes de otro vector (altitud). A cada valor de altitud lo quiero repetir 1247 veces, y de ahi continuar con el siguiente valor de altitud. Probé varias cosas, pero esto me pareció lo más coherente: altitud=read.csv("C:/Users/IER/Dropbox/Pasantia Castelar/YungasLigustroTS/altitud44.csv")

Dear r-list: I am using rpart to build a tree on a dataset. First I obtain a perhaps too large tree: > arbol.bsvg.02 <- rpart(formula, data = bsvg, subset=grp.entr, control=rpart.control(cp=0.001)) > arbol.bsvg.02 n= 100000 node), split, n, loss, yval, (yprob) * denotes terminal node 1) root 100000 6657 0 (0.93343000 0.06657000) 2) meses_antiguedad_svg>=10.5 73899 3658

#dsync -D -vv -u stak at nti.uff.br mirror maildir:/CX4/vdomains/ nti.uff.br/stak/mail dsync(root): Fatal: Mail locations must use the same virtual mailbox hierarchy separator (specify separator for the default namespace) My mail.conf *mail_location = mbox:/CX4/vdomains/%d/home/%n/mail:INBOX=/CX4/vdomains/%d/var/spool/mail/%n:INDEX=/tmp/imap-index/%d/%n *What that error means ? Before when i

Hello to everyone. I´m writing you because I´m feeling a bit frustrated with my work. My work consists in finding the relation between the amount of fires and the weather, so, my response variable is the amount of fires in a fire season and the explanatory variables are the temperature, the amount of precipitation and the some others…. my problem is this; I keep getting the wrong sign in the

Hi, If I have a DSL router with 300kpbs upload bandwith: Linux(eth0) ---(100Mbps)---> DSL ---(300kbps)---> INTERNET Which is the recomended value I must set in my QoS setup for the maximum outbound bandwith of my ethernet interface ? The max of 300kbps or something less like 290kbps ? Especify less than the maximum available bandwith is advisable to ensure that the QoS algoritms goes

Hi, I couldn''t find a better word for pop, but let''s go: There are 3 models: sale, pop and product sale: -total:integer pop: -quantity:integer -discount:float product: -name:string -price:float A sale can have many different products, each one has a quantity and discount associated. So each sale has_many :products through => pops And each product has_many :sales

Hi R community I need to interpolate precipitation data for a natural park. I have precipitation data from some climate stationts. (I know the table is not complete but I only need to show you X,Y, Altitude and PrepJul) X Y Altitude PrepJan PrepFeb PrepMar PrepAp PrepMay PrepJun PrepJul 597706 4093438 41 0 9 77,8 63,1 17,5 0 2,6 597535

I've been using liveice to feed an encoded stream to an icecast broadcast server. Now I want to feed two streams at the same time to the broadcast server. To acomplish that I am using a dual soundcard setup in a debian linux server. I've setup two different mountpoints in the icecast server wich will receive the stream and in turn will broadcast the streams to the listeners. This is, the

Hello to everyone. I'm trying to use the %in% to match to vectors in ff format. a<-as.ff(data[,1]) %in% fire$fecha > aff (open) logical length=3653 (3653) [1] [2] [3] [4] [5] [6] [7] [8] [3646] FALSE FALSE FALSE FALSE FALSE FALSE FALSE FALSE : FALSE [3647] [3648] [3649] [3650] [3651] [3652] [3653] FALSE FALSE FALSE FALSE FALSE FALSE

Hi R-experts, I am using the envfit function over an ordination of floristic data. The problem is that every time that I run it changes the results. Sometimes dramatically, selecting variables that the first time were not significant. I do not get what could be the problem or if is normal given the permutations are different. # the NMDS ordination gap_flor_NMDS_chord <- metaMDS(gaps_flor,

FWIW, majordomo bounces messages with a subject line of 'help'. ----- Forwarded message from owner-icecast@xiph.org ----- From: "Moises Soto" <glester@hotmail.com> To: <icecast@xiph.org> Subject: Help with Date: Wed, 9 Jan 2002 11:21:54 -0400 I've been using liveice to feed an encoded stream to an icecast broadcast server. Now I want to feed two streams at the

Hello to everyone. In MATLAB you can find a function that executes a text string. I have not find it in R. Does anyone knows a similiar function in R, i think a valuable tool. Thank you. Lucas. [[alternative HTML version deleted]]

Hello to everyone. I need your help. I´m trying to fit the same *glm.nb* to a different data set and i am getting these errors in some of the data. Sometimes, one data set has two of these errors when fitting the model. 1.- Error en while ((it <- it + 1) < limit && abs(del) > eps) { : valor ausente donde TRUE/FALSE es necesario 2.- Mensajes de aviso perdidos 1: In sqrt(1/i)

Dear R people. I´m facing a big problem. I need to create a matrix with 10.000 columns and 750.000 rows. matrix<- as.data.frame(matrix(data=0L, nrow=750000, ncol=10000) as you can see, the data frame has huge dimesions. I was able to find out about thr "L" in data, this way I´m telling that my data is "integer" class. By doing this object.size(matrix) tells me that this is