[R] bootstrap based confidence band

I'm trying to find a bootstrap based confidence band for a linear model. I have created a data set with X and Y X=runif(n,-1.25,1.25) e=rnorm(n,0,1) Y=exp(3*X)+5*sin((30*X)/(2*pi))+2*e fit=lm(Y~X) summary(fit)   I define a bootstrap function named PairedBootstrap which is not listed here. Than I t

[R] repeat matrix rows as a whole

Hi, If I have a matrix like 1  2  3  4 5  6  7  8, how can I repeat two rows as whole, to be like 1  2  3  4 5  6  7  8 1  2  3  4 5  6  7  8?   Since I have more two rows in a matrix and I need to repeat many times, I wonder whether there is a convenient command to do so.   Thanks! [[alte

[R] reorder a matrix

Hi, Here are part of my data, > pr16.5    origin dest ldco time.slot distortion 1   1    1    1 1    4.3 2   1    1    1 2    4.7 3   1    1    1     3    5.6 4   1    1    1 4    7.7 5   1    1    2 1    6.8 6 

[R] reorder a matrix

Hi, Here are part of my data, > pr16.5    origin dest ldco time.slot distortion 1   1    1    1 1    4.3 2   1    1    1 2    4.7 3   1    1    1     3    5.6 4   1    1    1 4    7.7 5   1    1    2 1    6.8 6  

[R] find the eigenvector corresponding to the largest eigenvalue

Hi, If I use the eigen() function to find the eigenvalues of a matrix, how can I find the eigenvector corresponding to the largest eigen value? Thanks! [[alternative HTML version deleted]] __ R-help@r-project.org mailing list https://stat.ethz.c

[R] Efficiently parallelize across columns of a data.table

I am trying to parallelize a task across columns of a data.table using foreach and doParallel. My data is large relative to my system memory (about 40%) so I'm avoiding making any copies of the full table. The function I am parallelizing is pretty simple, taking as input a fixed set of columns and

[R] Efficiently parallelize across columns of a data.table

Makes sense. Thanks for the clear explanation. Rebecca On Friday, August 19, 2016, Peter Langfelder > wrote: > Last time I looked (admittedly a few years back), on unix-alikes > (which you seem to be using, based on your use of top), > foreach/doParallel used forking. This means

[R] Error: OutOfMemoryError (Java): GC overhead limit exceeded

a): GC overhead limit exceeded" However, this seems not working any more, is there any other way to help avoiding the memory error issue in R? Thanks very much! Cheers, Rebecca 20 -- This message, and any attachments, i

Re: [R] Error: OutOfMemoryError (Java): GC overhead limit exceeded

Hello Jim, I already use “.xls” for the loading, but still have the memory issue…. Thanks, Rebecca From: jim holtman [mailto:jholt...@gmail.com] Sent: Monday, February 02, 2015 10:21 AM To: Jeff Newmiller Cc: Yuan, Rebecca; R help (r-help@r-project.org) Subject: Re: [R] Error

[R] randomForrest-imputation

Hi! How can I generate only positive data with randomForrest-imputation? I'm working with laboratory values which are always positive. Can anybody help out? Thanks! (P.S.: I needed to send this request again because there were problems by delivering it) __

[R] missRanger package

Hi, does anybody know where I need to insert the censoring in the missRanger package? Regards, Rebecca __ R-help@r-project.org mailing list -- To UNSUBSCRIBE and more, see https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read the posting guide

[R] Error : '.path.package' is defunct.

;) Attaching package: 'zoo' The following objects are masked from 'package:base': as.Date, as.Date.numeric Could you please help on this? Thanks! Rebecca -- This me

[R] Linking in R package documentation

name in another package? I thought I could use something like: \seealso{\code{\link[nlme]{ranef}}} but that gives a Page Not Found. Thank you in advance! Rebecca __ R-help@r-project.org mailing list https://stat.ethz.ch/mailman/listinfo/r-help PLEASE

[R] Correct nested design for GLM

cro is nest within Trt and Trt is nested with Year (0 + Micro/Trt/Year). I was hoping someone could confirm that this is actually what the model is running? Thank you very much for your time, Rebecca [[alternative HTML version deleted]] __ R-

[R] Zero inflated negative binomial

el 1 there is no Micro level 3. I cannot find a solution to this problem, is there any way I can solve this so I can run the zero inflated model? Thank you for your time Rebecca [[alternative HTML version deleted]] __ R-help@r-project.

[R] Problems with subsets in NLME

this list as well.) Thank you! Rebecca -- Rebecca Sela IOMS/Statistics Group Stern School of Business New York University __ R-help@r-project.org mailing list https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read the posting guide http://www.R-pro

[R] \dQuote in packages

UNKNOWN '\sideeffects' Warning in parse_Rd("./man/simpleREEMdata.Rd", encoding = "unknown") : ./man/simpleREEMdata.Rd:10: unknown macro '\item' Are \dquote, \sideeffects, and \item not supported in newer versions of R? Is there some und

Re: [R] \dQuote in packages

w to assemble a package in Windows.) Thank you for your help! Rebecca - Original Message - From: "Uwe Ligges" To: "Rebecca Sela" Cc: "r-help" Sent: Wednesday, July 8, 2009 6:11:33 PM GMT -05:00 US/Canada Eastern Subject: Re: [R] \dQuote in packages The d

Re: [R] \dQuote in packages

ted to the target variable)} } \references{Sela, Rebecca J., and Simonoff, Jeffrey S., \dQuote{RE-EM Trees: A New Data Mining Approach for Longitudinal Data}.} \keyword{datasets} Thanks again for your help! Rebecca - Original Message - From: "Uwe Ligges" To: "Rebecca

Re: [R] \dQuote in packages

Thank you! (That was easy to fix.) How does one deal with quoting (in \reference)? The following line causes problems: \references{Sela, Rebecca J., and Simonoff, Jeffrey S., \dQuote{RE-EM Trees: A New Data Mining Approach for Longitudinal Data}.} The error given is: Warning in parse_Rd

[R] RPART - printing full splitting rule number on tree plot

? Below is an example. library(rpart) temp <- as.data.frame(cbind(kyphosis, Start2=kyphosis$Start)) temp$Start2 <- temp$Start2+1 fit1 <- rpart(Age~ Kyphosis +Number + Start2, data=temp) plot(fit1) text(fit1, use.n=TRUE, digits=5) Thanks Dr Rebecca O&#

[R] Checking a (new) package - examples require other package functions

I refer to other functions from the package in examples? (I have the "Writing R-extensions" PDF, so it would be enough to point me to the right page, if the answer is in there and I just missed it.) Thanks! Rebecca __ R-help@r-project.org mail

Re: [R] Checking a (new) package - examples require other package functions

FALSE, domain = NA) 1: library(package, lib.loc = lib.loc, character.only = TRUE, verbose = FALSE) Execution halted I have attached the NAMESPACE file in case that is where the problem is. Thanks again for your help! Rebecca - Original Message - From: "Duncan Murdoch" To: &quo

[R] "Error: package/namespace load failed"

ompiled C or Fortran code, but I have none of that in my code. I imagine this is a simple problem (perhaps with my NAMESPACE file), but I don't know what it is. (The text of the NAMESPACE file is at the bottom of this e-mail.) Thank you in advance for your help! Rebecca NAMESPACE

[R] No CHM file

ckage was built under an older version of R? Or is there some other way to ensure that the Help files show up properly? Thanks in advance. Rebecca __ R-help@r-project.org mailing list https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read the posting

[R] generate individual panel keys/legends in lattice barchart

ter trying several iterations, I could not get it to work. any help you could give would be great, Regards, Rebecca. # code to reproduce my graph # my_subcat_data <- structure(list(subcat = structure(c(9L, 6L, 1L, 3L, 11L, 5L, 12L, 10L, 2L, 4L, 7L, 8L, 9L, 6L

Re: [R] Read many cvs files

I know that R can read cvs file, but could R read excel files? For example, in the workbook, I have 3 sheets. How could I read for those three sheets? Thanks, Rebecca -Original Message- From: r-help-boun...@r-project.org [mailto:r-help-boun...@r-project.org] On Behalf Of jim

[R] suggestions about import SAS results to R.

! Cheers, Rebecca -- This message, and any attachments, is for the intended r...{{dropped:5}} __ R-help@r-project.org mailing list https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do

Re: [R] suggestions about import SAS results to R.

- require(foreign) rm(list=ls()) sashome <- "C:/Program Files/SAS/SASFoundation/9.2/" checkfile <- read.ssd("C:/Documents and Settings/test","sales",sascmd = file.path(sashome, "sas.exe")) ----

Re: [R] suggestions about import SAS results to R.

Hello Ista, Thanks for mention this package, however, it will not support big file, where I got the error message: -- Error in read.sas7bdat("C:/Documents and Settings/test.sas7bdat") : big endian files are not supported -- Cheers

Re: [R] suggestions about import SAS results to R.

//Documents and Settings/test_20130103.xls",readOnly=TRUE); age20 <- sqlFetch(channel,"age=20",as.is=TRUE); odbcClose(channel) The error message shows that "External table is not in the expected format." -- Is this related to the warning message in

[R] How to cut the time format short

Hello all, I have a time format looks like 31JAN2002:00:00:00.000 How could I cut it to 31JAN2002 ? I tried with format() but not work. Thanks, Rebecca -- This message, and any attachments, is for the intended r

Re: [R] How to cut the time format short

Hello Ellison, Thanks very much! I use the substr() and it works! Have a great day! Cheers, Rebecca -Original Message- From: S Ellison [mailto:s.elli...@lgcgroup.com] Sent: Wednesday, January 16, 2013 1:35 PM To: Yuan, Rebecca; R help Subject: RE: How to cut the time format short

[R] A smart way to use "$" in data frame

5 23 datab$oldbalance[datab$oldstate== 'AR' ] 1234 5673 Could I have a way to quote data$balance[data$state == 'AR'] in general, where balance=oldbalance, state=oldstate when data=dataa, and balance = newbalance, state = newstate when data=datab ? Thank

Re: [R] A smart way to use "$" in data frame

Hello Greg, Thanks very much! This helps! Cheers, Rebecca From: Greg Snow [mailto:538...@gmail.com] Sent: Friday, January 18, 2013 5:17 PM To: Yuan, Rebecca Cc: R help Subject: Re: [R] A smart way to use "$" in data frame The important thing to understand is that $ is a shortcut

[R] plot two time series with different length and different starting point in one figure.

would have one data point from A and from Mar, 2012-Nov, 2012, it would have two data points from A and B, and in December 2012, it would have one data point from A. Thanks very much! Cheers, Rebecca -- This message, and any

Re: [R] plot two time series with different length and different starting point in one figure.

must specify uniquely valid column(s) To plot A and B in one plot is to compare the difference between them... Any other thoughts? Thanks, Rebecca -Original Message- From: PIKAL Petr [mailto:petr.pi...@precheza.cz] Sent: Tuesday

Re: [R] plot two time series with different length and different starting point in one figure.

NA NA Does R have a function to convert the date to some format of Sep,2012, therefore when I merge those two, they will not have those NA numbers... Thanks, Rebecca -Original Message- From: arun [mailto:smartpink...@yahoo.com] Sent: Tuesday, January 22, 2013 2:15 PM

Re: [R] plot two time series with different length and different starting point in one figure.

o not know how to do it. Thanks, Rebecca -Original Message- From: David Winsemius [mailto:dwinsem...@comcast.net] Sent: Tuesday, January 22, 2013 2:34 PM To: Yuan, Rebecca Cc: R help Subject: Re: [R] plot two time series with different length and different starting point in one figur

Re: [R] plot two time series with different length and different starting point in one figure.

Hello Arun, Thanks very much! In this way, it works! I convert both A and B to the same day of the month, and therefore there is no NA shown for different last business day and last calendar day of the month. You are very help! Cheers, Rebecca -Original Message- From: arun

[R] to check if a character string is in a group of character strings

Hello, How can I judge if a string is in a group of string? For example, I would like to have if (subpool in pool){ }else{ } Where > pool = c("s1","s2") > subpool = c("s1") How can I write the "subpool in pool" ri

Re: [R] to check if a character string is in a group of character strings

Hello Duncan, Thanks for this! This works! Best, Rebecca -Original Message- From: Duncan Murdoch [mailto:murdoch.dun...@gmail.com] Sent: Wednesday, January 23, 2013 11:23 AM To: Yuan, Rebecca Cc: R help Subject: Re: [R] to check if a character string is in a group of character

[R] fSeries not found in R

Hello all, When I tried to install fSeries in R, I got the following error messages: install.packages("fSeries",dependencies=T) Warning message: package 'fSeries' is not available (for R version 2.15.2) Is this package changing/merging to another packag

[R] Range difference of plot two arrays in one plot

oints from those two arrays are overlapping together, but the value 30 is not equal to 90. How could I draw those two arrays in one plot with their real values shown in the plot? i.e., x is in the lower part of the plot and y is in the upper part of the plot. Thanks very muc

Re: [R] Range difference of plot two arrays in one plot

Hello David, When I used the same ylim for both plots, I got what I need: x = c(0,0,0,10,20,30) y = c(40,50,60,70,80,90) plot(x, type='o', ylim=c(min(x,y),max(x,y))) par(new=T) plot(y, type='l', ylim=c(min(x,y),max(x,y))) Thanks, Rebecca -Original Message- F

[R] To convert a quarterly data to a monthly data.

to convert monthly data to quarterly data, not the way I need. Is there any library or function in R can do that? Thanks very much! Cheers, Rebecca -- This message, and any attachments, is for the intended r...{

Re: [R] To convert a quarterly data to a monthly data.

int(qtrly$rDate, qtrly[,names(qtrly)[i]] ,subset(mthly, rDate <= max(qtrly$rDate))$rDate )) where qtrly is the quarterly data and mthly is the monthly data. Cheers, Rebecca -Original Message- From: Rolf Turner [mailto:rolf.tur...@xtra.co.nz] Sent: Wednesday, February 20, 2013 3:14

[R] How to get the t-stat for arima()?

Hello all, fit = arima() and Summary(fit) will give some summary of the fit. However, the t-stats are not shown in the summary. How can I get the t-stats of it? Thanks, Rebecca -- This message, and any attachments, is for the

[R] How to look at the source code for predict()

Hello, I try to look at the source code of predict() it turns out that I cannot find it. I can see it with debug(library), but not efficient. Can someone help? Thanks, Rebecca -- This message, and any attachments, is for

Re: [R] How to get the t-stat for arima()?

Hello Rui, Thanks very much! This works! Cheers, Rebecca -Original Message- From: Rui Barradas [mailto:ruipbarra...@sapo.pt] Sent: Tuesday, March 19, 2013 6:26 PM To: Yuan, Rebecca Cc: R help Subject: Re: [R] How to get the t-stat for arima()? Hello, Sorry for the error, the sqrt(n

Re: [R] How to look at the source code for predict()

Hello all, I thought I found it, it is in the arima.R if I use arima to fit the model. But how could we see the source code of some function in R? Thanks, Rebecca From: Yuan, Rebecca Sent: Wednesday, March 20, 2013 2:38 PM To: R help Subject: How to look at the source code for predict

Re: [R] How to look at the source code for predict()

Hello MW, Thanks very much! I use > getS3method("predict","arima0") And be able to see the source code. Have a great day! Cheers, Rebecca -Original Message- From: R. Michael Weylandt [mailto:michael.weyla...@gmail.com] Sent: Wednesday, March 20, 2013 2:52

Re: [R] How to look at the source code for predict()

Hello Duncan, Thanks for this link! A very good article that has useful information! Best, Rebecca -Original Message- From: Duncan Murdoch [mailto:murdoch.dun...@gmail.com] Sent: Wednesday, March 20, 2013 2:53 PM To: Yuan, Rebecca Cc: 'R help' Subject: Re: [R] How to l

Re: [R] How to look at the source code for predict()

Hello Rui, This also works! Thanks so much! Cheers, Rebecca -Original Message- From: Rui Barradas [mailto:ruipbarra...@sapo.pt] Sent: Wednesday, March 20, 2013 3:34 PM To: Yuan, Rebecca Cc: 'R help' Subject: Re: [R] How to look at the source code for predict() Hello, In t

[R] How to store data frames into pdf file and csv file.

frame, such as > st = data.frame(est.aic, est.aic) > st est.aic est.aic.1 1 -1654.986 -1654.986 How could I save it to a .csv file? Thanks very much! Cheers, Rebecca -- This message

[R] Could I get the following stats from arima()?

How can I get the following stat from arima? Pr(>|t|) r2 adjust_r2 rmse Thanks, Rebecca -- This message, and any attachments, is for the intended r...{{dropped:5}} __ R-hel

Re: [R] How to store data frames into pdf file and csv file.

l.summary.out,collapse='\n'),family='mono') dev.off() I believe this time we try to use pdf to catch the graph via text, are we? Thanks, Rebecca From: jim holtman [mailto:jholt...@gmail.com] Sent: Thursday, March 21, 2013 9:44 AM To: Yuan, Rebecca Cc: R help S

[R] How could I specify "no interception term" in arima?

Hello, I cannot find it in ?arima, does someone know how to specify "no interception" term in arima? Thanks, Rebecca -- This message, and any attachments, is for the intended r...{

Re: [R] How to store data frames into pdf file and csv file.

Hello Jim, This graphics.off() helps me to release the previous pdf file that was had the error message. Thanks again for your help! Cheers, Rebecca From: jim holtman [mailto:jholt...@gmail.com] Sent: Thursday, March 21, 2013 10:07 AM To: Yuan, Rebecca Cc: R help Subject: Re: [R] How to store

Re: [R] How to store data frames into pdf file and csv file.

Hello Jim, Thanks very much! This time I am able to get the pdf file. Best regards, Rebecca From: jim holtman [mailto:jholt...@gmail.com] Sent: Thursday, March 21, 2013 10:04 AM To: Yuan, Rebecca Cc: R help Subject: Re: [R] How to store data frames into pdf file and csv file. Here is the

[R] easy way of paste

p)[1:5],sep="+") > a [1] "aa" "dummy1" "dummy2" "bb" "cc" But it will not give me the way I want. Thanks, Rebecca ---

Re: [R] easy way of paste

Hello Sarah, Yes, that is what I want. Thanks! Cheers, Rebecca -Original Message- From: Sarah Goslee [mailto:sarah.gos...@gmail.com] Sent: Thursday, March 21, 2013 11:21 AM To: Yuan, Rebecca Cc: R help Subject: Re: [R] easy way of paste It looks like you want collapse rather than

Re: [R] easy way of paste

Hello A.K., Thanks! I use collapse instead of sep and get the answer. Cheers, Rebecca -Original Message- From: arun [mailto:smartpink...@yahoo.com] Sent: Thursday, March 21, 2013 11:21 AM To: Yuan, Rebecca Cc: R help Subject: Re: [R] easy way of paste  dat2<-as.data.frame(matrix(1

Re: [R] How could I specify "no interception term" in arima?

Hello Rui, Thanks! I get it now. Should have read the documentation carefully! Cheers, Rebecca -Original Message- From: Rui Barradas [mailto:ruipbarra...@sapo.pt] Sent: Thursday, March 21, 2013 2:46 PM To: Yuan, Rebecca Cc: R help Subject: Re: [R] How could I specify "no interce

[R] error while extracting the p-value from adf.test

because the "-" between "p" and "value"? Thanks, Rebecca -- This message, and any attachments, is for the intended r...{{dropped:5}} __ R-help@r-projec

Re: [R] error while extracting the p-value from adf.test

Hello Sarah, Thanks! I will look into help next time before sending out the question. Thanks, Rebecca -Original Message- From: Sarah Goslee [mailto:sarah.gos...@gmail.com] Sent: Friday, March 22, 2013 10:22 AM To: Yuan, Rebecca Cc: R help Subject: Re: [R] error while extracting the p

Re: [R] error while extracting the p-value from adf.test

Hello José, This is a much convenient way of looking at the components! Thanks very much! Cheers, Rebecca -Original Message- From: Jose Iparraguirre [mailto:jose.iparragui...@ageuk.org.uk] Sent: Friday, March 22, 2013 11:25 AM To: Yuan, Rebecca; R help Subject: RE: error while

[R] predict.Arima error "'xreg' and 'newxreg' have different numbers of columns"

dep as the newxreg, but got error message(**): (**) Error in as.matrix(newxreg) %*% coefs : requires numeric/complex matrix/vector arguments In addition: Warning message: In chn.forecast[j] <- predict(fit, newxreg = t(indep), n.ahead = 1) : number of items t

[R] Assessing the fit of a nonlinear model to a new dataset

-statistic and p-values for each of the iterations using the trace function, but have not yet worked out how to do this. Any other suggestions are also welcome. Many thanks, Rebecca > model.wa <- nls(y ~ A*(x^B), start=list(A=107614,B=-0.415)) # create nls() > power model for WA data

[R] Can I open an existing pdf file and attach result to it?

Hello all, I would like to attach some results from R to an existing pdf file, can I do that through R? Thanks, Rebecca -- This message, and any attachments, is for the intended r...{{dropped:5

Re: [R] Assessing the fit of a nonlinear model to a new dataset

on of plausible models, and comparing my sum of squared residuals to those arising from that distribution. This seems to work. Thanks again. Cheers, Rebecca - Rebecca Lester Lecturer, Freshwater Ecology School of Life & Environmental Sciences Deakin University PO Box 423, Warrnambool, 3280

[R] The weak exogeneity test in R for the Error Correction Model?

google.com/books?id=ca5MkRbF3fYC Thanks very much! Cheers, Rebecca -- This message, and any attachments, is for the intended r...{{dropped:5}} __ R-help@r-project.org mailing

[R] sample size determination - comparing means for different size samples

Hmisc R package, specifically bsamsize to compare population proportions, but was wondering if there was an equivalent function for comparing means. -- Rebecca Greenblatt [[alternative HTML version deleted]] __ R-help@r-project.org mailing

[R] How to get a substring from a string

I get it out from the whole string? Thanks, Rebecca -- This message, and any attachments, is for the intended r...{{dropped:5}} __ R-help@r-project.org mailing list https://stat.ethz.c

Re: [R] How to get a substring from a string

This works! Thanks very much! From: Mohammad Tanvir Ahamed [mailto:mashra...@yahoo.com] Sent: Wednesday, January 15, 2014 9:48 AM To: Yuan, Rebecca; R help Subject: Re: [R] How to get a substring from a string Hi ! Try it st<-'yab~a+b+c' strsplit(st,"~")

[R] Format Sys.time()

Hello all, I know that the following code will give me the format of the Sys.time() till hour (%I), how can I get that till min? what shall I add between %I and %p? format(Sys.time(), "%Y_%m_%d_%I_%p") Thanks very much! Cheers

Re: [R] Format Sys.time()

Hello Petr, Thanks very much! This works! Cheers, Rebecca -Original Message- From: PIKAL Petr [mailto:petr.pi...@precheza.cz] Sent: Wednesday, January 15, 2014 10:01 AM To: Yuan, Rebecca; R help Subject: RE: Format Sys.time() Hi format(Sys.time(), "%Y_%m_%d_%I_%M%p&qu

Re: [R] How to get a substring from a string

Hello A.K., This gives more flexibility to substract the string, thanks very much! Cheers, Rebecca -Original Message- From: arun [mailto:smartpink...@yahoo.com] Sent: Wednesday, January 15, 2014 10:43 AM To: R help Cc: Yuan, Rebecca Subject: Re: [R] How to get a substring from a

[R] How to get the adjusted R squared from GLS() estimator

Hello, Although lm() gives a way to get the adjusted R squared by adjr2 <- summary(mdl)$adj.r.squared I cannot find a way to extract the adjusted R squared from gls(), any hint? Thanks, Rebecca -- This message, and

[R] Help required graphing factors with predicted model settings

Dear R wizards, Though I hate to do it after weeks of my code not working I need some help since I cant find an example which seems to work. I am trying to create a graph which show the probability of predation of a nest on one side (either 1 to 0) or (0% to 100%) on one side and grass height at t

Re: [R] Help required graphing factors with predicted model settings

;predict") : no applicable method for 'predict' applied to an object of class "mer" Can you see what I am doing wrong? What I am aiming for is a graph which looks like this. Thanks Rebecca On Thu, Oct 10, 2013 at 5:33 PM, Jim Lemon wrote: > On 10/10/2013 08:35 AM,

Re: [R] Help required graphing factors with predicted model settings

But with three lines for the three habitat types and grass length at the bottom On Thu, Oct 10, 2013 at 5:52 PM, Rebecca Stirnemann wrote: > Thanks Jim for helping, > > Your sample data actually looks like my dataset. The one I put up looks > strange for some reason so please igno

[R] How to altering the colour of confidence intervals in a lattice.

Hi R wizards, Can any one tell me how to code to remove the title panels in a lattice graph and to make my confidence intervals black instead of coloured? I have the following code: > library(effects) > library(lme4) > data(mao1, package="lme4") > fm1 <- lmer(frat ~ flandusenumb + ground.cover_lo

[R] Specify contrasts in R

problem using this contrast set-up in a mixed model? I am not sure how the 'block' effect would be interpreted for the mid-parent values, or indeed if it is sensible to include it as a random effect. Thank you in advance for your help, Rebecca Ross

[R] LMM and unequal variances

Hi all, I am running a LMM and have two fixed factors: Trt (with 2 levels) and Microhabitat (with 3 levels), and a randon factor: Group nested within Trt. I have unequal variances between the two Trt levels. Does this matter in a LMM? Thanks for your time Rebecca [[alternative HTML

[R] gsub wildcard

quot;states" with the period - not retaining the original value as I thought the wildcard character would. Any help is most appreciated! Thanks Rebecca [[alternative HTML version deleted]] __ R-help@r-project.org mailing list https://stat.ethz

Re: [R] gsub wildcard

d there are many more "[" that are unaffected. Rebecca On Mon, Aug 15, 2011 at 6:40 PM, Peter Langfelder < peter.langfel...@gmail.com> wrote: > On Mon, Aug 15, 2011 at 10:34 AM, Rebecca Gray wrote: > > Hello all, > > > > I have what I think is a simple qu

[R] gsub for numeric characters in string

; that is, I have not removed the names from the string so that now the state information is wrong. If anyone can help with the proper code to get rid of the names in tree while replacing the length with lengthplustate, I would appreciate it. I apologize if this was unclear. All the best,

[R] random sampling but with caveats!

Hi, I wonder if someone can help me. I have built a gam model to predict the presence of cold water corals and am now trying to evaluate my model by splitting my dataset into training/test datasets. In an ideal world I would use the sample() function to randomly select rows of data for me so fo

[R] plotting the ending point in a for loop

Hello, I would appreciate if someone could help me with this query. I would like to plot a line chart of all of the points in a "for" loop. I would also like to plot the final point with a symbol (to show where the random walk ends). Here is the code I am using: Brownian.fn <- function(mu, sigma

[R] R package fields: Thin-plate splines question

ction and subsequent method to separate the two plots (this would also require rearrangement of the data frame). However, I am stumped on the particulars. Does anyone have an idea of how to do this? I am including (hopefully reproducible) syntax for the two separate surfaces below. Thank y

[R] difference between foo$a[2] <- 1 and foo[2,"a"] <- 1

oo$a[1] <- 1 > foo$b[1] <- 2 > foo$a[2] <- 1 Error in `$<-.data.frame`(`*tmp*`, "a", value = c(1, 1)) : replacement has 2 rows, data has 1 > foo[2,"a"] <- 1 > foo a b 1 1 2 2 1 NA Thanks, Rebecca Hiller -- ETH Zürich Rebecca Hiller Institute of

[R] question on simple graph

I am having trouble generating a graph. I want to know the % of respondents who answered that they "strongly agree" or "agree" the "America owes R's ethnic group a better chance" (BTTRCHNC) and I want to organize it by racial group (RACESHRT). "BTTRCHNC" is organized ordinally from 1 through 5 wi

[R] Extracting variables from random effects formulas

k. Is there a way to get variable names out? Thanks in advance! Rebecca __ 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 provide commented

[R] Predicting a single observatio using LME

riable names to the left of | and to the right of |, the way one can with a normal formula.) Thanks in advance! Rebecca __ R-help@r-project.org mailing list https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read the posting guide http://www.R-proje

[R] "False convergence" in LME

way to get more information or to get the potentially wrong estimates from LME? (Also, the page in the NLMINB documentation, http://netlib.bell-labs.com/cm/cs/cstr/153.pdf, has errors in it, which makes it harder to check on what is happening.) Thank you in advance

[R] NaN as a parameter in NLMINB optimization

, and a few different (but very similar) likelihood functions. (They are complicated, but I can send them to someone if desired.) Is this something that can happen with NLMINB, perhaps because of a 0/0 in the gradient calculations? Or is it unique to my code? Thanks in advance! Rebecca

[R] LME prediction - object not subsettable?

the data for prediction back to the original dataset to which the lme model was fitted. Can anyone help? Thank you in advance! Rebecca __ R-help@r-project.org mailing list https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read the posting guide h

[R] Three questions about DSE

case with the attached code.) Is there a known reason why this sometimes occurs? Thank you in advance for all of your help! Rebecca -- Rebecca Sela Doctoral Candidate Statistics Group/IOMS Stern School of Business__ R-help@r-project.org mailing list

[R] DLM and matrices with 0 eigenvalues

as I can tell, the reason is that one of R's computed eigenvalues is very slightly negative (something like -5E-17). Is there a way to work around this? Thanks! Rebecca __ R-help@r-project.org mailing list https://stat.ethz.ch/mailman/listinfo/r

  1   2   >