Re: [R] Strange behavior using subset

(species 1, species 2, species 4, species 7, species 8)) cheers Ben Bolker __ R-help@stat.math.ethz.ch mailing list https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read the posting guide http://www.R-project.org/posting-guide.html

Re: [R] Interpreting the eigen value of a population matrix (2nd try)

/B6VHC-4JMM5XY-1/2/a698149bc3798c273766cfacdf40bba5 (accessed August 29, 2007). I've written a little bit of generic delta-method code, but I don't know if it's this generic. Ben Bolker __ R-help@stat.math.ethz.ch mailing list https://stat.ethz.ch

Re: [R] Optimization problem

(Hope this gets threaded properly. Sorry if it doesn't.) Gabor: Lac and Lacfac being the same is irrelevant, wouldn't produce NAs (but would produce something like a singular Hessian and maybe other problems) -- but they're not even specified in this model. The bottom line is that you

Re: [R] (no subject)

Series 326:283-293 (although the paper doesn't include any R code). cheers Ben Bolker __ R-help@stat.math.ethz.ch mailing list https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read the posting guide http://www.R-project.org/posting-guide.html

Re: [R] AIC and logLik for logistic regression in R and S-PLUS

Ben Bolker __ R-help@stat.math.ethz.ch 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.

[R] memory allocation glitches

.) and look for obvious problems? The package uses malloc/realloc rather than Calloc/Realloc -- does it make sense to go through the code replacing these all and see if that fixes the problem? cheers Ben Bolker __ R-help@stat.math.ethz.ch mailing list https

Re: [R] help with scatterplot3d

(0.19,0.295,0.24,0.22,0.24,0.28,0.205,0.26, 0.21,0.255,0.275,0.225,0.235,0.265,0.135) my.3dplot$points3d(x=xvals,y=yvals,z=zvals,col=blue,type=h,pch=16) Look carefully at ?seq and you may understand what you've done wrong ... Ben Bolker __ R-help

Re: [R] Artifacts in pdf() of image()

/R/Rhelp02a/archive/30581.html http://finzi.psych.upenn.edu/R/Rhelp02a/archive/71592.html does it happen when you print? On Ubuntu 7.04 I see a grid in evince; not in xpdf or acrobat reader; and in gv (but I can turn it off with 'a') Ben Bolker

Re: [R] question about logistic models (AIC)

log-likelihood is negative, therefore -2L + 2k is positive), although continuous probability densities or neglected normalization coefficients can lead to negative AICs -- but smaller (more negative, if AIC0) is still better. Ben Bolker __ R-help

Re: [R] t-distribution

Bill.Venables at csiro.au writes: for the upper tail: 1-pt(1.11, 9) [1] 0.1478873 wouldn't pt(1.11, 9, lower.tail=FALSE) be more accurate? __ R-help@stat.math.ethz.ch mailing list https://stat.ethz.ch/mailman/listinfo/r-help

Re: [R] stop criteria when L-BFGS-B needs finite values of 'fn' in optim

of 'fn' See ?try : e.g. myfun - try(optim(...)) if (class(myfun)==try-error) { ... whatever ... } else { ... success ... } or tryCatch this is FAQ 7.32 ... cheers Ben Bolker __ R-help@stat.math.ethz.ch mailing list

Re: [R] beta regressions in R

Walter R. Paczkowski dataanalytics at earthlink.net writes: Does anyone know of a package or function to do a beta regression? RSiteSearch(\\\beta regression\\\) (all those extra s seem necessary to get the equivalent of beta AND regression ...) betareg package on CRAN

Re: [R] the large dataset problem

Eric Doviak edoviak at earthlink.net writes: Dear useRs, I recently began a job at a very large and heavily bureaucratic organization. We're setting up a research office and statistical analysis will form the backbone of our work. We'll be working with large datasets such the SIPP as well

Re: [R] generating symmetric matrices

Gregory Gentlemen gregory_gentlemen at yahoo.ca writes: Greetings, I have a seemingly simple task which I have not been able to solve today. I want to construct a symmetric matrix of arbtriray size w/o using loops. The following I thought would do it: [snip] p - 6 Rmat - diag(p) vals

Re: [R] R codes for g-and-h distribution

filame uyaco filams0704 at yahoo.com writes: hi! I would like to ask help how to generate numbers from g-and-h distribution. This distribution is like normal distribution but span more of the kurtosis and skewness plane. Has R any package on how to generate them? Someone else

Re: [R] Constructing bar charts with standard error bars

John Zabroski wrote: On 7/25/07, Ben Bolker [EMAIL PROTECTED] wrote: Thanks a lot! I tried all three and they all seem very dependable. Also, I appreciate you rewriting my solution and adding elegance. Is there a way to extend the tick marks to the ylim values, such that the yscale ymax

Re: [R] Constructing bar charts with standard error bars

package library(gplots) with(testdata, barplot2(xbar,names.arg=group,main=a=4.0, xlab=Group,ylab=xbar,plot.ci=TRUE, ci.u=xbar+se,ci.l=xbar)) P.S. I hope you're not hoping to infer a statistically significant difference among these groups ... cheers Ben

Re: [R] Alternative to xyplot()?

lines (type r) are a little harder. How about: plot(y~x,col=grps) invisible(mapply(function(z,col) {abline(lm(y~x,data=z),col=col)}, split(data.frame(x,y),grps),1:3)) cheers Ben Bolker __ R-help@stat.math.ethz.ch mailing list https

Re: [R] how to use mle with a defined function

Lin Pan linpan1975 at yahoo.com writes: Hi all, I am trying to use mle() to find a self-defined function. Here is my function: test - function(a=0.1, b=0.1, c=0.001, e=0.2){ # omega is the known covariance matrix, Y is the response vector, X is the explanatory matrix odet =

Re: [R] Standard Probability Distributions.

David Barron mothsailor at googlemail.com writes: Try RSiteSearch to look for specific distributions. also try http://wiki.r-project.org/rwiki/doku.php?id=tips:stats-distri:0verviews=binomial __ R-help@stat.math.ethz.ch mailing list

Re: [R] Dominant eigenvector displayed as third (Marco Visser)

pair 6th, too). Ben Bolker __ R-help@stat.math.ethz.ch 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

Re: [R] LanguageR pvals.fnc error message

seem to do any fitting, just post-processing of the lmer object ... (by the way, that's the languageR package -- library(fortunes); fortune(Maechler.*package)) Ben Bolker __ R-help@stat.math.ethz.ch mailing list https://stat.ethz.ch/mailman/listinfo/r

Re: [R] Asteriscs in a plot to represent approximate size of p-values

Judith Flores juryef at yahoo.com writes: Hi, I need to place double and triple asterics (or stars) to highlight very low p-values. I am using points, for example: points(ssdx,ssdy,pch=8,cex=.9) but this allows me to place only one asterisc, how can I place 2 or 3 asteriscs?

Re: [R] extract index during execution of sapply

Christian Bieli christian.bieli at unibas.ch writes: Hi there During execution of sapply I want to extract the number of times the function given to supply has been executed. I came up with: mylist - list(a=3,b=6,c=9)

Re: [R] How to use mix to estimate the parameters for mixture gamma distribution?

=recover) to start diagnosing, or change print.level in your mix() call. But your best bet is to find a reproducible example for us to look at. cheers Ben Bolker __ R-help@stat.math.ethz.ch mailing list https://stat.ethz.ch/mailman/listinfo/r-help

Re: [R] how to obtain optimization with constraints

on the log scale (so that log(z)- 0.5*log x -0.5*log y) = 0), then use constrOptim to set linear inequality constraints. (This is assuming all the parameters are positive ...) Ben Bolker __ R-help@stat.math.ethz.ch mailing list https://stat.ethz.ch/mailman

Re: [R] How to compare GLM and GAM models

very much. Yuanchang Xie Since they can't possibly be nested I would suggest AIC. Ben Bolker __ R-help@stat.math.ethz.ch mailing list https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read the posting guide http://www.R-project.org

Re: [R] besselK

Robin Hankin r.hankin at noc.soton.ac.uk writes: Hello AFAIK, R has no capability for evaluating Bessel functions for complex arguments. Bessel functions for complex arguments are difficult to evaluate numerically. Some Bessel functions require cut lines or consideration of

Re: [R] specify constraints in maximum likelihood

adschai at optonline.net writes: Hi,I know only mle function but it seems that in mle one can only specify the bound of the unknowns forming the likelihood function. But I would like to specify something like, a = 2b or a = 2b where 'a' and 'b' could be my parameters in the likelihood

Re: [R] VGAM Pareto

livia yn19832 at msn.com writes: I would like to fit a Pareto Distribution and I am using the following codes fit=vglm(ycf1 ~ 1, pareto1(location=alpha), trace=TRUE, crit=c) fitted(fit) But the fitted values turn out to be the same for each observation. I guess the problem is with

Re: [R] pie chart in lattice - trellis class

-on-maps somewhere ... searching the R Graphics Gallery for pie also produces the hexbin pie plot (which doesn't use grid either ...) Ben Bolker __ R-help@stat.math.ethz.ch mailing list https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read

Re: [R] Eigen values

multivariate normal would be simplest.) Hints: mvrnorm (MASS), var, sample, eigen . Ben Bolker __ R-help@stat.math.ethz.ch mailing list https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read the posting guide http://www.R-project.org/posting

Re: [R] Problem with numerical integration and optimization with BFGS

Deepankar Basu basu.15 at osu.edu writes: For my model, the likelihood function for each observation is the sum of three integrals. The integrand in each of these integrals is of the following form: A*exp(B+C*x-D*x^2) (where D is positive) Being very lazy, I tried Mathematica's

Re: [R] problem with read.table

use quote= to disable the effect of ' cheers Ben Bolker __ R-help@stat.math.ethz.ch 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

Re: [R] comparing fit of cubic spline

AIC() on each of the fitted models and see what happens? In general R can compute log-likelihoods and AICs for most fitted model objects ... If it doesn't work, you can post your simple example to the list and get more advice then ... Ben Bolker __ R

Re: [R] plot problems

= as.numeric(names(t3)) lines(x2,t2,type=b,col=blue) lines(x3,t3,type=b,col=darkgreen) good luck, Ben Bolker __ R-help@stat.math.ethz.ch mailing list https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read the posting guide http://www.R-project.org

Re: [R] Anova Test

CrazyJoe keizer_61 at hotmail.com writes: I am very new to R. I am trying to perform an Anova Test and see if it differs or not. Basically, i have 4 tests and 1 control. Tester Test1 Test2 Test3 Test4 Control 20 25 1510 17 You can't make any

Re: [R] Simple question about function with glm

? [ hint: ?boxcox in the MASS package might be a better way to go ] Ben Bolker __ R-help@stat.math.ethz.ch mailing list https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read the posting guide http://www.R-project.org/posting-guide.html

Re: [R] spatial/ patches and gaps/ SADIE

, but I wonder if it would be worth writing an R interface ... Ben Bolker __ R-help@stat.math.ethz.ch 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

Re: [R] summary and min max

Sebastian P. Luque spluque at gmail.com writes: Hi, I came across a case where there's a discrepancy between minimum and maximum values reported by 'summary' and the 'min' and 'max' functions: By default summary only lists 3 significant digits ... see ?summary Ben Bolker

Re: [R] software comparison

than SAS. If anyone has a whole zoo of R versions lying around they could go back and try this code on them ... Ben Bolker scanstr = function(con,target,verbose=FALSE) { s = readLines(con,n=1) while (length(s)0 length(grep(target,s))==0) { if (verbose) cat(s,\n) s = readLines(con

Re: [R] software comparison

Thomas Lumley tlumley at u.washington.edu writes: On Mon, 16 Apr 2007, Philippe Grosjean wrote: Hello, [snip] No -- there is a new x.y.0 twice per year. Checking the r-announce archives shows that 2.0.0 came out in October 2004 and 1.9.1 in June 2004. Describing 1.9.1 as the

Re: [R] graphics - wireframe

=line will get you a wireframe (see ?rgl.material). http://wiki.r-project.org/rwiki/doku.php?id=tips:graphics-3d:graphics-3d good luck Ben Bolker __ [EMAIL PROTECTED] mailing list https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read

Re: [R] Simulated annealing using optim()

otalora at ugr.es writes: Finally, I tried to write my gr function to compute new random test configurations. Here I found a problem: In the simulated annealing algorithm, changes between succesive iterations during the optimization must be progressively smaller. The distance from the

Re: [R] Wikibooks

for the answers on the wiki.) Just my two cents -- and I've been delinquent in my wiki'ing recently too ... Ben Bolker __ R-help@stat.math.ethz.ch mailing list https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read the posting guide http://www.R

Re: [R] Problem adjusting x-labels with bargraphCI

you can use them in subsequent graphics calls. Ben Bolker __ R-help@stat.math.ethz.ch 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

Re: [R] fitting a mixed exponential distribution

like this: mixexplik - function(p,lambda1,lambda2) { log(p*dexp(x,lambda1) + (1-p)*dexp(x,lambda2)) } good luck, Ben Bolker __ R-help@stat.math.ethz.ch mailing list https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read the posting guide

Re: [R] Reg. strings and numeric data in matrix.

Mallika Veeramalai mallikav at burnham.org writes: I would like to know, is it possible to consider a matrix with some columns having numeric data and some other's with characters (strings) data? How do I get this type of data from a flat file. It's called a data frame. See the

Re: [R] No years() function?

Sérgio Nunes snunes at gmail.com writes: Hi, I'm trying to aggregate date values using the aggregate function. For example: aggregate(data,by=list(weekdays(LM),months(LM)),FUN=length) I would also like to aggregate by year but there seems to be no years() function. Should there be

Re: [R] Non : Confidence intervals for p**2 ??

(See section 6 of chapter 7 at http://www.zoo.ufl.edu/emdbook for more details if you like). Sorry to follow up on my own post, but the URL is wrong: http://www.zoo.ufl.edu/bolker/emdbook __ R-help@stat.math.ethz.ch mailing list

Re: [R] Non : Confidence intervals for p**2 ??

://www.zoo.ufl.edu/emdbook for more details if you like). Ben Bolker __ R-help@stat.math.ethz.ch 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

Re: [R] Angle of Bar Plot

)[as.numeric(factor(z[,2]))],density=5) I will point out that this is pretty ugly though ... Ben Bolker __ R-help@stat.math.ethz.ch mailing list https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read the posting guide http://www.R-project.org/posting

Re: [R] Google, hard disc drives and R

of boxplots, and why they broke the data into discrete age groups? Given that they had a relatively large data set (several percent of 100,000 disk drives), they could have done some cool visualization stuff ... Ben Bolker __ R-help@stat.math.ethz.ch

Re: [R] Any Limitations for the dataframes?

more specifics on the problem (approx. number of observations and elements per observation, and what you plan to try to do with them) you may get more useful feedback about whether what you plan to do is feasible and/or sensible. Ben Bolker signature.asc Description: OpenPGP digital signature

Re: [R] Wiki for Graphics tips for MacOS X

well enough even if it is not optimal. cheers Ben Bolker __ R-help@stat.math.ethz.ch 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

Re: [R] would you please navigate me?

a.eslami a.eslami at Mail.sbu.ac.ir writes: Hello My name is Aida Eslami. I am a M.S.c student of statistics at Shahid Beheshti University , Tehran, Iran. The subject of my thesis is Analysis of Masked Data. I have some problems in writing of my program (optimization). Would you

Re: [R] how to write randomforest in r

(preferably someone at your own institution) who can help you get started with R basics. Reading the posting guide wouldn't hurt either. good luck, Ben Bolker __ R-help@stat.math.ethz.ch mailing list https://stat.ethz.ch/mailman/listinfo/r-help

Re: [R] maximum likelihood, 1st and 2nd derivative

. See ?optim. Blatant plug: www.zoo.ufl.edu/bolker/emdbook/chap7A.pdf pp. 3-4 might be helpful too. Ben Bolker __ R-help@stat.math.ethz.ch mailing list https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read the posting guide http://www.R

Re: [R] nonparametric significance test for one sample

in the stats4 package ... Ben Bolker __ R-help@stat.math.ethz.ch 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

Re: [R] optimizing inverse student-t

and k are fixed scalars) Ben Bolker PS library(fortunes) fortune(Yoda) Evelyn Hall: I would like to know how (if) I can extract some of the information from the summary of my nlme. Simon Blomberg: This is R. There is no if. Only how. -- Evelyn Hall and Simon `Yoda' Blomberg R-help

Re: [R] Simulation with R

that this isn't a homework problem? Ben Bolker __ R-help@stat.math.ethz.ch 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

Re: [R] fmincon equivalent in R

to your problem, but you'll very much have to roll your own. (If any R-helpers out there would like to correct me or provide pointers to examples of Lagrange multiplier implementations in R, that would be great ...) Ben Bolker __ R-help

Re: [R] 3D histogramm

, as well as a variety of strong opinions about why you shouldn't plot your data as 3d histograms in the first place. Hints: demo(hist3d) in the rgl package example(scatterplot3d) in the scatterplot3d package Ben Bolker __ R-help

Re: [R] scaling y-axis to relative frequency in multiple hi stogram (multhist)

ahimsa campos-arceiz ahimsa at camposarceiz.com writes: Hi, I'm plotting a multiple histogram using the function multhist {package plotrix}, something like: library(plotrix) mh - list(rnorm(200, mean=200, sd=50), rnorm(200, mean=250, sd=50)) multhist(mh) In this graph y-axis

Re: [R] scaling y-axis to relative frequency in multiple histogram (multhist)

ahimsa campos-arceiz wrote: Dear Ben, thank you very much for your prompt reply. but I'm afraid I'm missing something: when I apply your function pmulthist I obtain exactly the same results as with multhist (with the y axis representing frequencies rather than a probability). I was

[R] low-variance warning in lmer

may not have looked in the right place ...) thanks Ben Bolker self-contained but unnecessarily complicated simulation code/demonstration: --- library(lme4) library(lattice) simfun - function(reefeff,ntreat=2,nreef=12, nreefpertreat=3

Re: [R] Stochastic SEIR model

see also: http://www.zoo.ufl.edu/bolker/eid/ __ R-help@stat.math.ethz.ch 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,

Re: [R] OPTIM--non finite finite different [13]

) instead of what you have written here. - A way of doing further debugging is to use cat() statements to print out parameter values and negative log-likelihoods for each function call. good luck Ben Bolker __ R-help@stat.math.ethz.ch mailing list

Re: [R] Fitting a survival curve

-squares fit through the data -- again, this is not quite right, but in the absence of information about sample size it's a little hard to know how to weight the different points. Hope that helps. Ben Bolker __ R-help@stat.math.ethz.ch mailing list

Re: [R] Generating a double-exponential jump diffusion process

be dir - ifelse(runif(nt)up,1,-1) jump - rexp(nt,scale=scale) cumsum(start+dir*jump) hope that helps Ben Bolker __ R-help@stat.math.ethz.ch mailing list https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read the posting guide http://www.R

Re: [R] Better y-axis labels, and x-axis scaling

that... warning, totally untested ... (the only way I can think of to get pretty and equally spaced ticks is to use ONLY powers of ten ...) Ben Bolker __ R-help@stat.math.ethz.ch mailing list https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read

Re: [R] graphics not find source

on cutting and pasting) * type version and cut and paste the results * tell us what distribution and version of Linux you're using You might get some help after that. Ben Bolker __ R-help@stat.math.ethz.ch mailing list https://stat.ethz.ch/mailman/listinfo

Re: [R] one problem about how to hold graphic with R

the command to hold the graphic with R? I remember with Matlab you may use hold on. Thanks. [snip] Is this what you were looking for? HTH G Actually, it sounds more like what yang wants is par(new=TRUE) Ben Bolker __ R-help

Re: [R] glmmPQL in 2.3.1

to interpret. use lmer (in the lme4 package) with the family argument set, to do real (or restricted) rather than quasi-likelihood estimation? Ben Bolker __ R-help@stat.math.ethz.ch mailing list https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do

Re: [R] Rmix package and different distributions

] Please note that this is completely untested! See Venables and Ripley's mixture-modeling example in MASS for more details of this kind of custom fitting (their example is a mixture of two normals, but it should still be helpful) Ben Bolker __ R-help

Re: [R] Measurements of 3000 criminals

-criminal subsample, and for a sample of 1000 Cambridge men. It would be an interesting challenge for a historian of science ... Ben Bolker __ R-help@stat.math.ethz.ch mailing list https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read the posting

Re: [R] Book for Maximum Likelihood Methods in R

/workshopnew/Documentation/sample%20likelihood.doc As is often the case, it would also help if you gave more detail about what you wanted to do -- probably the majority of modern statistical methods use maximum likelihood in one format or another ... cheers Ben Bolker

Re: [R] problem with mode of marginal distriubtion of rdiri chlet{gtools}

to the full posterior modes. (Think about the two-parameter normal distribution for example; for a sample {x_i}, the mode of the posterior variance is near sum((x_i-\bar x)^2)/n while the marginal mode is near sum((x_i-\bar x)^2)/(n-1) ...) [haven't checked this at all carefully.] Ben Bolker

Re: [R] [Q] How to fit data to a straight line

)) [1] 11.3995 I probably shouldn't have answered this -- please read the posting guide before you post again ... sincerely, Ben Bolker __ R-help@stat.math.ethz.ch mailing list https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read the posting

Re: [R] The old question: R vs MATLAB vs Mathematica vs ...

interfaces between them and R. cheers Ben Bolker __ R-help@stat.math.ethz.ch 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

Re: [R] MLE Methods

. Ben Bolker __ R-help@stat.math.ethz.ch 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] unexpected behavior of boxplot(x, notch=TRUE, log=y)

bogdan romocea br44114 at gmail.com writes: A function I've been using for a while returned a surprising [to me, given the data] error recently: Error in plot.window(xlim, ylim, log, asp, ...) : Logarithmic axis must have positive limits After some digging I realized what was

Re: [R] How can I generate this numbers

Ben Bolker __ R-help@stat.math.ethz.ch 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.

[R] linear terms within a nonlinear model

of implementing the linear structure within myfun(). Looked at Jim Lindsey's gnlm package but haven't yet been able to figure it out. Does anyone have any ideas or tips? thanks Ben Bolker __ R-help@stat.math.ethz.ch mailing list https://stat.ethz.ch

Re: [R] help on dirichlet distribution

and optimize with optim/nlmin/etc. Ben Bolker __ R-help@stat.math.ethz.ch 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

Re: [R] acos(0.5) == pi/3 FALSE

Ben Bolker __ R-help@stat.math.ethz.ch 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] rgl: exporting to pdf or png does not work

for rgl.snapshot() [although it would still be nice to have rgl.postscript() working ...] Ben Bolker __ R-help@stat.math.ethz.ch mailing list https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read the posting guide http://www.R-project.org/posting-guide.html

Re: [R] plotting curves

an example of the code you used and simultaneously to post examples of the output to a web page somewhere, if that's possible. I would just say, e.g., curve(x^3-2*x^2+x-1,from=-10,to=10) there is an n= argument that specifies how many points to use, but it doesn't seem necessary. Ben

Re: [R] maximizing a likelihood function containing an integral

is not reproducible/self-contained -- since I don't know what z1, z2, n1, rho are, I can't run it and see if that's the only problem. cheers Ben Bolker __ R-help@stat.math.ethz.ch mailing list https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do

Re: [R] Computing skewness and kurtosis with the moments package

the uncertainties in the skew and kurtosis are likely to be much larger. from _Numerical Recipes_: if the difference between n and n−1 ever matters to you, then you are probably up to no good anyway - e.g., trying to substantiate a questionable hypothesis with marginal data. cheers Ben Bolker

Re: [R] About the Skew Student distribution

this could be quite ugly). good luck Ben Bolker finance __ R-help@stat.math.ethz.ch 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

Re: [R] generating loglogistic distribution in R

a shifted, scaled version of the beta. Hence a + (b-a)*rbeta(n,shape1,shape2) cheers Ben Bolker __ R-help@stat.math.ethz.ch mailing list https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read the posting guide http://www.R-project.org

Re: [R] Handling realisations in geoRglm

directly would probably be more useful. Ben Bolker __ R-help@stat.math.ethz.ch 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

Re: [R] CircStats help

the sines and cosines of the individual angles and then taking the arc-tangent. With this size problem you can work through the details and understand how it all works ... Ben Bolker __ R-help@stat.math.ethz.ch mailing list https://stat.ethz.ch/mailman

Re: [R] Building a random walk vector

(runif(49,-2,2))) works. Ben Bolker __ R-help@stat.math.ethz.ch 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

Re: [R] maximum likelihood

methodology question than as an R question. Please read the posting guide ... Doing this without specifying any parametric forms would be tricky. You may be able to do this by searching for autoregressive model methods in RSiteSearch ... Ben Bolker

Re: [R] negative binomial lmer

neg binomial, which is Poisson-gamma) distribution, which might make this all work rather well in lmer: www.cefe.cnrs.fr/esp/TBElston_Parasitology2001.pdf cheers Ben Bolker __ R-help@stat.math.ethz.ch mailing list https://stat.ethz.ch/mailman

Re: [R] scatter plot with axes drawn on the same scale

Sundar Dorai-Raj sundar.dorai-raj at pdf.com writes: Try: plot(x, y, asp = 1) --sundar or eqscplot from the MASS package. __ R-help@stat.math.ethz.ch mailing list https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read the posting

Re: [R] String frequencies in rows

Mario Falchi mariofalchi at yahoo.com writes: I’m trying to evaluate the frequency of different strings in each row of a data.frame : INPUT: ID G1 G2 G3 G4 … GN 1 AA BB AB AB … Something like z - data[,-1] table(z,row(z)) ? Ben Bolker

Re: [R] Sweave and tth

environment! Ben Bolker __ R-help@stat.math.ethz.ch 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.

  1   2   3   >