similar to: Variance Inflation factor

Dear All, This question sounds very simple but I don't know where I am wrong. I just want to pad leading zeros in some string, for example, "123" becomes "00123". What is wrong if I do following? > sprintf("%05s", "123") [1] " 123" It didn't return "00123", instead it padded with 'blank'. Thank you for your help

Hi All, I have two Rs, one has been installed in Windows system and another one has been installed under UNIX system. Is there any environmental variable or function to tell me which R I am using? The reason that I need to know it is under different system, the data path could be different. I want to do something like if it is R under Windows path =

Hi all, I have a dataset. Each time I want to sample N(i) elements from it and I want to repeated sampling M times. N(i) is varied from time to time. For example, dataset = 1:50; a = list(); M = 1000; I want to do something like for(i in 1:M) { a[[i]] = sample(dataset,

Hi All, I have a data frame, whose colnames like "A.x, B.x, C.x, A.y, B.y, C.y". There could be many columns like this pattern. I want to compare data in columns A.x with A.y, B.x with B.y and C.x with C.y etc. Suppose my data frame is d, names(d) = c("A.x", "B.x", "C.x", "A.y", "B.y", "C.y"); If I want to D1 =

Hi All, I have a self-cooked package and save it to a zip file after running make, say named xxx.zip. After installing it to R by running "Install packages from local zip files" under Packages menu in R (Windows), I realized I needed to change some source codes and re-make it. My question is that if I need re-install that package or R would automatically update the R

Hi All, I have a data frame like a = data.frame(date = c(20081201, 20081202, 20081201), product = c("a b c d e", "a b c g h t", "d e h a c e h g"), sales = c(1, 2, 3)) Now I want to aggregate the sales by part of the a$product. 'Product' is the product name, a string separated by a space. The key in my aggregate function is

Hi All, I have a question about package installation in R. We have developed a package, say 'ABC'. We have installed it in two machines, A and B by running 'Install Package(s) from local zip file'. Everything was fine. Right now, suppose that package got damaged in machine A and our zipped file is gone, My question is that may I directly copy ../library/ABC from machine B to

I need to apply a yearly inflation factor to some wages and supply some simple sums by work category. I have gone at it with a brute force "for" loop approach which seems okay as it is a small dataset. It looks a bit inelegant and given all the warnings in the Intro to R, etc, about using loops I wondered if anyone could suggest something a bit simpler or more efficent? Example:

Hi All, Suppose I have the following codes: x = data.frame(A = rnorm(20), B = rnorm(20), C = rnorm(20)) a = list() a[["1.1"]] = x a[["1.2"]] = x b = list() b[["1.1"]] = c("A", "B") b[["1.2"]] = c("B", "C") Now I want to apply b to a like this, for each element of 'a', only select the corresponding columns

Hi All, I have two questions regarding install.packages(). Q1: may I run it in non-interactive mode, which means just install the packages I want rather than letting me to choose which mirror etc? Q2: I ran a R code in batch mode, for example, R CMD BATCH a.r In a.r, I have some statements like packages = library()$result[, "Package"]; f(! "plyr" %in% packages) {

Hi All, I know in R there is function named 'step', which does the stepwise regression and choose the model by AIC. However, if I want to choose a model per this logic: 1. Run a full model (linear regression, f = lm(y ~., data = ZZZ), for example) 2. Pick up the variable with biggest p value, delete it from the module and get a new regression model. 3. Repeat step 2

Hi All, Does somebody know how to know the detail of the line types? For example, lty = 1, means what kind of line?, lty = 2, means what kind of line? Thanks. HXD [[alternative HTML version deleted]]

Hi all, I am planning to parse some information on a website which includes lots of Chinese characters. Does someone know how to read/display Chinese in R? Thanks. url = "http://www.teec.org.cn/html/renwujieshao/" x = readLines(url) I tried encoding = 'UTF-8' already but it didn't help. My R version is $platform [1] "i386-pc-mingw32" $arch [1] "i386"

Hello, Seeking a way to get the variance inflation factor VIF for a model of regression in complex survey, I have understood that without this package (SURVEY) RGui VIF obtained as follows: fit <- lm(mpg~disp+hp+wt+drat, data=mtcars) vif(fit) But I want to know if survey, Vif is obtained so vif( svyglm(api00~ell+meals+mobility, design=dstrat)) Thank you, happy day

Hi All, The following is a snippet of my code. It works fine but it is very slow. Is it possible to speed it up by using different data structure or better solution? For 40000 runs, it takes 8 minutes now. Thanks a lot fun_activation = function(s_k, s_hat, x_k, s_hat_len) { common = intersect(s_k, s_hat); if(length(common) != 0) { index = match(common, s_k);

Hi group: I have a data set, which has severe colinearity problem. While running linear regression in R and SPSS, I got different models. I am wondering if somebody knows how to make the two software output the same results. (I guess the way R and SPSS handling singularity is different, which leads to different models.) Thanks. [[alternative HTML version deleted]]

hello all - i was searching for theoretical articles/vector equations regarding variance inflation factor (or generalization) for the linear mixed effects model (repeated measures data) sincerely, tom -- View this message in context: http://n4.nabble.com/variance-inflation-factor-for-linear-mixed-effects-tp1009143p1009143.html Sent from the R help mailing list archive at Nabble.com.

Hi everyone, Running the vif() function from the car package like ---------------------------------------------------- > reg2 <- lm(CARsPur~Delay_max10+LawChange+MarketTrend_20d+MultiTrade, data=data.frame(VarVecPur)) > vif(reg2) Delay_max10 LawChange MarketTrend_20d MultiTrade 1.010572 1.009874 1.004278 1.003351

Hi All, Suppose I have data like [[alternative HTML version deleted]]

I have the following script, how can I implement to achieve that calculate the VIF. Thanks. U1.7km<-c(15:24) R<-c(1.2,0.2,3.6,2.5,4.8,6.3,2.3,4.1,7.2,6.1) Hm<-c(1:10) mod<-nls(R~a*(U1.7km^b)*(Hm^c), start=list(a=2.031, b=0.800, c=-0.255), trace=T) summary(mod) coef(mod) coef(summary(mod)) -- View this message in context: