search for: resp

...I'm trying to fit a model using variance components in R, but if very new on it, so I'm asking for your help. I have imported the SPSS database onto R, but I don't know how to convert the commands... the SPSS commands I'm trying to convert are: VARCOMP RATING BY CHAIN SECTOR RESP ASPECT ITEM /RANDOM = CHAIN SECTOR RESP ASPECT ITEM /METHOD = MINQUE (1) /DESIGN = CHAIN SECTOR RESP ASPECT ITEM SECTOR*RESP SECTOR*ASPECT SECTOR*ITEM CHAIN*RESP CHAIN*ASPECT CHAIN*ITEM RESP*ASPECT RESP*ITEM SECTOR*RESP*ASPECT SECTOR*RESP*ITEM CHAIN*RESP*ASP...

I am running samba 3.0.20a in my FC4 linux machine, its running for two days. In /tmp direcotry I see lots of temporary files. Is it safe to remove these files? also I would like to why its creating these many files? SMBclose.8.req SMBreadX.58.resp SMBsesssetupX.61.req SMBtrans2.3.resp SMBwriteX.60.req SMBclose.8.resp SMBreadX.59.req SMBsesssetupX.61.resp SMBtrans.24.req SMBwriteX.60.resp SMBclose.9.req SMBreadX.59.resp SMBsesssetupX.62.req SMBtrans2.4.req SMBwriteX.61.req SMBclose.9.resp...

Hey, everyone, I am using proportional odds model for ordinal responses in dose-response experiments. For some samll data, SAS can successfully provide estimators of the parameters, but the built-in function polr() in R fails. Would you like to tell me how to make some change so I can use polr() to obtain the estimators? Or anyone can give me a hint about the con...

...working with some data where I need to select the first occurrence of a value after the occurrence of another value. The data has two columns, one with a time and one with occurence of certain events. The column of data I want to select from looks like this (and each of these events have a corresponding time in another column). Stat71 OutMag FirstResp InMag MagUp OutMag MagDwn Resp Resp Resp InMag MagUp OutMag InMag OutMag InMag OutMag InMag OutMag InMag MagDwn OutMag Resp MagUp InMag MagDwn OutMag Resp MagUp Using "which" I have selected all the individual occurrence...

...ssary to match up the diagnostic error messages. -- David Favro lists.samba.org at meta-dynamic dot com http://meta-dynamic.com The patch: --- socket.c.~1.116.~ 2006-02-03 13:46:38.000000000 -0500 +++ socket.c 2006-02-20 07:18:11.000000000 -0500 @@ -364,6 +364,7 @@ resp->ai_protocol); if (s == -1) { + rsyserr(FERROR, errno, "socket(%d,%d)", resp->ai_family, resp->ai_protocol ); /* See if there's another address that will work... */ continue;...

...make a trellis boxplot with the HSuccess (y-axis) in each Rayos (beach sections) (x-axis), for each Aeventexhumed (A, B, C) - nesting event. I am not able to do so and keep receiving: Error in eval(expr, envir, enclos) : object 'Rayos' not found Please advise, Jean require(plyr) resp <- read.csv("ABC Arribada R File Dec 12 Jean Jang.csv") envir <- read.csv("Responses Environ. Arribada Dec. 12.csv") resp <- resp[!is.na(resp$Aeventexhumed), ] resp$QuadratEvent <- paste(resp$QuadratID, resp$Aeventexhumed, sep="") resp$QuadratEvent &lt...

...priority field */ u8 crn; - u32 pi_bytesout; /* DataOUT PI Number of bytes */ - u32 pi_bytesin; /* DataIN PI Number of bytes */ + __virtio32 pi_bytesout; /* DataOUT PI Number of bytes */ + __virtio32 pi_bytesin; /* DataIN PI Number of bytes */ u8 cdb[VIRTIO_SCSI_CDB_SIZE]; } __packed; /* Response, followed by sense data and data-in */ struct virtio_scsi_cmd_resp { - u32 sense_len; /* Sense data length */ - u32 resid; /* Residual bytes in data buffer */ - u16 status_qualifier; /* Status qualifier */ + __virtio32 sense_len; /* Sense data length */ + __virtio32 resid; /* Residual byte...

...priority field */ u8 crn; - u32 pi_bytesout; /* DataOUT PI Number of bytes */ - u32 pi_bytesin; /* DataIN PI Number of bytes */ + __virtio32 pi_bytesout; /* DataOUT PI Number of bytes */ + __virtio32 pi_bytesin; /* DataIN PI Number of bytes */ u8 cdb[VIRTIO_SCSI_CDB_SIZE]; } __packed; /* Response, followed by sense data and data-in */ struct virtio_scsi_cmd_resp { - u32 sense_len; /* Sense data length */ - u32 resid; /* Residual bytes in data buffer */ - u16 status_qualifier; /* Status qualifier */ + __virtio32 sense_len; /* Sense data length */ + __virtio32 resid; /* Residual byte...

...riables, so I should be able to aggregate up my data frame and run glm(..., family="binomial",weights=wts) and get the same coefficients (maybe this is my incorrect assumption, but I can't see why it would be). Anyways, here's a minimum working example below: > d = data.frame( RESP=c(rep(1,5),rep(0,5)), INDEP=c(1,1,1,1,0,0,0,0,0,0) ) > glm( RESP ~ INDEP, family="binomial", data=d ) Call: glm(formula = RESP ~ INDEP, family = "binomial", data = d) Coefficients: (Intercept) INDEP -1.609 21.176 Degrees of Freedom: 9 Total (i.e. Null);...

I'd like to convert a data.frame from (long format): gen maxh resp 1 12.3 y . . . . . . where resp is a factor with levels "y" and "n" to the wide format: gen maxh.y maxh.n . . . . . . I've done it as follows: maxh.y <- split(maxh, resp)$y gy <- split(gen,...

Hi, I would have thought that these two constructions would produce the same result but they do not. Resp <- rbinom(10, 1, 0.5) Stim <- rep(0:1, 5) mm <- model.matrix(~ Stim) Xb <- mm %*% c(0, 1) ifelse(Resp, log(pnorm(Xb)), log(1 - pnorm(Xb))) pnorm(as.vector(Xb), lower.tail = Resp, log.p = TRUE) > ifelse(Resp, log(pnorm(Xb)), log(1 - pnorm(Xb))) [1] -0.6931472 -1.8410216 -0.6931472...

...ly understand how to use the function "apply" instead of a "for loop", particularly for a function with multiple entries. I have a big data file and would like to apply a function in multi thread to accelerate the processus. I have a data frame containing values of* CO2 in ppm (resp[i,6])* that I want to convert in umol of CO2 emitted by stem volume biomass (CO2v) and stem area (CO2s). (tree respiration) The loop that I have is calculating the CO2 fluxes for each row. *Script :* for (i in 1:nrow(resp)) { k=resp[i,5] CO2umol[i]<-resp[i,6]*((Press*infoch[k,11]*1e-6)/(R*(re...

...priority field */ u8 crn; - u32 pi_bytesout; /* DataOUT PI Number of bytes */ - u32 pi_bytesin; /* DataIN PI Number of bytes */ + __virtio32 pi_bytesout; /* DataOUT PI Number of bytes */ + __virtio32 pi_bytesin; /* DataIN PI Number of bytes */ u8 cdb[VIRTIO_SCSI_CDB_SIZE]; } __packed; /* Response, followed by sense data and data-in */ struct virtio_scsi_cmd_resp { - u32 sense_len; /* Sense data length */ - u32 resid; /* Residual bytes in data buffer */ - u16 status_qualifier; /* Status qualifier */ + __virtio32 sense_len; /* Sense data length */ + __virtio32 resid; /* Residual byte...

...priority field */ u8 crn; - u32 pi_bytesout; /* DataOUT PI Number of bytes */ - u32 pi_bytesin; /* DataIN PI Number of bytes */ + __virtio32 pi_bytesout; /* DataOUT PI Number of bytes */ + __virtio32 pi_bytesin; /* DataIN PI Number of bytes */ u8 cdb[VIRTIO_SCSI_CDB_SIZE]; } __packed; /* Response, followed by sense data and data-in */ struct virtio_scsi_cmd_resp { - u32 sense_len; /* Sense data length */ - u32 resid; /* Residual bytes in data buffer */ - u16 status_qualifier; /* Status qualifier */ + __virtio32 sense_len; /* Sense data length */ + __virtio32 resid; /* Residual byte...

...<- runif(n, 0, 1) f0 <- function(x) 2 * sin(pi * x) f1 <- function(x) exp(2 * x) f2 <- function(x) 0.2*x^11*(10*(1-x))^6+10*(10*x)^3*(1-x)^10 f <- f0(x0) + f1(x1) + f2(x2) g<-exp(f/4) y<-rpois(rep(1,n),g) mean.y <- mean(y) gam2 <- gam(y~ s(x2), poisson) # to plot on the response scale val.for.pred <- data.frame(x2=seq(min(x2), max(x2), length.out=100)) pred.2.resp <- predict.gam(gam2, val.for.pred ,type="response", se.fit=TRUE) lower.band <- pred.2.resp$fit - 2*pred.2.resp$se.fit upper.band <- pred.2.resp$fit + 2*pred.2.resp$se.fit pred.2.resp &...

Hi! I would like to be able to create formulas automatically. For example, I want to be able to create a function that takes on two values: resp and x, and then creates the proper formula to regress resp on x. My code: fit.main <- function(resp,x) { form <- expression(paste(resp," ~ ",paste(x,sep="",collapse=" + "),sep="")) z <- lm(eval(form)) z } main <- fit.main("y",c(&q...

dear R wizards:? I am wrestling with reshape.? I have a long data set that I want to convert into a wide data set, in which rows are firms and columns are years. > summary(rin) firm fyear sim1 Min. :1004.00 Min. :1964.0 Min. : -1.00000 1st Qu.:1010.00 1st Qu.:1979.0 1st Qu.: -0.14334 Median :1016.00 Median :1986.0 Median : 0.00116 Mean

...hole script is as follows: setwd("E:/") ##### (1) Read data ### dat <- read.table("Komfort_Tatsaechliche_ID_Versuchsreihe_1.txt", header=TRUE, sep="\t", dec=",") # tab-getrennte Textdatei names(dat) <- c("id", "resp", "amp") dat$id <- factor(dat$id) # kategoriale Variablen als Faktor definieren dat$amp <- as.numeric(gsub("^([0-9]),([0-9]+) Nm$", "\\1.\\2", dat$amp)) ##### (2) Explorative plots ##### library(lattice) xyplot(resp ~ amp|id, dat, type=c("g&quo...

Hi all, I have one array of predictors, one observation per row, and one array of responses, also arranged one observation per row. I arrange these into a data.frame and call lm() with a pasted-together formula. I would like to call lm() with a number of responses in excess of 100, but for some reason, 39 seems to be a limit. Why do I get an "invalid variable names" e...

Why when I call resp.getHeader(''Location''); after my ajax call does the browser get redirected to the location header. I just want to get the value of it. --~--~---------~--~----~------------~-------~--~----~ You received this message because you are subscribed to the Google Groups "Ruby on Rai...