search for: nausea

> I am using a random intercept model with SITEID as random and NAUSEA as > outcome. > > I tried using a dataset without missing values and changed my model > statement accordingly but still get the same error. Follwoing in an excerpt. > > anal.data <- read.table("nausea.txt", header=T, sep="\t") > > list(names(anal.data)...

...I have a list of patient identifiers and diagnoses from inpatient admissions. I would like to reorganize the list, presently in a long format to a wide format in reshape, but in the absence of a "time" element, I am uncertain how to do this - any help greatly appreciated. ID Dx A nausea A diabetes A kidney failure A heart attack A fever B fever B pneumonia B heart attack B nausea B cough C kidney failure C nausea C foot pain I want this to be ID nausea diabetes kidney_failure heart_attack fever pneumonia A 1 1...

R-experts: I would like to explore the partial proportional odds models of Peterson and Harrell (Applied Statistics 1990, 39(2): 205-217) for a dataset that I am analyzing. I have not been able to locate a R package that implements these models. Is anyone aware of existing R functions, packages, etc... that might be used to implement the partial proportional odds models? Brant Inman

I am using a random intercept model with SITEID as random and NAUSEA as outcome. Thanks. Anamika I tried using a dataset without missing values and changed my model statement accordingly but still get the same error. Follwoing in an excerpt. > anal.data <- read.table("nausea.txt", header=T, sep="\t") > list(names(anal.data)) [[1]] [1]...

Hi:   I am using the following code to read a data set in txt in R and using the follwoing model. But it seems to give me an error ' expected variable name error pos 134022'. Any help is greatly appreciated.   Code: anal.data <- read.table("nausea.txt", header=T, sep="\t") list(names(anal.data)) attach(anal.data) n.samples <- dim(anal.data) [1] # number of data points n.params <- dim(anal.data)[2] - 2  # number of parameters to estimate n.sites <- length(unique(SITEID))                                                 ...

...that comes out shows my parameters as nos with 1 decimal place after. I would want to have the parameters with 5 places after decimal. How would I specify that in my code for R2WinBUGS. trial.data <- read.table("trial_data.txt", header=T) bugs.output <- list() for(i in 1:5){ nausea <- as.integer(trial.data[i,]) bugs.output[[i]] <- bugs( data=list(nausea=nausea, N=63), inits=list( list(alpha0=0,tau=1), list(alpha0=0.9, tau=1), list(alpha0=0.8,tau=1) ), model.file="B-LN_model_tria...

...error message (please see output below) no matter what shape 1 or 2 I use. It works perfect with the cauchy or normal priors. Do you know if there is a catch there somewhere? Thanks logpriorfun <- function(beta,shape1,shape2){ sum(dbeta(beta,shape1,shape2, log=T)) } posterior <- MCMClogit(nausea~nsaid*diuretic, data=w, verbose=2000,burnin = 1000, mcmc = 10000, user.prior.density=logpriorfun,shape1=1,shape2=1) user.prior.density(beta.start) == -Inf. Error in MCMClogit(nausea ~ nsaid * diuretic, data = w, verbose = 2000, : Respecify and call MCMClogit() again. -- View this mess...

...NSAID*Diuretic) 0.85407 0.30603 2.791 0.00526 ** --- Signif. Codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1 Odds ratio is 2.35 [ln(0.85407)] times higher when NSAID is present in addition to Diuretic. ------------------------------- Odds ratio of Nausea when on Diuretic is exp(0.25847)= 1.29 and the odds ratio of Nausea when on NSAID is exp(0.13630)=1.14 Normally when we want to see the odds ratio of Nausea when a patient is on both drugs we multiply 1.29*1.14= 1.48 (is this correct? do we multiply or do we add?) But since the interaction term is...

...e I am using > > library("R2WinBUGS") > > > > #Set working directory > > setwd("E://Achaudhuri/Trial") > > trial.data <- read.table("trial_data.txt", header=T) > > bugs.output <- list() > > for(i in 1:5){ > > nausea <- as.integer(trial.data[i,]) > > bugs.output[[i]] <- bugs( > > data=list(nausea=nausea, N=63), > > inits=list( > > list(alpha0=0,tau=1), > > list(alpha0=0.9, tau=1), > > list(alpha0=0.8,tau=...

...the following code to generate Beta-Binomial data. Now I want to add a covariate to the equation. I would then like to use the simulated data to run the Beta-Binomial model with covariates on it. Appreciate any help. set.seed(111) k<-20 n<-60 x<-NULL p<-rbeta(k,3,3)# so that the mean nausea rate is alpha/(alpha+beta) for(i in 1:k) x<-cbind(x,rbinom(300,n,p[i])) Thanks Anamika [[alternative HTML version deleted]]

Hi: I have asked this question on Cross-Validated. So it might be a cross posting but havent received any responses to it. I am trying to see which distribution will best fit the data I am working on. The dataset is as following: Site Nausea headache Abdominal Distension 1 17 5 10 2 12 8 7 ..... So each site has total # adverse events for each type/category and have equal # patients per site, say 60 and...

HI All: Does anyone know the code behind the qbeta function in R? I am using it to calculate exact confidence intervals and I am getting 'NaN' at places I shouldnt be. Heres the simple code I am using: k<-3 > x<-NULL > p<-rbeta(k,3,3)# so that the mean nausea rate is alpha/(alpha+beta) > min<-10 > max<-60 > n<-as.integer(runif(3,min,max)) > for(i in 1:k) + x<-cbind(x,rbinom(5,n[i],p[i])) > > # Exact Confidence Interval > > l_cl_exact<-qbeta(.025, x, n-x+1) Warning message: In qbeta(p, shape1, shape2, lower.tail, lo...

...JRClient to build logistic regression model in the following manner : Rconnection c = new Rconnection(); c.eval("KSN<-read.table(\"/Users/amine/Documents/Research/ Tools/R/D2R1.txt\",header=T,sep=\",\")"); c.eval("result <- glm(Nausea ~ Kaletra*Sustiva, family = binomial(link = logit), data =KSN)"); I can obtain the coef. using: REXP coeffExp = c.eval("coefficients(result)"); double[] coeff = coeffExp.asDoubleArray(); how could I transfer the p-values, i.e., Pr (>|z|), for each coef., in java vecto...

...because I need a solution that retain certain ways of doing things because they give me the information i need. I want to examine the effect of increasing my sampling interval on my data. Example: what if instead of sampling every hour I sampled every two, oh yeah, how about every three?.. etc ad nausea. How I want to do this is to take the data I have now, add an index to it, that contains counters. Those counters will look something like 1,2,1,2,.. for the first one, 1,2,3,1,2,3.. for the next one. I have a lot of them, like say a thousand... Then for each column in the index my loops should s...

Hi, Any problems with erejction? try VIjAGRA - http://www.porutagundes.com only from the soles of my boots. The blood rushed to my head-as well it might bringing with it a surge of nausea and a sensation of great unease.

Hello ... >From the icecast website (www.icecast.com) i found the following info > (i copy-pasted it): { Note that icecast-1 is deprecated and unsupported, and may possibly induce nausea in small children, hair loss in men, and estrogen fluctuations in women. If, for some reason, you absolutely must have ancient versions of icecast (possibly for historic preservation), download icecast 1.3.12. Note that this is deprecated and unsupported. } I have a project whereby it is require...

Yes, I know this is lame, but my location limits me to using a dial-up ISP. I am running asterisk with a T100P and a TDM400 card. I currently have dial-on-demand setup on the same box, using diald and an external modem. To prevent DOD from trying to dial out during an external call, I have an AGI script kill diald before initiating a call, and then restart it at the end of the call by looking

Hi, My application has a table of People that has information about them, and I want to have a subset of those people as Users - people that can login to the application. What is the best way to do this? Add columns to People that only Users would use? I''ve read the threads on single table inheritance, but it seemed a little bit wasteful: 10,000 people vs 50 users For those more

Hi to all, Here is a question which I cannot solve. Appreciate so much for any suggestions! I have a squared region which is irregularly divided into many rectangular patches. Each patch is associated with a value, and two patches possibly share a common value. I hope to fill in each patch a pattern according to its value. For instance, if a patch has value 1, then I fill in that patch

Hi Assorgia, Thanks for the flash animation. I think Icecast1 might be the solution. I am just a little bit scared by the message "Note that icecast 1 is deprecated and unsupported, and may possibly induce nausea in small children, hair loss in men, and estrogen fluctuations in women" that appears on the download page. I have a nice blond hair of "Principe Azzuro" that I don't want to lose! Just kidding. Anyway, I am trying a little more with Icecast2 but if it turns impossible (what...