Re: [ESS] Emacs, ESS and Rmarkdown: Is this the way compile is supposed to work?

For package cl is depreciated you can read this thread. I don't think this is causing any issues with markdown. https://www.reddit.com/r/emacs/comments/zslyem/package_cl_is_deprecated_debugging/ On Sun, Apr 23, 2023, 17:16 Kevin Zembower via ESS-help < ess-help@r-project.org> wrote: > Tyler,

[R] [R-pkgs] {peruse}

in-depth examples and extensive explanation of how to use the tool. CRAN link: https://CRAN.R-project.org/package=peruse<https://cran.r-project.org/package=peruse> Best, Jacob Goldsmith [[alternative HTML version deleted]] ___ R-pa

Re: [R] Seeking implementation of my algorithm 'spdspds' - a novel algorithm for solving Linear Programming Problems with O(L^1.5) computational complexity

" not having any access to the resources " You have internet (as is proven by sending this email) you even have R cli these days on mobile (I do anyway). Learn R , code, let me know how it turns out On Wed, Jun 10, 2020 at 8:46 PM Keshava PRASADa Halemane < k.prasa...@gmail.com> wrote: >

[R] [R-pkgs] Major Update to asciiSetupReader Package

please see here: https://cran.r-project.org/web/packages/asciiSetupReader/news/news.html Best, Jacob [[alternative HTML version deleted]] ___ R-packages mailing list r-packa...@r-project.org https://stat.ethz.ch/mailman/listinfo/r-packages

Re: [R] Proposed changes to the R Language

It has some good ideas, but R (my personal assesment) is not build for superspeed, but for superease of use.)) On Fri, Dec 28, 2018 at 3:01 PM Angus C Hewitt (DHHS) < angus.hew...@dhhs.vic.gov.au> wrote: > Hi Team, > > Please advise if there are any plans in the pipeline for change the R >

Re: [R] Codes to conduct network meta-analysis, gemtc package

ailing list where this probably belongs: https://stat.ethz.ch/mailman/listinfo/r-sig-meta-analysis Best, Wolfgang >-Original Message- >From: R-help [mailto:r-help-boun...@r-project.org] On Behalf Of Bert >Gunter >Sent: Thursday, 15 February, 2018 16:49 >To: shaju jacob >Cc:

[R] Codes to conduct network meta-analysis, gemtc package

Hi,I want codes that could help me do NMA from the data. I would appreciate if you could give a dummy table with which I could practice the codes that would help read the data and calculateOdds ratio (OR) Mean difference (MD) forest plot Density plot of posterior samples Estimates of ranks

Re: [R] Problems with "predict" function

Your messages about masking come from attaching your data set to the R session. In general, that is bad practice as it leads to confusing code. It is typically better to use the “data” argument in things like lm() to accomplish this task. As near as I can tell, your second set of predictions

Re: [R] Help with dates

Your dates have an incomplete year information with 34. R assumes that 00-68 are 2000 to 2068 and 69 to 99 are 1969 to 1999. See ?strptime and the details for %y. You can either append “19” to the start of your year variable to make it completely express the year or check if the date is in

[R] [R-pkgs] New CRAN Package Announcement: asciiSetupReader

the package with SPSS setup files. The process is the same as for SAS setup files. Please let me know if you if you find any bugs or problems in the package. https://github.com/jacobkap/asciiReader/issues Jacob [[alternative HTML version deleted

Re: [R] seek non-black box alternative to randomForest

Barry, This is mostly a mailing list about R - you have have more luck with statistical questions on www.stat.stackexchange.com. That said - the editor is wrong. The limitations of trees that random forests “solves” is overfitting. The mechanism by which a random forest classifier is built

[R] [R-pkgs] New package: fastDummies

tables or data.frames. This package is much faster than model.matrix() in handling data with greater than 1,000 rows. Please submit all question or issues here (https://github.com/jacobkap/ fastDummies/issues) Best, Jacob [[alternative HTML version d

Re: [R] Could not load Package 'rgl'

Christofer, This SO thread may be helpful: http://stackoverflow.com/questions/33634871/installing-rgl-package-in-r-mac-osx-el-captian-fixed On Apr 12, 2017, at 12:22 PM, Christofer Bogaso > wrote: Hi again, I could not load the

Re: [R] Too strange that I cannot install several packages

Bruce, `sample` doesn’t appear to be an R package. `resample` installed for me. `apply` and `sapply` aren’t packages either. `sample`, `apply` and `sapply` are all functions, however. > On Apr 10, 2017, at 12:08 PM, BR_email wrote: > > Hi Rers: > Is there anything I can

[R] problem with extract raster with polygon

ot; $`svn rev` [1] "71607" $language [1] "R" $version.string [1] "R version 3.3.2 (2016-10-31)" $nickname [1] "Sincere Pumpkin Patch" "Documentation for package ‘raster’ version 2.5-8" from Raster help pages Thanks for any help -- Jacob

[R] Validating Minitab's "Expanded Gage R Study" using R and lme4

I'm trying to validate the results of an "Expanded Gage R Study" in Minitab using R and lme4, but I can't get the numbers to match up in certain situations. I can't tell whether my model is wrong, my data is bad, or something else is going on. For instance, here's some data for which the results

[R] [R-pkgs] Announcing caesar

like to know more please visit the repository on GitHub ( https://github.com/jacobkap/caesar). Best, Jacob [[alternative HTML version deleted]] ___ R-packages mailing list r-packa...@r-project.org https://stat.ethz.ch/mailman/listinfo/r-packages

[R] extract minimal variables from model

-0.053701293 0.031099018 3 59 -0.041381647 -0.001334588 4 64 -0.029062001 -0.025298582 5 76 -0.016742356 -0.040792965 6 93 -0.004422710 -0.047817737 minimalvariablesfrommodel20161120(z1, ChickWeight)

Re: [R] package.skeleton fails

option like "--no-site-file" Thank you, Jacob -Original Message- From: Uwe Ligges [mailto:lig...@statistik.tu-dortmund.de] Sent: Sunday, August 28, 2016 7:28 AM To: Strunk, Jacob (DNR) <jacob.str...@dnr.wa.gov>; r-help@r-project.org Subject: Re: [R] package.skeleton fails

[R] package.skeleton fails

d the error appears happen within the subfunction '.fixPackageFileNames' Thanks for any assistance you might be able to provide. Jacob [[alternative HTML version deleted]] __ R-help@r-project.org mailing list -- To UNSUBSCRIBE and more, see http

[R] package.skeleton, environment argument causes error

zz", list0[wrong]) ok <- grepl("^[[:alnum:]]", list0) if (any(!ok)) list0[!ok] <- paste0("z", list0[!ok]) list1 <- tolower(list0) list2 <- make.unique(list1, sep = "_") changed <- (list2 != list1) list0[changed] <-

Re: [R] interaction between clusterMap(), read.csv() and try() - try does not catch error

Ok - got it, I can handle that. Thank you Luke! Jacob L Strunk ___ From: luke-tier...@uiowa.edu [luke-tier...@uiowa.edu] Sent: Monday, August 08, 2016 12:17 PM To: Strunk, Jacob (DNR) Cc: r-help@r-project.org Subject: Re: [R] interaction between clusterMap

[R] interaction between clusterMap(), read.csv() and try() - try does not catch error

which fails: try does > Not catch error > x0=clusterMap(clus0,function(...)try(read.csv(...)),c("c:/temp/badcsv.csv","c:/temp/goodcsv.csv"),SIMPLIFY=F) Error in checkForRemoteErrors(val) : one node produced an error: Error in read.table(file = file, header = header, sep

Re: [R] bug (?) with lapply / clusterMap / clusterApply etc

raded to 3.2.4revised and observed the same bug. I have googled for this issue and have not encountered any other individuals having a similar problem. I have attempted to reboot my machine without effect (aside from the obvious). Any suggestions would be great

[R] bug (?) with lapply / clusterMap / clusterApply etc

ine without effect (aside from the obvious). Any suggestions would be greatly appreciated! With regards, Jacob L Strunk Forest Biometrician (PhD), Statistician (MSc) and Data Munger __ R-help@r-project.org mailing list -- To UNSUBSCRIBE and more, see https://s

[R] confirm family==binomial and link==logistic

ctions?) rather than "deconstructing" an object (his term), as I do here. Accordingly, is there a smarter way to perform the check that I perform inside toyfunction? Thanks for any insight Jacob A. Wegelin Assistant Professor C. Kenneth and Dianne Wright Center for Clinical and Translation

[R] ggplot2 geom_boxplot limits

digest_0.6.8 proto_0.3-10 munsell_0.4.2 Jacob A. Wegelin __ 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 http://www.R-project.org/posting

Re: [R] ggplot2 geom_boxplot limits

*and* to zoom in? myplot+coord_cartesian(ylim=mylimits) # zooms in myplot+coord_cartesian(ylim=mylimits) + coord_flip() # flips but does not zoom myplot + coord_flip()+coord_cartesian(ylim=mylimits) # zooms but does not flip Jacob Wegelin Op 20-jul.-2015 20:29 schreef Jacob Wegelin jacobwege

[R] nlme splines model.frame.default error

(some=eks)),random=~1|ID)) print(mylme$call) Of course, with the workaround the resulting mylme does not know where the data came from. Why should a workaround be necessary? Is there something inherently misguided about the desired syntax (above)? Thanks for any insight Jacob Wegelin

Re: [R] nlme splines model.frame.default error

On 2015-06-02 Tue 14:20, Jacob Wegelin wrote: I want to use a specialized function to compute knots on the fly to pass to splines::bs within a call to nlme::lme. The specialized function will compute knots from the x-axis variable (the x argument in a call to splines::bs). The syntax works

[R] Probability distribution of fitted gaussian distribution.

sd 0.001996421 0.001403953 Reference http://cran.r-project.org/web/packages/fitdistrplus/fitdistrplus.pdf ~ Page 15 -- *Jacob Varughese* __ R-help@r-project.org mailing list -- To UNSUBSCRIBE and more, see https://stat.ethz.ch/mailman

[R] lme4 2 factor factorial model with random factors

Using lme4 how does one define a 2 factor factorial model with both factors being random? Specifically I am just trying to recreate the results from Montgomery's Design of Experiments book (7th edition), example 13.2. In this example there are 2 random factors and I want to include the

Re: [R] LogLikelihood of a Distribution Given Fixed Parameters

,1.66,0.25) and F(y) = plnorm(y,1.66,0.25). Note that instead of using 1-F(cens.time) you can use plnorm(cens.time,1.66,0.25,lower=TRUE) and that instead of taking logs explicitly you can set log=TRUE in the calls to dlnorm() and plnorm(). cheers, Rolf Turner On 25/04/14 07:27, Jacob Warren (RIT

[R] LogLikelihood of a Distribution Given Fixed Parameters

I'm trying to figure out if there is a way in R to get the loglikelihood of a distribution fit to a set of data where the parameter values are fixed. For example, I want to simulate data from a given alternate lognormal distribution and then I will fit it to a lognormal distribution with null

[R] ggplot2 Histogram with density curve

to a variable 'H' and use 'H' in place of '0.2' above why do I receive an error saying Error in eval(expr, envir, enclos) : object 'H' not found Jacob Trochelman OnPoint Contract Engineering (OCE) Services Engineering GE Aviation 513.243.6449 P Please consider the environment before printing

[R] scales percent precision

function (x) { x - round_any(x, precision(x)/100) str_c(comma(x * 100), %) } environment: 0x10c0f9350 find(precision) character(0) find(round_any) character(0) Thanks for any insights Jacob Wegelin sessionInfo() R version 2.15.3 (2013-03-01) Platform: x86_64-apple-darwin9.8.0/x86_64

Re: [R] scales percent precision

, Feb 27, 2014 at 8:49 AM, Jacob Wegelin jacobwege...@fastmail.fm wrote: scales::percent appears not to be documented. ?scales::percent_format where it tells you that it takes its argument, multiplies it by 100 and then attaches a percent sign to it. For most situations, the data should

Re: [R] scales percent precision

Incidentally, ?scales::percent brings up exactly the same text as ?scales::percent_format On 2014-02-27 Thu 14:47, Jacob Wegelin wrote: But percent_format() does not take the argument, multiply it by 100, and paste on a percent sign, as we see here: ?scales::percent_format percent_format

[R] plus/minus +/- in factor; not plotmath not expression

to do with the solution. Jacob A. Wegelin Assistant Professor Department of Biostatistics Virginia Commonwealth University 830 E. Main St., Seventh Floor P. O. Box 980032 Richmond VA 23298-0032 U.S.A. CTSA grant: UL1TR58 URL: http://www.people.vcu.edu/~jwegelin

[R] Plotting histogram of RT data

Hey there, I'm relatively new to R, and am currently working my way through some basic tutorials. I have a large data set that I've been able to import into the program. I'm using a script in which x=readdir (the directory). I am trying to create a histogram of the data in the 17th column, which

Re: [R] scope, lme, ns, nlme, splines

Dec 2012, Duncan Murdoch wrote: On 06/12/2012 2:28 PM, Jacob Wegelin wrote: I want to fit a series of lme() regression models that differ only in the degrees of freedom of a ns() spline. I want to use a wrapper function to do this. The models will be of the form y ~ ns(x, df=splineDF) where

[R] scope, lme, ns, nlme, splines

), data= junk) thislm cat(finished thislm\n) mystring-paste( thislme-lme( fixed= y ~ ns(x, df= , splineDF , ), random= ~ 1 | idvar , data= junk) , sep=) eval(parse(text=mystring)) thislme } KLUDGEfit() Thanks for any insight Jacob

Re: [R] Performing gage RR study in R w/more than 2 factors

On Mon, Nov 19, 2012, at 16:31, Bert Gunter wrote: I believe that you need to consult a local statistician, as there are likely way too many statistical issues here that you do not fully understand. Alternatively, try posting to a statistical list like stats.stackexchange.com, as I think most

Re: [R] Performing gage RR study in R w/more than 2 factors

On Mon, Nov 19, 2012, at 18:26, David Winsemius wrote: My guess is that you do not understand the meaning of a random factor. I certainly did not when I first encountered it. All my training had been with ordinary regression and analysis of variance. These are methods for what in mixed

[R] Performing gage RR study in R w/more than 2 factors

Hi everyone, I'm fairly new to R, and I don't have a background in statistics, so please bear with me. ;-) I'm dealing with 2^k factorial designs, and I was just wondering if there's any way to analyze more than two factors of a gage RR study in R. For example, Minitab has an expanded gage RR

[R] Logistic curve fitting with y values 1 (R version 2.15.2, OS is OS X 10.6.8)

Hello, I'm trying to fit a logistic curve to data but I'm having a hard time discovering how. Every tutorial I've come across either assumes the logistic curve has 0y1 or assumes I have multiple categories of data I simply have two vectors, a and b, of equal length with no missing data, and I

[R] lme(y ~ ns(x, df=splineDF)) error

, aes(x,y)) + geom_point(shape=1) + facet_wrap( ~ ID ) + geom_line( aes(x, fullfit), color=red) ) invisible(mylme) } Jacob A. Wegelin Assistant Professor Department of Biostatistics Virginia Commonwealth University 830 E. Main St., Seventh Floor P. O. Box 980032 Richmond VA

[R] lme( y ~ ns(x, df=splineDF)) error

) + geom_line( aes(x, fullfit), color=red) ) invisible(mylme) } Jacob A. Wegelin Assistant Professor Department of Biostatistics Virginia Commonwealth University 830 E. Main St., Seventh Floor P. O. Box 980032 Richmond VA 23298-0032 U.S.A. __ R

[R] Unbalanced Design Power Analysis

I have an unbalanced design I would like to run a power analysis on. What I have been able to find has pointed me to using the pwr.f2.test function as described below. My problem is that I don't know how to appropriately define the numerator and demoninator df. If someone can help here is some

[R] rJava load failure, 64-bit R, 64-bit Java, R 2.14.1 (works fine on R 2.12.2 64-bit, same computer)

issues: memory.limit() [1] 16381 Thank you! Jacob L Strunk [[alternative HTML version deleted]] __ R-help@r-project.org mailing list https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read the posting guide http://www.R-project.org/posting

[R] unable to load Hmisc in R 2.14.0

] survival_2.36-10 loaded via a namespace (and not attached): [1] cluster_1.14.1 grid_2.14.0lattice_0.20-0 psych_1.1.11 tools_2.14.0 Jacob A. Wegelin Assistant Professor Department of Biostatistics Virginia Commonwealth University 830 E. Main St., Seventh Floor P. O. Box 980032 Richmond VA

[R] plotting groups via density and different colors

I have a data set that looks like this: dene - data.frame(length = c(35,32,33,34,41,40,46,35,41,40,45,36,38,37,39,40,42,42,42,43,44), sex=c(1,1,1,1,2,2,2,1,2,2,2,1,2,2,2,2,2,2,2,2,2)) I would like to plot the density (frequency of occurrence) of each length class but I want to have different

[R] Retaining ID # with factor.scores procedure

included in the output are item responses, theta values, and standard errors for each individual. I need to be able to match the theta values with additional variables in another data set. Any help would be appreciated. Jacob T. Lutz, M.A. Doctoral Candidate School Psychology Program Ball State

[R] coxph() - unexpected result using Crawley's seedlings data (The R Book)

Rsquare= 0.022 (max possible= 0.993 ) Likelihood ratio test= 1.35 on 2 df, p=0.5097 Wald test= 1.32 on 2 df, p=0.5178 Score (logrank) test = 1.33 on 2 df, p=0.514 Anyone have an idea why this is occurring? Kind Regards Jacob __ R

Re: [R] coxph() - unexpected result using Crawley's seedlings data (The R Book)

Hi, sorry about that; here is the full output - data set, structure, model and result. Cheers Jacob seedlings cohort death gapsize status 1 September 7 0.5889 1 2 September 3 0.6869 1 3 September12 0.1397 1 4 September 1 0.1921 1 5 September

Re: [R] coxph() - unexpected result using Crawley's seedlings data (The R Book)

All, I rerun once again and managed to reproduce the results from the text book. Made no changes to the code. Could it be some problem with convergence? Anyhow, now it works! Cheers Jacob ps. I find The R Book very useful ds. 28 jun 2011 kl. 15.48 skrev Robert A LaBudde: Did you create

[R] Ghost values after subsetting

I am using subset to select the data I want to use for my analysis and find that after I subset my data frame on one column I get ghost values in the other columns. here is an example: table(data$Dags) 2008/04/12 2008/04/13 2008/04/16 2008/04/17 2008/04/19 2008/05/06 103

Re: [R] Ghost values after subsetting

It is a factor, thank you for pointing me in the right direction. Jacob On Thu, Jan 13, 2011 at 11:12, Sarah Goslee sarah.gos...@gmail.com wrote: Hi Jacob, You don't give us enough information to answer your question. Specifically, what is your dataframe? str(data) would be helpful

[R] ancova help

on 2 and 9864 DF, p-value: 2.2e-16 but this has not told me anything about where the differences in length that are attributable to sexual development lie. any suggestions? Thank you Jacob [[alternative HTML version deleted]] __ R-help@r

[R] ancova help

on 2 and 9864 DF, p-value: 2.2e-16 but this has not told me anything about where the differences in length that are attributable to sexual development lie. any suggestions? Thank you Jacob [[alternative HTML version deleted]] __ R-help@r

[R] plot xyz data in 2D

Hello, I want to make an 2D plot of .xyz data. So plot on x-axis and y-axis and use a color scale for the z. For every x-location along a cross section of the soil at several depths a resistivity must be displayed. This must result in a picture/graph which shows the resistivity for that cross

[R] ess - History Expansion

Hello R world, Can't seem to get history expansion working in R under emacs 23.2 with ess on mac osx. (Vincent Goulet's package). Even '!!' does not get me the previous command. What could be wrong? Thanks, vj __ R-help@r-project.org mailing

[R] ggplot2 boxplot: horizontal, univariate

important: How do I make a boxplot that is rotated (or transposed) from the default, so that the x axis carries the information? Thanks Jacob A. Wegelin Assistant Professor Department of Biostatistics Virginia Commonwealth University 730 East Broad Street Room 3006 P. O. Box 980032 Richmond VA

[R] fixed portion of lme4 model in newdata

=NEWDATA_frame) print(NEWDATA_matrix) Thanks for any insights Jacob A. Wegelin Assistant Professor Department of Biostatistics Virginia Commonwealth University 730 East Broad Street Room 3006 P. O. Box 980032 Richmond VA 23298-0032 APPENDIX sessionInfo() R version 2.10.1 (2009-12-14) i386

[R] R abrupt exit

on the topic so far. Any suggestions that the whole OS may have conflicts? I have only one copy of the following file on my computer system? /usr/lib/R/lib/libguide.so Thanks. jacob # options(KMP_DUPLICATE_LIB_OK

[R] fisher.test gives p1

] en_US.UTF-8/en_US.UTF-8/C/C/en_US.UTF-8/en_US.UTF-8 attached base packages: [1] stats graphics grDevices utils datasets methods base loaded via a namespace (and not attached): [1] tools_2.10.1 Thanks for any insight Jacob A. Wegelin Assistant Professor Department of Biostatistics

[R] match.call to obtain the name of a function

(that is, to subtract) the meat, lettuce, and mayonnaise. Thanks for any insights Jacob Wegelin __ 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

[R] conditioned xyplot, many y variables

as the axis tick labels. Thanks for any insights Jacob A. Wegelin Assistant Professor Department of Biostatistics Virginia Commonwealth University 730 East Broad Street Room 3006 P. O. Box 980032 Richmond VA 23298-0032 U.S.A.__ R-help@r-project.org

[R] parallel xyplot, different y scales, conditioned

. With relation=same, xyplot # automatically gets rid of that repetition. And we see an example of the fact that # scales$alternating=3 only works when relation=same. Thanks for any insight Jacob A. Wegelin Assistant Professor Department of Biostatistics Virginia Commonwealth University

[R] xyplot 3 panels 3 different Y variables

(lattice) xyplot( Y ~ Day | outcome, data=JUNKlong, layout=c(1,3), type=o, scales=list(x=list(alternating=3), y=list(relation=free, alternating=3)), ylab=) Am I making this harder than it needs to be? Thanks Jacob A. Wegelin Assistant Professor Department of Biostatistics Virginia Commonwealth

[R] wilcox.test loop through variable names

://tolstoy.newcastle.edu.au/R/e8/help/09/11/4520.html), if we are calling a regression function such as lm() we can avoid using paste(), as shown above. But is there a way to avoid paste() and eval() when one uses t.test() or wilcox.test()? Thanks Jacob A. Wegelin Department of Biostatistics Virginia

Re: [R] wilcox.test loop through variable names

On Sun, 15 Nov 2009 14:33 -0500, Jacob Wegelin jacobwege...@fastmail.fm wrote: Often I perform the same task on a series of variables in a dataframe, by looping through a character vector that holds the names and using paste(), eval(), and parse() inside the loop. For instance: rm

[R] loop through variable names

() for dichotomous. But the use of paste(), eval(), and parse() seems awkward. Is there a more elegant way to approach this? Thanks Jacob A. Wegelin Department of Biostatistics Virginia Commonwealth University Richmond VA 23298-0032 U.S.A. __ R-help@r-project.org

[R] superscript troubles

(bty = n,topright,paste(r2=,round(summary(lat_x)$r.squared,digits=3),, P=,round(coefficients(summary(lat_x))[2,4], digits=3))) Thank you Jacob -- Jacob Kasper http://twitter.com/Protect_Oceans 66°04' N 23°07' W Coastal Marine Management Master's Student University Centre of the Westfjords

[R] variable labels to accompany data.frame

: • If we try to add a value to a field in somedata$DATA which is not permitted by the permissible values listed for this field in somedata$VARIABLE, we get an error. Has anyone already thought this through, maybe defined a class and associated methods? Thanks Jacob A. Wegelin

[R] twoord.plot y lab size

(1985,2010), xlab=Year,ylab=Individuals, rylab=# Stations) -- Jacob Kasper http://twitter.com/Protect_Oceans 66°04' N 23°07' W Coastal Marine Management Master's Student University Centre of the Westfjords Sundstræti 1437 Devens Rd Ísafjörður, 400

[R] Mixing LaTeX and R Code in Loops and Functions in Sweave

I have a question about mixing LaTeX and R code in loops and functions in Sweave. Here's my problem: I want to do something like this: = # some R code for a loop for(i in 1:10) { @ My LaTex code here would describe what I do in the loop. For instance, in this loop print out the numbers one

[R] xyplot does not find variable in data

) ) ) , subset=Variety==Victory ) # But the following returns an error: xyplot( yield ~ nitro , data=Oats , scales=list( x=list( at=unique(nitro) ) ) ) Thanks for any insight Jacob A. Wegelin

[R] lattice auto.key drop unused levels

is the elegant or proper R thinking way to do this? That is, I want to get a key that only mentions the levels of Block that are used in the plot. Thanks Jacob A. Wegelin Assistant Professor Department of Biostatistics Virginia Commonwealth University 730 East Broad Street Room 3006 P. O. Box 980032

[R] princomp data input format

I am trying do to a PCA analysis using princomp. I get a result, but I wonder if I have the data in the correct format. My data contains many stations where fish were sampled as well as environmental information for each station (lat, lon, depth, temp and year). the format is like this: species 1|

[R] italics help in plot

Part of my script reads: speciesName - names(data)[i] plot(year,depth, xlab=Year, ylab=Depth(m),main=expression(italic(paste(speciesName))) ) Unfortunately, this just plots *speciesName *on my graph, not the name of the species in italics. Any suggestions on how to resolve this? Thank you Jacob

[R] symbols(x,y, circles=sqrt(N)) with lattice xyplot

How would I create the following plot using lattice? symbols( combPsummary$pastRate, combPsummary$finRate, circles=sqrt(combPsummary$N) ) The idea is to plot finRate vs pastRate using circles whose areas are proportional to the number of people in each group. The following attempt does not

Re: [R] Simulating points from GLM corresponding to new x-values

Hi Cliff, Excellent, just what I needed -- wish I had thought of that myself. Thanks a lot! Jacob On 13 Aug 2009, at 01:49, Clifford Long wrote: Hi Jacob, At the risk of embarrassing myself, I gave it a shot. I'll throw this out on the list, if for no other reason than perhaps someone

[R] Simulating points from GLM corresponding to new x-values

-values, which is not what I need. Ideally I whould like to simulate new values for GLM objects both with family=gaussian and with family=binomial. Thanks in advance, Jacob Jacob Nabe-Nielsen, PhD, MSc Scientist -- Section for Climate Effects

Re: [R] Simulating points from GLM corresponding to new x-values

() or simulate(), and it is also unclear whether simulate() can be applied to glms (with family=gaussian or binomial). Any suggestions for how to proceed? Jacob On 12 Aug 2009, at 13:11, Clifford Long wrote: Would the predict routine (using 'newdata') do what you need? Cliff Long Hollister

[R] lattice: control size of axis title and axis labels

Was: Re: Adjusting x/y text labels for a bwplot using cex.lab The purpose of this email is merely to explain how to control, separately, two different text sizes on the axes in an R plot, in traditional or lattice graphics. I did not find an explicit exposition of this, with examples, in a

[R] lattice: simultaneously control aspect outer whitespace

to be able to specify both the aspect ratio and the amount of whitespace outside the plot, and then let the software compute the necessary dimensions of the graphic device that are necessary to accomodate these specs. Thanks for any ideas Jacob A. Wegelin Assistant Professor Department of Biostatistics

[R] multiple lty on same panel in xyplot

into that kind of detail? (I use lattice, not traditional graphics, because the plot will eventually be conditioned on a third categorical variable.) Thanks for any suggestions. Jacob A. Wegelin Assistant Professor Department of Biostatistics Virginia Commonwealth University 730 East Broad Street

Re: [R] multiple lty on same panel in xyplot

On Wed, Aug 5, 2009 at 1:30 PM, Jacob Wegelin jacob.wege...@gmail.comwrote: I would like to use lattice graphics to plot multiple functions (or groups or subpopulations) on the same plot region, using different line types lty or colors col to distinguish the functions (or groups

Re: [R] multiple lty on same panel in xyplot

On Wed, 5 Aug 2009, Deepayan Sarkar wrote: On 8/5/09, Jacob Wegelin jacob.wege...@gmail.com wrote: I would like to use lattice graphics to plot multiple functions (or groups or subpopulations) on the same plot region, using different line types lty or colors col to distinguish the functions

[R] specify lattice black-and-white theme

-28 Jacob A. Wegelin Assistant Professor Department of Biostatistics Virginia Commonwealth University 730 East Broad Street Room 3006 P. O. Box 980032 Richmond VA 23298-0032 U.S.A. E-mail: jwege...@vcu.edu URL: http://www.people.vcu.edu/~jwegelin [[alternative HTML version deleted

Re: [R] session logging

of your session should do it? Jacob Wegelin wrote: Consider all the text that one sees on the console during an R session. Is there a way, within R, to make all this text--both the output and the messages--automatically get copied to a single text file, in addition to seeing it on the console

[R] session logging

grid_2.8.1 Thanks for any insights. Jacob A. Wegelin Assistant Professor Department of Biostatistics Virginia Commonwealth University 730 East Broad Street Room 3006 P. O. Box 980032 Richmond VA 23298-0032 U.S.A. E-mail: jwege...@vcu.edu URL: http://www.people.vcu.edu/~jwegelin

[R] session logging

for any insights. Jacob A. Wegelin Assistant Professor Department of Biostatistics Virginia Commonwealth University 730 East Broad Street Room 3006 P. O. Box 980032 Richmond VA 23298-0032 U.S.A. E-mail: jwege...@vcu.edu URL: http://www.people.vcu.edu/~jwegelin [[alternative HTML version

Re: [R] session logging

of Illinois fax: 217-244-6678Urbana, IL 61801 On Aug 3, 2009, at 3:01 PM, Jacob Wegelin wrote: Consider all the text that one sees on the console during an R session. Is there a way, within R, that all this text--i.e., both the output and the messages--can

[R] obtain names of variables and data from glm object

() function doesn't quite solve this problem: deparse(formula(myglm)) [1] y ~ x deparse(formula(myglm)[2]) [1] y() deparse(formula(myglm)[3]) [1] x() Ideally the elegant method would, in this example, return the character strings x, y, and DAT. Thanks for any insights. Jake Jacob A. Wegelin

[R] facet_grid problem: don't understand error message

trouble producing a reproducible example or solving the problem. Has anyone seen this before? Many thanks, Jacob Etches p - qplot(year,prop.excl,data=subset(excl,!is.na(ftf) exclusion %in% c(disab_inc_rec,family_death,fs_ch,ft2_ch,retir_inc_rec)), geom=c(line

[R] course in ecological data analysis

Dear all, For my PhD study I'm looking for relevant courses/workshops (short term) in ecological data anlysis with R in Europe. After 2 days searching I'm convinced that google is probably not the right medium to find this information. If anyone can help me I will be most grateful. Best regards

[R] ggplot2: specifying legend titles

I am trying to specify a legend title to be other than the variable name, but I find that the legend splits because scale_shape() takes effect but scale_colour() does not. Can someone spot my error? Here's some toy code that produces the problem on my system (R2.8.1, windows, today's

[R] mixed exponential distribution

Good morning, Is there anyway to do Mixed Exponential Distribution in R? I am trying to load some lag-weighted empirical survival distribution into R and run a mixed exponential on that data. Thanks, Jacob Fazekas Jacob Fazekas Assistant Actuary Auto-Owners Insurance Company 517-703-2543

  1   2   >