similar to: Customized R Regression Output?

Dear all, I have some trouble with special characters while building my R package. I tried to follow the usual LATEX format, but could not fix the problem: For instance, for "greater than or equal", I tried "\geq", but R says that this is an unknown macro. Could anyone direct me how to solve this issue? Best Ozgur ----- ************************************ Ozgur ASAR

Hello, I have a dataframe (Lx) with 5 Lb, and 5 Lw variables. I want to create several variables according to the loop presented below (data attached). Lx <- read.csv("Lx.csv", header=T, sep=",") for (i in 1:20) { Lx$sb1[i] <- Lx$Lb1[i+1]/Lx$Lb1[i] Lx$sb2[i] <- Lx$Lb2[i+1]/Lx$Lb2[i] Lx$sb3[i] <- Lx$Lb3[i+1]/Lx$Lb3[i] Lx$sb4[i] <-

Dear all, Actually I could not decide who to contact, then decided to post here. If not appropriate sorry for that. I have written some functions in R which might be supposed to available in the base package but not available. I am wondering that how can I submit those functions to that package? Best Ozgur ----- ************************************ Ozgur ASAR Research Assistant Middle East

Hi, Your mistake seems to be in sum(v[1:x]) You create "x" as a vector but your treat it as a single number. v[1:x] expects "x" to be a single number and only considers its first element which is 1. If I understand your query correctly, the following might handle your problem: sum.vec <-NULL for (x in 1:1000){ t <- rbinom(1000, 1, 0.5) v <- replace(t,t==0,-1)

hello, I try to run simulation of levene's test to find the p-value but the error of replacement has length zero occur, could anyone help me to fix this problem? asim <- 1000 pv<-rep(NA,asim) for(i in 1:asim) {print(i) set.seed(i) g1 <- rnorm(20,0,2) g2 <- rnorm(20,0,2) g3 <- rnorm(20,0,2) x <- c(g1,g2,g3) group<-as.factor(c(rep(1,20),rep(2,20),rep(3,20))) library(Rcmdr)

Hi group, I've been thinking of calculating euclidean distance between each column of a data frames that each consists of standardized numerical columns. However, I don't know if there's a way of summarizing the overall distance by some kind of metrics. If anyone know a proper way of doing so and/or a package I would greatly appreciate your suggestions. Thanks very much! Kel --

Dear list users, I have a data frame as below specified. From the 1st of May to the 30th of September of several years (e.g. from 2004 to 2011) I have a frequency of accidents. I need the mean of accidents divided by weeks (i.e. the mean of accidents from the 1st to the 7th of May of all the years, from the 8th to the 14th of May, ..., from the 29th to the 31st of May, from the 1st to the 7th of

Dear all, I want to concatenate the elements of two vectors such as a<-c("a1","a2") b<-c("b1","b2") and obtain "a1b1", "a1b2","a2b1","a2b2" I tried the paste and paste0 functions, but they yielded elementwise concatenation such as "a1b1","a2b2" I am wondering that is there an efficient

Dear all, I want to manipulate a character string such as ex<-"cbind(data$response1,data$response2)" in R in two ways: 1) extracting the "response1" portion of ex 2) replacing "$" with "." I am wondering that is it possible efficiently doing these in R? Best Ozgur -- View this message in context:

Dear R users, I'd like to remove some descriptions when I use "filter". > filter(1:10, rep(1, 3)) Time Series: Start = 1 End = 10 Frequency = 1 [1] NA 6 9 12 15 18 21 24 27 NA That is, I want only this [1] NA 6 9 12 15 18 21 24 27 NA Thank you in advance. Kathie -- View this message in context:

Hi Do I have to load/download the whole package over and over again everytime when I use the R-Project? It often prompt me the error message "could not find function"- but I'm rather sure I have used the function before so I'm not entirely sure it couldn't find it again when I shut down the computer or so. Is it to do with directory or clashing with other software? Any

Dear Researchers, I am looking for to read a SAS macro in R. Although I searched it on web, I couldn’t find anything. Can you help me or direct me? Thank you for your interest and patience. Best. Ozgur [[alternative HTML version deleted]]

Hi Friends, I'm trying to model the consumer decisions (Click-Through Rate and Conversion) in Search Engine Advertising using a hierarchical Bayesian binary logit. The input data is the weekly CTRs and Avg. Position for each search keyword. CTR is modeled as (for each keyword i and week j): Pij = exp(C + Bi x Positionij + A1 x Lengthi + A2 x Brandi + A3 x ProductSpecifici) / [1 + exp(C +

Hi, I'm trying to get maximum likelihood estimates of \alpha, \beta_0 and \beta_1, this can be achieved by solving the following three equations: n / \alpha + \sum\limits_{i=1}^{n} ln(\psihat(i)) - \sum\limits_{i=1}^{n} ( ln(x_i + \psihat(i)) ) = 0 \alpha \sum\limits_{i=1}^{n} 1/(psihat(i)) - (\alpha+1) \sum\limits_{i=1}^{n} ( 1 / (x_i + \psihat(i)) ) = 0 \alpha \sum\limits_{i=1}^{n} (

HI >i have a little problem please help me to solve it >this is the code in R: >> beta0 [1] 64.90614 > beta1 [1] 17.7025 > beta [1] 17 64 >her beta<- c(beta0, beta1) thank you in advance hafida -- View this message in context: http://r.789695.n4.nabble.com/Decimal-number-tp4639428.html Sent from the R help mailing list archive at Nabble.com.

Dear Users, I am new to both of things, so do not blame me too much... I am busy with semiparametric regression and use WinBUGS to sample posteriors. The code to call Winbugs is as follows: data <- list("y","X","n","m") #My variables inits.beta <- rep(0,K) inits.beta0 <- 0 inits <-

Hi all-- We are planning an intervention study for adolescent alcohol use, and I am planning to use simulations based on a hurdle model (using the hurdle() function in package pscl) for sample size estimation. The simulation code and power code are below -- note that at the moment the "power" code is just returning the coefficients, as something isn't working quite right. The

Hi All, I have tried every fix on my try or tryCatch that I have found on the internet, but so far have not been able to get my R code to continue with the "for loop" after the lmer model results in an error. Here is two attemps of my code, the input is a 3D array file, but really any function would do.... metatrialstry<-function(mydata){ a<-matrix(data=NA, nrow=dim(mydata)[3],

I use "while" loop but it produces an errro. I have no idea about this. Error in "[<-"(`*tmp*`, i, value = numeric(0)) : nothing to replace with The problem description is The likelihood includes two parameters to be estimated: lambda (=beta0+beta1*x) and alpha. The algorithm for the estimation is as following: 1) with alpha=0, estimate lambda (estimate beta0

Hello, I've written a simple (although probably overly roundabout) function to test whether two regression slope coefficients from two linear models on independent data sets are significantly different. I'm a bit concerned, because when I test it on simulated data with different sample sizes and variances, the function seems to be extremely sensitive both of these. I am wondering if