[R] Rcolorbrewer Package

Hi all, I am trying to download the Rcolorbrewer package from Cran http://cran.r-project.org/web/packages/RColorBrewer//index.html It seems the files have been removed. Does anyone know where can I download the package? Thanks. Hannah [[alternative HTML version deleted]] _

Re: [R] Overlay Histogram

e.y=theme_blank()) grid.arrange(xhist, empty, none, yhist, ncol=2, nrow=2, widths=c(3, 1), heights=c(1, 3)) 2012/8/7 li li > Please ignore my previous message. The limits are correct. > Thanks again. > Hannah > > 2012/8/7 li li > >> Hi John, >>

[R] Density

Dear all, Given a set of observations X, I want to evaluate the kernel density estimator at these observed values. If I do the following, I do not get the those estimated values directly. Can anyone familiar with this give an idea on how to find out the estimated density values for X? > X <- rn

Re: [R] Density

amp;M University > College Station, TX 77843-4352 > > > > -Original Message- > > From: r-help-boun...@r-project.org [mailto:r-help-bounces@r- > > project.org] On Behalf Of li li > > Sent: Wednesday, August 08, 2012 9:03 PM > > To: r-help > > Subject: [R]

[R] Import Data from Excel

Dear all, I want to import just part of an excel data file into R. I would like to have the data imported without rownames or colume names. I used read.delim("clipboard", header=F). Somehow even though I added the argument "header=F", I still have the row names V1, V2, ..., Does anyone kno

[R] nls()

Dear all, I have a question on the four parameter logistic regression. Not sure about how to set the start values for the parameters. Can anyone familiar with this give some hint? Thanks so much. Hannah For example, suppose I have the following data values, > temp mean_st

[R] Apply a function according to factor levels.

Dear all, I have a question on applying a function to the data according to factor levels. For example, for the data below, what is the best way to apply a function to "values" according to different levels of "samples" (1,2,3,4,5)? values ind sample 1 0.03325 1 1 2 0.03305 1

Re: [R] Apply a function according to factor levels.

t; #1 1 0.033556 > #2 2 0.011120 > #3 3 0.017750 > #4 4 0.030220 > #5 5 0.062060 > A.K. > > - Original Message - > From: li li > To: r-help > Cc: > Sent: Thursday, August 16, 2012 9:31 PM > Subject: [R] Apply a function acco

[R] Linear mixed model using R

Dear all, I am trying to use R to fit mixed models. Take the following example, where ind is a random effect and sample is fixed. I wanted to fit Model 1: values = ind + sample Model 2: values =ind * sample Model 3: values=ind(sample) + sample Tried to use the below for mod1, but it

Re: [R] Linear mixed model using R

Thanks, Bert. 2012/8/17 Bert Gunter > 1. Post on R-sig-mixed-models, not here. > > 2. Models 2 and 3 make no sense (to me, anyway). What do you think > they mean? (Don't answer here -- explain on the mixed models list). > > -- Bert > > On Thu, Aug 16, 2012 at 9:13 P

[R] mixed model in r

Dear all, I am really confused with the syntax for mixed model using lmer in lme4 library. My model is Value = Sample + Run (nested with Sample ) Here sample is a fixed effect and Run is a random effect nested with sample. Can anyone give some hint on the correct syntex for this model.

[R] Latex question

Dear all, This may not be the right place for latex questions, but I do not know where else to ask the question. I have two subfigures. I want the first subfigure to have a caption on the side of it and the second subfigure to have a caption below it. Thank you. Hannah [[

[R] qmvnorm function

Dear all, I need to use the "qmvnorm" function in mtvnorm package. Here is the error message I got > qmvnorm(0.05, tail="upper", sigma=var)$quantile Error in uniroot(pfct, interval = interval) : f() values at end points not of opposite sign There is no problem for 50th quantile. > qmvnorm(0.

Re: [R] qmvnorm function

gma=var)$quantile [1] 0.8935925 > ÔÚ 2012Äê3ÔÂ14ÈÕ ÏÂÎç6:37£¬Rui Barradas дµÀ£º > Hello, > > li li-13 wrote > > > > Dear all, > >I need to use the "qmvnorm" function in mtvnorm package. > > > > Here is the error message I got > >> qmvnor

[R] R help

Dear all, I need to generate numbers from multivariate normal with large dimensions (5,000,000). Below is my code and the error I got from R. Sigma in the code is the covariance matrix. Can anyone give some idea on how to take care of this error. Thank you. Hannah > m <- 500

Re: [R] R help

. Hannah ÔÚ 2012Äê2ÔÂ19ÈÕ ÉÏÎç5:01£¬Petr Savicky дµÀ£º > On Sat, Feb 18, 2012 at 06:00:53PM -0500, li li wrote: > > Dear all, > > I need to generate numbers from multivariate normal with large > dimensions > > (5,000,000). > > Below is my code and the error I got fr

Re: [R] R help

error code=12) *** error: can't allocate region *** set a breakpoint in malloc_error_break to debug R(173,0xa0a7b540) malloc: *** mmap(size=43072) failed (error code=12) *** error: can't allocate region *** set a breakpoint in malloc_error_break to debug ÔÚ 2012Äê2ÔÂ19ÈÕ ÉÏÎç10:5

[R] Some Help Needed

Dear all, I need to do some calculation where the code used are below. I get error message when I choose k to be large, say greater than 25. The error message is "Error in integrate(temp, lower = 0, upper = 1, k, x, rho, m) : the integral is probably divergent". Can anyone give some help on res

[R] means by group

Hi all, I have a 100 by 100 matrix and I divided this matrix into 100 groups, each is a 10 by 10 submatrix. I want find out the means of each group. I know we can use apply function for mean by margins. Is there a function in R for means by groups also? Thanks. Hannah [[alter

[R] Some help needed

Dear all, For the following code, I have the error message "Error in uniroot(f1star, lower = -10, upper = 0, tol = 1e-10, lambda = lam[i], : f() values at end points not of opposite sign". It seems the problem occurs when lambda is equal to 0.99. However, there should be a solution for "f1

[R] heatmap

Dear all, I have a data set which is 100 \times 100 matrix. I made a heatmap of the data. The xlim and ylim are both shown to be c(0,1). I want to make both the xlim and ylim to be c(0,100). So the heatmap correspond to the 100 \times 100 locations. Thank you. Hannah [[alternative

[R] Integration of two dimension function

Dear all, I have a function f <- function(x,y,c){as.numeric(x*y < c) } I need to solve the value of c so that when I take the integral of the function f from 0.2 to 0.8 with respect to x and from 0 to 1 with respect to y, the integral equal some prefixed value, say 0.025. It involves two dim

[R] Questions on RWeka classifiers?

Hi, I am using some classifiers in RWeka packages and met a couple problems. (1) J48 implements C45 classifier, the C45 should be able to handle missing values in both training set and test set. But I found the J48 classifier can not be evaluated on test set with missing values--it just

[R] Questions on RWeka classifiers

Hi, I am using some classifiers in RWeka packages and met a couple problems. (1) J48 implements C45 classifier, the C45 should be able to handle missing values in both training set and test set. But I found the J48 classifier can not be evaluated on test set with missing values--it just

[R] how to check the model adequacy for a cox proportional hazard model fit?

Anybody has any suggestion? Thanks! [[alternative HTML version deleted]] __ R-help@r-project.org mailing list https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read the posting guide http://www.R-project.org/posting-guide.html and p

[R] How can C++ read the R object written into socket with saveRDS or save

Hi, all, Recently, I met one issue when using socket between R & C++ to transmit R object. Would you pls help give me some suggestions? Many thanks! [Background]: I create a socket connection between R & C++ binary first, and then, want to use saveRDS() or save() in R to save the object into con

<    1   2   3