Re: [R] Help about nlminb function

to be worse (but not crazily worse) (it should give higher values). Ben Bolker __ 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, minimal, self-contained, reproducible code.

Re: [R] How to apitalize leading letters & else of personal names?

RockO gmail.com> writes: > > > Dear R world, > > Do you know about a function that would capitalize in the correct manner > first and family names? > I found in the cwhmisc only the CapLeading function, but it just does not do > the job, taking care only to capitalize the first letter of a wor

Re: [R] Integrate two function in R

Alaios yahoo.com> writes: > > Hello I have two function in R > like > g(x)=2x-3 > and s(x)=5x^2+2 > > and I want to find the integrafl of the g(x)*s(x) inside the interval of [a,b] > Analytically or numerically? It sounds like you want the answer analytically, in which case R can't do it

Re: [R] extracting IDs from topGO results

My guess is that you want to answer this question on the Bioconductor mailing list instead ... Ben Bolker __ 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/po

Re: [R] Integrate two function in R

Ben Bolker gmail.com> writes: > Analytically or numerically? > It sounds like you want the answer analytically, in which case R > can't do it, but it is an easy integral > > g(x)*s(x) = 10*x^3-15*x^2+4*x-6 > indefinite integral = 10/3*x^4 -15/3*x^3 + 4/2*x^2 -

Re: [R] Solution to differential equation

Mike Marchywka hotmail.com> writes: [snip] > > did you see my earlier post with link to wolfram integrator? Where i also > > requested anyone wanting to get rid of a copy of G&R Integral Tables to > > contact me off list since a dog really did eat mine? I think it came up > > with "F" or hype

Re: [R] Solution to differential equation

Ben Bolker gmail.com> writes: > Mike Marchywka hotmail.com> writes: [snip] > The gsl package has this function, apparently -- it agrees with > Mathematica/Wolfram Alpha's Hypergeometric2F1 for a single set of > inputs (2,3,4,0.5), although apparently the algori

Re: [R] dotchart for matrix data

David Winsemius comcast.net> writes: > > > On Dec 18, 2010, at 7:01 AM, e-letter wrote: > > > Readers, > > > > I am trying to use the function dotchart. The data is: > > > >> testdot > > category values1 values2 values3 values4 > > 1a 10 27 56 709 > > 2b

Re: [R] R.matlab memory use

Stefano Ghirlanda gmail.com> writes: > I am trying to load into R a MATLAB format file (actually, as saved by > octave). The file is about 300kB but R complains with a memory > allocation error: > > > library(Rcompression) > > library(R.matlab) > Loading required package: R.oo > Loading required

Re: [R] transposing panel data

Try the melt and cast functions in the reshape (or reshape2) packages (reshape2 may be faster). It may take a bit of work to get the syntax right (try it out on a small subset of your data). Ben Bolker __ R-help@r-project.org mailing list https://stat.eth

Re: [R] R.matlab memory use

d more memory than I > have (I have 3GB on this machine). > > Thanks, > Stefano > > On Sun, Dec 19, 2010 at 10:54 PM, Ben Bolker gmail.com> wrote: > > Stefano Ghirlanda gmail.com> writes: > > > >> I am trying to load into R a MATLAB format file (actually, a

Re: [R] logistic regression or not?

array chip yahoo.com> writes: [snip] > I can think of analyzing this data using glm() with the attached dataset: > > test<-read.table('test.txt',sep='\t') > fit<-glm(cbind(positive,total-positive)~treatment,test,family=binomial) > summary(fit) > anova(fit, test='Chisq') > First, is this still

Re: [R] Coding a new variable based on criteria in a dataset

RaoulD gmail.com> writes: > > > Hi, > > I'm a bit stuck and need some help with R code to code a variable F_R based > on a combination of conditions. > > The first condition would code F_R as "F" and would be based on the > min(Date) and Min(Time) for each combination of UniqueID & Reason. T

Re: [R] logistic regression or not?

package) is another way of handling analysis of proportions. > > -------- > *From:* Ben Bolker > *To:* r-h...@stat.math.ethz.ch > *Sent:* Tue, December 21, 2010 5:08:34 AM > *Subject:* Re: [R] logistic regression or not? > > array chip ya

Re: [R] Sweave for "big" data analysis

Duncan Murdoch gmail.com> writes: > As others have said, there are packages that provide caching. > > I haven't used them, because I like to keep my projects as > self-contained as possible: adding a dependency on one of those > packages is undesirable[1]. What I do in the case where there a

Re: [R] Retrieving Factors with Levels Ordered

Berwin A Turlach maths.uwa.edu.au> writes: > > G'day H.T. > > On Sat, 1 Jan 2011 00:41:10 -0500 (EST) > "H. T. Reynolds" udel.edu> wrote: > > > When I create a factor with labels in the order I want, write the > > data as a text file, and then retrieve them, the factor levels are no > > long

Re: [R] problem with postscript command

By the way, PostScript might not support transparency ... ? For your previous question (why does read.table not work the way you expected with respect to importing \u2264)? -- Reread the detailed answers you got previously that distinguish the difference between how the character is stored within R

Re: [R] Example in " Applied Spatial Analysis with R" giving error

markdx gmail.com> writes: > > > Example from page 29, Chapter 2.3 from "Applied Spatial Analysis with R" > > I am new to R...just trying to replicate the example from the book. > > > m <- matrix(c(0,0,1,1), ncol = 2, dimnames = list(NULL, + c ("min", > > "max"))) > Error in +c("min", "max") :

Re: [R] changing method of estimation in GLM

ATANU gmail.com> writes: > can anyone tell me how can i control the method of estimation (i.e. scoring > method or Newton raphson method) in glm and compute deviance function ? I don't think you can; you would have to write your own, although you can take advantage of the framework of the cur

Re: [R] Please, need help with a plot

Victor F Seabra uol.com.br> writes: > > >Please, I wonder if someone knows how to add the >less than or equal to symbol in the plot generated by the code below: >var1<-c('age <= 3','age <= 7','age <= 10','age <= 11','age <= 20','age <= >25','age <= 30','age <= 45','age <= 50','a

Re: [R] using " plot" with time series object - " axes = FALSE" option does not appear to work

} else axes <- TRUE This would make sense if axes were *not* in the explicit list of arguments to the function, and the default was supposed to be TRUE. Because axes *is* in the explicit list of arguments, this overrides it. I would cont

Re: [R] Greetings. I have a question with mixed beta regr ession model in nlme.

S and AD Model Builder (possibly among others). good luck Ben Bolker __ 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, minimal, self-contained, reproducible code.

Re: [R] how to keep keep matching column in output of merge

Nevil Amos gmail.com> writes: > How do I keep the linking column[s] in a merge()? > I need to use the values again in a further merge. simple reproducible example please? (e.g. make up a couple of 4-row datasets that show what you want, what you are trying, and what you are getting instead).

Re: [R] Cost-benefit/value for money analysis

Graham Smith gmail.com> writes: > > I assume this has a "proper" name, but I don't know what it is and wondered > if anyone knew of a package that might do the following, or something > similar. > > As an example, assume I have borrowed and read 10 books on R , and I have > subjectively given

Re: [R] How to save graphs out of ACF ?

bafin.de> writes: > > Hi, > > I want to save the autocorrelation plots resulting out of ACF (acf(ts)), not just by using the "Save as" > command in the R Gui but using some sort of code, which allows me to chose the format and the path. > Thank you, > > Mihai for example: a <- acf(runif(10

Re: [R] How to 'explode' a matrix

lt;- matrix(1:4,2,2) kronecker(m,matrix(1,2,2)) cheers Ben Bolker __ 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, minimal, self-contained, reproducible code.

Re: [R] Different LLRs on multinomial logit models in R and SPSS

[1] 5.22 > 204.2904-199.0659 [1] 5.2245 Ben Bolker __ 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, minimal, self-contained, reproducible code.

Re: [R] Problems with glht function for lme object

anord zooekol.lu.se> writes: > > > Dear all, > > I'm trying to make multiple comparisons for an lme-object. The data is for > an experiment on parental work load in birds, in which adults at different > sites were induced to work at one of three levels ('treat'; H, M, L). The > response is 'f

Re: [R] Anova with repeated measures for unbalanced design

Frodo Jedi yahoo.com> writes: > > Dear all, > I need an help because I am really not able to find over > internet a good example > in R to analyze an unbalanced table with Anova with > repeated measures. > For unbalanced table I mean that the questions are > not answered all by the same > nu

Re: [R] GLMM with lme4 and octopus behaviour

Samaritan gmail.com> writes: > [snip] You might want to ask follow-up questions on the R-sig-mixed-models list > At the most basic level, I'm testing the effect of sleep deprivation on > various behaviours (e.g. amount of time spent awake, amount of time spend > expressing difference texture

Re: [R] Using anova() with glmmPQL()

Toby Marthews ouce.ox.ac.uk> writes: > > Dear R HELP, > > ABOUT glmmPQL and the anova command. Here is an example of a > repeated-measures ANOVA focussing on the way > starling masses vary according to (i) roost situation and > (ii) time (two time points only). [snip] > lmeres=lme(fixed=stm

Re: [R] Help with logistic model with random effects in R

t's OK, but occasionally it could cause problems. * A rule of thumb is that you should not try to fit a model with more than k=N/10 parameters, where for binary data 'N' is the minimum of (number of successes, number of failures). * 6 groups is at the lower edge of feasibility for fi

Re: [R] MCMC object indexing

trix, e.g.: z <- cbind(a=runif(20),b=runif(20),c=runif(20)) library(coda) m <- mcmc(data=z) t(apply(m,2,quantile,c(0.025,0.975))) If that sort of thing doesn't work for you, please post an example of what doesn't work ... Ben Bolker

Re: [R] faster mvrnorm alternative

fantomas gmail.com> writes: > > > Hello, > does anybody know another faster function for random multivariate normal > variable simulation? I'm using mvrnorm, but as profiling shows, my algorithm > spends approximately 50 % in executing mvrnorm function. > Maybe some of you knows much faster fu

Re: [R] Ordering box plots

Stuart Luppescu ccsr.uchicago.edu> writes: > > I want box plots by group to display in order of increasing > mean (or median) of each group but can't seem to figure it out and > couldn't find anything on R-seek, either. > ?reorder ... ? __ R-help

Re: [R] Implementing step-wise linear regression

Tal Galili gmail.com> writes: > > Hello Troy. > > A tiny question (without answering your question), why did you choose to do > it this way instead of using > ?step > or > ?stepAIC > > ? [snip snip] > > My questions are: > > > > Should I be using var.test to run the F-test to decide which

Re: [R] Linear mixed model: question about t-values

Kostenko, Olga nioo.knaw.nl> writes: > I have a question about the output of linear mixed model fitted in R > using nlme package. In particular, what are the t-values that are given > in an output, how are they calculated and based on what test? I guess it > cannot be a simple Student t-test, oth

Re: [R] Help Derivate for Nonlinear Growth Models

utely need the derivative to use nls(), although it helps a lot. I would also suggest going to Google scholar and searching for 'nls "growth curve" Bates" to find some papers that have used this approach. If you need to post again, please read the posting guide and show

Re: [R] Checking package licences including dependencies?

Allan Engelhardt cybaea.com> writes: > > I only recently discovered options("available_packages_filters" = > list(add = TRUE, "license/FOSS")) [cf. help("available.packages", > package="utils") in R 2.10.0 or later] which goes nicely with my > options("checkPackageLicense" = TRUE) [new in

Re: [R] Axes=F and plotting dual y axes

Johnson, Cedrick W. cedrickjohnson.com> writes: > > Howdy. Been running into a bit of trouble with plotting. Seems that > axes=F is not "working". Whenever I plot (either a dataframe or xts/zoo > series) and I set axes=F along with xlab/ylab="" I still get the default > axes printed in my cha

Re: [R] read.delim()

Doran, Harold air.org> writes: > > Thank you, Phil. Unfortunately, there are quotes used properly elsewhere. > - Original Message - > From: Phil Spector stat.berkeley.edu> > To: Doran, Harold > Cc: r-help r-project.org r-project.org> > Sent: Wed Jul 28 18:29:32 2010 > Subject: Re: [R

Re: [R] Crash report: projection pursuit & predict

Horace Tso pgn.com> writes: > I'm using 2.11.1 under the 32-bit Windows XP. I've also tested the > codes under a 64-bit Ubuntu and the same > problem occurred. By the way, right before Rgui shuts down, the \ > pop-up window suggests it's a runtime > error in C++. > > Here are the codes that ca

Re: [R] Fry Plots

cassandra browne hotmail.com> writes: > > > Does anyone know how to do a fry plot in R? I have 600-800 points per image, and I really don't want to attempt > that manually. > > Thank you! > Cassie library(sos) ## you'll have to install this first findFn("fry plot") finds us

Re: [R] Is profile.mle flexible enough?

ent (I had a respectful disagreement with mle's author about whether this was a good idea or not, so I went off and wrote my own extensions ...) If I have a chance I will try to think about problem #1. My version (mle2) is, I think, more robust and featureful than stock mle, but it doe

Re: [R] Is profile.mle flexible enough?

Peter Dalgaard gmail.com> writes: > Issue #2 is stickier. I think I must say that the idelology of mle() is > that the user passes a likelihood function. If the likelihood function > depends on global data, and the user changes the global data, the user > deserves what he or she gets... It is pre

Re: [R] Help -normal distribution

n.a.s gmail.com> writes: > > Hi , can any one help in this problem > > According to some study, the height for Northern European adult males is > normally distributed with an average of 181 centimeter and a standard > deviation of 7.3 centimeter. Suppose such an adult male is randomly chosen. >

Re: [R] Including a text file (JAGS model) in an R package

Prof Brian Ripley stats.ox.ac.uk> writes: > > Here is how I do it: > > write(" > model { >for (year in 1:N) { > D[year] ~ dpois(mu[year]) > log(mu[year]) <- b[1] + step(year - changeyear) * b[2] >} >for (j in 1:2) { b[j] ~ dnorm(0.0, 1.0E-6) } >changeyear ~ dunif

Re: [R] Modelling poisson distribution with variance structure

Karen Moore tcd.ie> writes: > > I'm dealing with count data that's nested and has spatial dependence. > I ran a glmm in lmer with a random factor for nestedness. Spatial dependence > seems to have been accommodated by model. However I can't add a variance > strcuture to this model (to accommodat

Re: [R] pdf indifferent to lwd after some small number< 1 (e.g. 0.18)

Chuen Tan yahoo.com> writes: > > Hi, > > To plot a figure, I used lwd=0.18 in the function lines() to plot > very fine lines. However, using png() or > pdf() to save the figure, it appears that lwd has no effect on > pdf() beyond some small number (e.g. 0.18) and > the lines looks much thicker

Re: [R] Latex errors when checking package

Christian Hennig stats.ucl.ac.uk> writes: > LaTeX errors when creating PDF version. > This typically indicates Rd problems. > LaTeX errors found: > ! Font T1/ptm/m/n/10=ptmr8t at 10.0pt not loadable: Metric (TFM) file not > found I'm guessing you're running Ubuntu. If so,

Re: [R] Confidence Intervals for logistic regression

manchester.ac.uk> writes: > > On 07-Aug-10 09:29:41, Michael Bedward wrote: > > Thanks for that clarification Peter - much appreciated. > > > > Is there an R function that you'd recommend for calculating > > more valid CIs ? > > Michael > > It depends on what you want to mean by "more valid"!

Re: [R] Invalid plot type '1'

int that Erik was making was that if you have a question about a particular function, the help page for that function (in this case the results of ?plot or help("plot")) is the definitive documentation. Ben Bolker __ R-help@r-project.org

Re: [R] R Base Code

can probably adapt the existing pt/qt functions to get where you need, adjusting normalizing constants as appropriate ... PS: when asking questions about field-specific topics on the general R list, it's a good idea to give a little more background -- even references, in case someone is i

Re: [R] Hi!

leepama hanmail.net> writes: > I found some amazing package called "lqa" [enthusiasm snipped] > Can I use this package in linaer model(not GLM)? It's helpful to use informative subject lines in your posts to the list. A linear model is a special case of a generalized linear model with fami

Re: [R] Creating vectors

clips10 lancaster.ac.uk> writes: > I didn't really know what to post as the topic subject, but I have a vector, > for instance (2,2,4,6,2,4,4,6,8,6) [... snip to make gmane happy ...], > so my new vector would be (1,1,2,3,1,2,2,3,4,3). x <- c(2,2,4,6,2,4,4,6,8,6) as.numeric(factor(x)) [1] 1 1

Re: [R] Kalman filter

FMH yahoo.com> writes: > > Dear All, > > Could anyone give me a hand to suggest few packages in R to running Kalman > prediction and filtration ? Teach a person to fish ... install.packages("sos") library(sos) findFn("kalman") ## perhaps this could be added to the posting guide?

Re: [R] different outcomes of P values in SPSS and R

Leo Vorthoren nioo.knaw.nl> writes: > I have been using generalized linear models in SPSS 18, in order to build > models and to calculate the P values. When I was building models in Excel > (using the intercept and Bs from SPSS), I noticed that the graphs differed > from my expectations. When I r

Re: [R] detecting a key press

Martin Teicher hms.harvard.edu> writes: > I'm relatively new to r. I'd like to have a user respond by > pressing a 1 or a 2 and determining their choice and > the time of response. Previous postings have indicated that > keyboard responses can be processed using > scan and readline but both se

Re: [R] what is scaling in plot

elaine kuo gmail.com> writes: > Pls kindly advise what scaling is in plot. > Sometime it could be negative but sometimes it might be positive > .(I guess it is the proportion between the plot and the margin) Your question is unclear. Please give more context and/or details. __

Re: [R] different outcomes of P values in SPSS and R

gspot.com/2008/05/obtaining-same-anova-results-in-r-as-in.html > > Best Wishes, > > Martin H. Teicher > Dept of Psychiatry > McLean Hospital / Harvard Medical School > Belmont MA 02478 > > > On Aug 13, 2010, at 10:32 PM, Ben Bolker wrote: > >> Leo Vorthoren ni

Re: [R] how to display delta AIC of all models

to show the delta AIC of all models > rather than the final model with the least delta AIC. > > Please kindly advise if it is possible to exhibit all delta AIC of different > combinations of the independent variables. The MuMIn and AICcmodavg pac

Re: [R] AIC for model selection in Package â� �MuMInâ��

elaine kuo gmail.com> writes: [re: MuMIn package] > However, the criteria in this package for model selection are AICc > (second-order AIC) and QAIC (quasi-AIC). > Please kindly share if it is possible to carry AIC rather than the metrics > above using this package. > Thank you. I think you'r

Re: [R] AIC for model selection in Package â??Mu MInâ??

UE, eval = TRUE, rank = "AIC") > print(dd) > > Elaine > > On Mon, Aug 16, 2010 at 9:21 PM, Ben Bolker <mailto:bbol...@gmail.com>> wrote: > > elaine kuo http://elaine.kuo.tw> gmail.com > <http://gmail.com>> writes: >

[R] navel-gazing

month <- "2010-August" list <- "r-help" ##list <- "r-sig-ecology" ##list <- "r-sig-mixed-models" ## month <- "2010q3" n <- 50 baseurl <- "https://stat.ethz.ch/pipermail/"; library(RCurl) z <- getURL(paste(baseurl,list,"/",month,"/author.html",sep="")) zz <- strsplit(z,"")[[1]] namefun <- function(x

Re: [R] how to selection model by BIC

Xin__ Li mail.wbs.ac.uk> writes: > > Hi All: > the package "MuMIn" can be used to select the model based on AIC or AICc. > The code is as follows: > > data(Cement) > lm1 <- lm(y ~ ., data = Cement) > > dd <- dredge(lm1,rank="AIC") > print(dd) > > If I want to select the model by BIC, what cod

Re: [R] plotting functions of chi square

maiya gmail.com> writes: > Thanks, but that wasn't what I was going for. Like I said, I know how to do a > simple chi-square density plot with dchisq(). > > What I'm trying to do is chi-square / degrees of freedom. Hence > rchisq(10,i)/i). > > How do I do that with dchisq? By explicitl

Re: [R] GLMM random effects

> Natasha gmail.com> writes: For future reference, you may have better luck *either* on the r-sig-mixed-models help list (for advice about GLMMs) *or* on r-sig-ecology (for advice about ecological studies) [it's not considered polite to cross-post: try looking at the archives of both lists to d

Re: [R] (no subject)

what operating system you are using (which may well be important in dealing with this sort of bug), as well as the versions of R and frailtypack. good luck Ben Bolker __ R-help@r-project.org mailing list https://stat.ethz.ch/mailman/listinfo/r-help P

Re: [R] error heatmap and stack overflow

>michy har.mrc.ac.uk> writes: > > Hello, > Im trying to create a heatmap with a dataset (38 x 15037) but get the > error below: > > Error: protect(): protection stack overflow > Execution halted > > or > > Error: C stack usage is too close to the limit > Execution halted > > I tried to

Re: [R] lmer() causes segfault

cture, that's a clue -- it says it's some more particular characteristic of your data that triggers the problem) and (4) post to to *either* the R-sig-mac or the R-sig-mixed-models list, where the post is more likely to come to the attention of those who can help diagnose/fix ...

Re: [R] sendmailR-package-valid code needed

Barry Rowlingson lancaster.ac.uk> writes: > > On Mon, Aug 23, 2010 at 3:55 PM, Velappan Periasamy gmail.com> wrote: > > Hello Chris Campbell , > > I just posted my embryonic 'Rmail' package, which does a form of SMTP authentication (maybe not the version you want), to http://www.mathserv.mcm

Re: [R] R Package about Variable Selection for GLMM (Generalized Linear Mixed Model)?

be an appropriate place for the comment, "tell us what problem you are trying to solve, not just which methods you want to use", since sometimes people on the list are aware of different (and/or better, at least according to some criteria) methods for solving the

Re: [R] sendmailR-package-valid code needed

arrell Jr Professor and ChairmanSchool of Medicine > Department of Biostatistics Vanderbilt University > > On Mon, 23 Aug 2010, Ben Bolker wrote: > >> Barry Rowlingson lancaster.ac.uk> writes: >> >>> >>> On Mon, Aug 23, 20

Re: [R] 3D stariway plot

mauede alice.it alice.it> writes: > > > Please, is there an R function /package that allows for 3D stairway plots like the attached one ? > In addition, how can I overlay a parametric grid plot?? Not exactly, that I know of, but maybe you can adapt library(rgl) demo(hist3d) to do what

Re: [R] generate random numbers from a multivariate distribution with specified correlation matrix

rusers.sh gmail.com> writes: > rmvnorm()can be used to generate the random numbers from a multivariate > normal distribution with specified means and covariance matrix, but i want > to specify the correlation matrix instead of covariance matrix for the > multivariate > normal distribution. > Do

Re: [R] sendmailR-package-valid code needed

veepsirtt gmail.com> writes: > Warning: dependency ‘caTools’ is not available I'm not quite sure why this is complaining when you supposedly just loaded it ... do you really have it installed? > trying URL 'http://www.math.mcmaster.ca/~bolker/R/src/contrib/Rmail_1.0.zip' > Error in download.f

Re: [R] : Automatic Debugging

ArtemisIrvine hotmail.com> writes: Try options(error=recover) ? __ 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, mini

Re: [R] sendmailR-package-valid code needed

veepsirtt gmail.com> writes: > > > I could not install Rmail Package . > I got the following errors. > Then how to do.please > > library("caTools") > > install.packages("Rmail",contriburl="http://www.math.mcmaster.ca/~bolker/R/src/contrib";) try install.packages("caTools") install.packages

Re: [R] How to obtain seed after generating random number?

> On Aug 24, 2010, at 10:11 AM, Bogaso Christofer wrote: > > Dear all, I was doing an experiment to disprove some theory therefore > performing lot of random simulation. Goal is to show the audience that > although something has very rare chance to occur but it doesn't mean that > event would be

[R] several odfWeave questions

in the vignette.) My test files are posted at http://www.math.mcmaster.ca/~bolker/junk/odftest_in.odt and http://www.math.mcmaster.ca/~bolker/junk/odftest.odt thanks, Ben Bolker > > sessionInfo() > R version 2.11.1 (2010-05-31) i486-pc-linux-gnu locale: [1] LC_CTYPE=en_US

Re: [R] package MuMI

Marino Taussig De Bodonia, Agnese imperial.ac.uk> writes: > I am using the package "MuMI" to run all the possible combinations > of variables in my full model, and select > my best models. When I enter my variables in the original model I > write them like this > > lm(y~ a +b +c +a:b) > > How

Re: [R] package MuMIn

variables above: central$year:central$RC2 and central$RC2:central$year are both present, as are central$hunter:central$year and central$year:central$hunter . > > Can you tell me why, if they are the same thing, they are present twice and have different values? > > Thanks a lot, >

Re: [R] Problems loading rjags in R 2.11.0 under OS X 10.5.8

r version of JAGS ... Please update, if you can -- if not, or if that's going to be a major hassle, at least post the results of sessionInfo() Ben Bolker __ R-help@r-project.org mailing list https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read

Re: [R] Problems loading rjags in R 2.11.0 under OS X 10.5.8

Ben Bolker gmail.com> writes: > > Ruben Garcia Berasategui binus.ac.id> writes: > > > > > Dear list members, > > I'm trying to use JAGS 2.1.0 from within R ... > Sorry, didn't see the subject line (&

Re: [R] R.matlab package help

michael gmail.com> writes: > > C:\Users\FAN'S~1\AppData\Local\Temp\tpe2b4012b_f9ed_402d_af0f_f21ebd8116a6.mat > -V6 B I bet the problem is with the single quote (') in the path. __ R-help@r-project.org mailing list https://stat.ethz.ch/mailman/list

Re: [R] AIC using nls function

Bert Gunter gene.com> writes: > > John: > > 1. As always, and as requested (see posting guide), a small > reproducible example might help. Bert is right that things aren't well defined. However, AIC is still *widely* used for nonlinear models. For the sloppy folks among us, here are some u

Re: [R] R crashes when communicating with JAGS

Florian Weiler johnshopkins.it> writes: > First I create some data: > N <- 1000 > x <- rnorm(N, 0, 5) > > Then I specify a model in JAGS, storing it in the directory with the > extension .bug: > model { > for (i in 1:N) { > x[i] ~ dnorm(mu, tau) ## the model > } >

Re: [R] R crashes when communicating with JAGS

JAGS from a Windows shell (see the JAGS manual for examples), so that you can decide whether rjags is part of the problem, or whether it's a pure JAGS issue. > > On Tue, Aug 31, 2010 at 12:38 AM, Ben Bolker wrote: > >> Florian Weiler johnshopkins.it> writes: >> &

Re: [R] Cohen's f

Iasonas Lamprianou yahoo.com> writes: > > Dear all, > does anyone know how I can use any package in R to compute ANOVA's omega squared or Cohen's f? > > Dr. Iasonas Lamprianou > I don't know whether any of these do what you need, but have you browsed through the results of: library(sos) f

Re: [R] Retrieving predictions after fitting a model

alexander pswarayi ualberta.ca> writes: > > Hi there > > I am trying to retrieve the predicted values after fitting a model (dose > response model). How do I do it. What is the code for it? For estimated > parameters, one uses : summary(filename). And for the predicted values, what > do I use?

Re: [R] Using packages built in linux in windows

rajesh j gmail.com> writes: > I thought building a binary took care of that issue. How do I port a linux > library to windows? I think we need a little clarification here. When you say "a linux library", do you mean a chunk of compiled (non-R) code that the R package (sic) needs to use in or

Re: [R] problems with unique

n.vialma libero.it libero.it> writes: > > Hi everyone, > I have a question concerning the unique function. I'm using it with two > different data frame,which are similar in structure, but as in one case > unique > works properly in the other doesn't work and I don't understand why. > Someone

Re: [R] diamond/triangle "plot"

David Winsemius comcast.net> writes: > > > On Sep 1, 2010, at 7:50 PM, John K. Williams wrote: > > > Hi, does anyone have any idea how I might make a plot in R similar > > to this: > > http://assets.sbnation.com/assets/272612/BATTER-1B-PENA.png > > > > Specifically looking to plot 3 numbers

Re: [R] Linear models (lme4) - basic question

James Nead yahoo.com> writes: > > Sorry, forgot to mention that the processed data will be used as input for a > classification algorithm. So, I need to adjust for known effects before I can > use the data. > > > I am trying to adjust raw data for both fixed and mixed effects. > The data tha

Re: [R] Linear models (lme4) - basic question

() and residuals() for your example and see if it's doing what you want. > > > ---- > *From:* Ben Bolker > *To:* r-h...@stat.math.ethz.ch > *Sent:* Thu, September 2, 2010 2:06:47 PM > *Subject:* R

Re: [R] specify the covariance matrix for random effect

e, which are coded slightly more transparently. Followups should probably go to r-sig-mixed-mod...@r-project.org good luck Ben Bolker __ R-help@r-project.org mailing list https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read the postin

Re: [R] specify the covariance matrix for random effect

Ben Bolker gmail.com> writes: > > Qiu, Weiyu ualberta.ca> writes: > > > > > Hi, > > > > I'm doing a generalized linear mixed model, and I currently use an > > R function called "glmm". > > If you could get by with an AR

Re: [R] How can I fixe convergence=1 in optim

Peng, C gmail.com> writes: > > > To change the default maximum number of iterations (mxit =100 for derivative > based algorithm), add mxit = whatever number you want. > that's "maxit" i.e. optim(...,control=list(maxit=...)) __ R-help@r-project.

Re: [R] How can I fixe convergence=1 in optim

[forwarding back to r-help for archiving/further discussion] On 10-09-05 08:48 PM, Sally Luo wrote: > Prof. Bolker, > > Thanks for your reply and the helpful info. > > I still have a few questions. > > 1. I also tried to use different methods other than "BFGS" without > changing their def

Re: [R] R package to identify model

Ken Knoblauch inserm.fr> writes: > > Peng, C gmail.com> writes: > > > > > > > what is ESP package? Thanks. > > I've heard that It's only available over from a repository > accessible through a next-generation > wifi system call oui-ja. > > (Beware humor travels poorly over the internet >

<    5   6   7   8   9   10   11   12   13   14   >