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

2019-02-05 Thread Jacob
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

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

2016-03-22 Thread jacob
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

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

2016-03-23 Thread jacob
. I upgraded 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 greatly appreciated! With regards, Jac

[R] R abrupt exit

2010-04-02 Thread jacob
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_DUPLICAT

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

2017-04-21 Thread Jacob Kaplan
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

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

2017-09-15 Thread Jacob Kaplan
e 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 de

[R] package.skeleton, environment argument causes error

2016-08-18 Thread Jacob Strunk
list0[wrong] <- paste0("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 <- (l

[R] nlme splines model.frame.default error

2015-06-02 Thread Jacob Wegelin
o error: mylme<- with(toydat, lme(fixed=why~bs(eks,knots=customKnotsFn(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

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

2015-06-02 Thread Jacob Wegelin
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] ggplot2 geom_boxplot limits

2015-07-20 Thread Jacob Wegelin
tools_3.2.1 gtable_0.1.2 reshape2_1.4.1 [9] Rcpp_0.11.6 stringi_0.5-5 grid_3.2.1 stringr_1.0.0 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

Re: [R] ggplot2 geom_boxplot limits

2015-07-20 Thread Jacob Wegelin
tally *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 Wege

[R] Probability distribution of fitted gaussian distribution.

2015-01-03 Thread Jacob Varughese
28654032 $method [1] "mle" $sd mean 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 UNSU

Re: [R] Proposed changes to the R Language

2018-12-29 Thread jacob bogers
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 > langu

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

2016-02-12 Thread Jacob Wegelin
de functions?) 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 T

[R] extract minimal variables from model

2017-01-06 Thread Jacob Wegelin
1 -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, ChickWeig

[R] [R-pkgs] Announcing caesar

2017-01-23 Thread Jacob Kaplan
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] Validating Minitab's "Expanded Gage R&R Study" using R and lme4

2017-02-15 Thread Matt Jacob
I'm trying to validate the results of an "Expanded Gage R&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] problem with extract raster with polygon

2017-03-20 Thread Jacob Strunk
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] nested calls, variable scope

2008-07-17 Thread Jacob Wegelin
ice_0.17-6 loaded via a namespace (and not attached): [1] boot_1.2-32 grid_2.6.2 Thanks for any insight. Jacob A. Wegelin [EMAIL PROTECTED] Assistant Professor Department of Biostatistics Virginia Commonwealth University 730 East Broad Street Room 3006 P. O. Box 980032 Richmond VA 23298-0032

[R] nested calls, variable scope

2008-07-17 Thread Jacob Wegelin
0.999375-9 lattice_0.17-6 loaded via a namespace (and not attached): [1] boot_1.2-32 grid_2.6.2 > Thanks for any insight. Jacob A. Wegelin [EMAIL PROTECTED] Assistant Professor Department of Biostatistics Virginia Commonwealth University 730 East Broad Street Room 3006 P

[R] dput vs unclass to see what a factor really is composed of

2008-07-31 Thread Jacob Wegelin
I used read.dta() to read in a Stata 9 dataset to R. The "Sex01" variable takes on two values in Stata: 0 and 1, and it is labeled "M" and "F" respectively, analogous to an R factor. Thus, read.dta reads it in as a factor. Now, I wanted to see what this variable *really* is, in R. For instance, so

Re: [R] nested calls, variable scope

2008-08-01 Thread Jacob Wegelin
, > > Bert Gunter > Genentech > > -Original Message- > From: [EMAIL PROTECTED] [mailto:[EMAIL PROTECTED] > On > Behalf Of Jacob Wegelin > Sent: Thursday, July 17, 2008 1:58 PM > To: [EMAIL PROTECTED] > Subject: [R] nested calls, variable scope > > Below i

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

2009-08-08 Thread Jacob Wegelin
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 quick

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

2009-09-14 Thread Jacob Wegelin
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 reall

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

2010-02-26 Thread Jacob Wegelin
ordering a chicken sandwich and then telling the waitress to hold (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

[R] fisher.test gives p>1

2010-03-04 Thread Jacob Wegelin
R version 2.10.1 (2009-12-14) x86_64-apple-darwin9.8.0 locale: [1] 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 ins

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

2009-07-26 Thread Jacob Wegelin
would, in this example, return the character strings "x", "y", and "DAT". Thanks for any insights. Jake Jacob A. Wegelin Assistant Professor Department of Biostatistics Virginia Commonwealth University 730 East Broad Street Room 3006 P. O. Box 980032 Richmond V

[R] session logging

2009-08-03 Thread Jacob Wegelin
[1] lme4_0.999375-28 Matrix_0.999375-21 lattice_0.17-17 foreign_0.8-29 loaded via a namespace (and not attached): [1] boot_1.2-34 grid_2.8.1 Thanks for any insights. Jacob A. Wegelin Assistant Professor Department of Biostatistics Virginia Commonwealth University 730 East Broad Street

[R] session logging

2009-08-03 Thread Jacob Wegelin
me4_0.999375-28 Matrix_0.999375-21 lattice_0.17-17 foreign_0.8-29 loaded via a namespace (and not attached): [1] boot_1.2-34 grid_2.8.1 Thanks for any insights. Jacob A. Wegelin Assistant Professor Department of Biostatistics Virginia Commonwealth University 730 East Broad Stree

Re: [R] session logging

2009-08-03 Thread Jacob Wegelin
4558University 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, withi

Re: [R] session logging

2009-08-04 Thread Jacob Wegelin
t the beginning 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"

[R] multiple lty on same panel in xyplot

2009-08-05 Thread Jacob Wegelin
is kind of thing, with examples? I have not yet found an answer to this question in Paul Murrell's book. Does Deepayan Sarkar's _Lattice_ go into that kind of detail? (I use lattice, not traditional graphics, because the plot will eventually be conditioned on a third categorical variable.)

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

2009-08-05 Thread Jacob Wegelin
On Wed, Aug 5, 2009 at 1:30 PM, Jacob Wegelin 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 (or groups)

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

2009-08-05 Thread Jacob Wegelin
On Wed, 5 Aug 2009, Deepayan Sarkar wrote: > On 8/5/09, Jacob Wegelin 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"

[R] specify lattice black-and-white theme

2009-08-05 Thread Jacob Wegelin
tached): [1] Matrix_0.999375-21 boot_1.2-34grid_2.8.1 lme4_0.999375-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://

[R] lattice: simultaneously control aspect & outer whitespace

2009-08-07 Thread Jacob Wegelin
would be nice 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 Departm

[R] fixed portion of lme4 model in newdata

2010-04-15 Thread Jacob Wegelin
The following does not conform to NEWDATA but instead reverts to SomeDATA: NEWDATA_matrix<-model.matrix( object=MyModel, data=NEWDATA_frame) print(NEWDATA_matrix) Thanks for any insights Jacob A. Wegelin Assistant Professor Department of Biostatistics Virginia Commonwealth University 730 East Broad Str

[R] italics help in plot

2009-10-06 Thread Jacob Kasper
esolve this? Thank you Jacob [[alternative HTML version deleted]] __ R-help@r-project.org mailing list https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read the posting guide http://www.R-project.org/posting-guide.html and provide c

[R] princomp data input format

2009-10-07 Thread Jacob Kasper
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] lattice auto.key drop unused levels

2009-10-10 Thread Jacob Wegelin
$Block<-factor(as.character(junk$Block)) xyplot(yield~nitro, group=Block, data=junk, auto.key=T) What 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 P

Re: [R] lattice auto.key drop unused levels

2009-10-10 Thread Jacob Wegelin
=T) whereas dropping levels in the "group" argument does not create the desired effect xyplot(yield~nitro, subset=(Block=="I" | Block=="II"), data=Oats, group=Block [, drop=T], auto.key=T) Jacob Wegelin >  -Peter Ehlers > > Jacob Wegelin wrote: >

[R] xyplot does not find variable in data

2009-10-12 Thread Jacob Wegelin
x=list( at=unique(Oats$nitro) ) ) , subset=Variety=="Victory" ) # But the following returns an error: xyplot( yield ~ nitro , data=Oats , scales=list( x=list( at=unique(nitro)

[R] twoord.plot y lab size

2009-10-22 Thread Jacob Kasper
(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 14

[R] variable labels to accompany data.frame

2009-10-28 Thread Jacob Wegelin
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 Assistant Professor Department of Biostatistics Virginia Commonwealth University 730 East Broad Street Room 3006 P. O. Box 980032 Ri

[R] superscript troubles

2009-11-02 Thread Jacob Kasper
(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 & Marin

[R] loop through variable names

2009-11-11 Thread Jacob Wegelin
ing on its type--maybe lm() for continuous but wilcox.test() 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.

[R] wilcox.test loop through variable names

2009-11-15 Thread Jacob Wegelin
depending on its type. The use of paste(), eval(), and parse() seems awkward. As Gabor Grothendieck showed (http://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() an

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

2009-11-15 Thread Jacob Wegelin
On Sun, 15 Nov 2009 14:33 -0500, "Jacob Wegelin" 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 ins

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

2009-05-26 Thread Etches Jacob
aving 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&q

[R] course in ecological data analysis

2009-04-03 Thread Capelle, Jacob
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] xyplot 3 panels 3 different Y variables

2010-02-04 Thread Jacob Wegelin
t;]<-"Y.2" names(JUNK2)[names(JUNK2)=="Ratio"]<-"Y.3" JUNKlong<- reshape(JUNK2, dir="long", varying=2:4) JUNKlong$outcome<-factor( JUNKlong$time, levels=1:3, labels=c("Creatinine", "Estrogen", "Ratio") ) JUNKlong$t

[R] parallel xyplot, different y scales, conditioned

2010-02-05 Thread Jacob Wegelin
( alternating=3) , y=list( relation="sliced" , alternating=3 , rot=0 , log=T ) ) ) # In comparing these plots, we also

[R] conditioned xyplot, many y variables

2010-02-07 Thread Jacob Wegelin
) , scales=list( ppp=list( alternating=3) , y=list( relation="sliced" , alternating=3 , rot=0 , log=T ) )

[R] [R-pkgs] {peruse}

2021-03-10 Thread Jacob Goldsmith
-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

2020-06-11 Thread jacob bogers
" 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: > Frien

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

2012-12-19 Thread Jacob Wegelin
PRETTYNEWDATA$fitlme<-predict(thislme,level=0, newdata=PRETTYNEWDATA) thislme } Jake Wegelin On Thu, 6 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.

[R] Plotting histogram of RT data

2013-03-31 Thread Jacob Fake
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 i

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

2012-09-26 Thread Jacob Wegelin
ta=longdat)") thecommandstring eval(parse(text=thecommandstring)) } else { stop(paste("WhichApproach=", WhichApproach, " not valid.")) } mylme longdat$fullfit<-predict(mylme) library(ggplot2) print( ggplot( longdat, aes(x,y)) + geom_point(shape=1)

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

2012-11-09 Thread Jacob B.
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 0http://www.apsnet.org/EDCENTER/ADVANCED/TOPICS/ECOLOGYANDEPIDEMIOLOGYINR/DISEASEPROGRESS/Pages/NonlinearRegression.aspx, which assum

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

2012-11-19 Thread Matt Jacob
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 R&R study in R. For example, Minitab has an "expanded gage R&R"

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

2012-11-20 Thread Matt Jacob
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 mo

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

2012-11-20 Thread Matt Jacob
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 mi

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

2012-12-06 Thread Jacob Wegelin
sep="") eval(parse(text=mystring)) thislme } KLUDGEfit() Thanks for any insight 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] scales percent precision

2014-02-27 Thread Jacob Wegelin
path. Nor does round_any(). percent(0.010101) [1] "1.01%" percent function (x) { x <- round_any(x, precision(x)/100) str_c(comma(x * 100), "%") } find("precision") character(0) find("round_any") character(0) Thanks for any insights Jacob

Re: [R] scales percent precision

2014-02-27 Thread Jacob Wegelin
: Hi: On Thu, Feb 27, 2014 at 8:49 AM, Jacob Wegelin 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 be relative f

Re: [R] scales percent precision

2014-02-27 Thread Jacob Wegelin
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

2013-12-02 Thread Jacob Wegelin
ogle search initially led me to posts about expression() and plotmath, those eventually had nothing 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-0

[R] ggplot2 boxplot: horizontal, univariate

2010-06-18 Thread Jacob Wegelin
not from the y axis? But, more 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

[R] ess - History Expansion

2010-06-21 Thread Vinu Jacob
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 lis

[R] plot xyz data in 2D

2010-09-21 Thread Jacob Oosterwijk
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] ancova help

2010-10-19 Thread Jacob Kasper
iple R-squared: 0.1934,Adjusted R-squared: 0.1933 F-statistic: 1183 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 vers

[R] ancova help

2010-10-19 Thread Jacob Kasper
iple R-squared: 0.1934,Adjusted R-squared: 0.1933 F-statistic: 1183 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 vers

[R] "Ghost" values after subsetting

2011-01-13 Thread Jacob Kasper
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

2011-01-13 Thread Jacob Kasper
It is a factor, thank you for pointing me in the right direction. Jacob On Thu, Jan 13, 2011 at 11:12, Sarah Goslee 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] coxph() - unexpected result using Crawley's seedlings data (The R Book)

2011-06-28 Thread Jacob Brogren
trata(cohort)cohort=September2.0491 0.4880.379211.074 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 Any

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

2011-06-28 Thread Jacob Brogren
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

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

2011-06-28 Thread Jacob Brogren
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: > D

[R] plotting groups via density and different colors

2011-07-18 Thread Jacob Kasper
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 colo

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

2012-09-06 Thread Jacob Wegelin
andstring)) } else { stop(paste("WhichApproach=", WhichApproach, " not valid.")) } mylme longdat$fullfit<-predict(mylme) library(ggplot2) print( ggplot( longdat, aes(x,y)) + geom_point(shape=1) + facet_wrap( ~ ID ) + geom_line( aes(x, fullfit),

[R] unable to load Hmisc in R 2.14.0

2011-11-10 Thread Jacob Wegelin
tasets methods base other attached packages: [1] 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

[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)

2012-02-07 Thread Strunk, Jacob
of size 2.9 Gb Error: package/namespace load failed for 'rJava' There should not be any memory issues: > memory.limit() [1] 16381 Thank you! Jacob L Strunk [[alternative HTML version deleted]] __ R-help@r-project.org mailing list

[R] Help

2008-05-15 Thread Jacob, Mody
Dear All: I am new in Linux and I have inherited a Linux redhat system with "R Version 2.1.0 (2005-04-18)" installed. Can somebody help me with step by step instructions to download the new version and install it? Thanks mjacob The contents of this communication, including any attachments,

[R] Graphics after invoking R from the command line

2007-11-28 Thread Mithun Jacob
I've tried running graphics commands like plot by invoking R at the command line but the graphics window does not appear. I'm using R-2.6.0 on Windows XP and am using the cmd shell. Here's a sample session: R --slave --save --file=- x<-c(1,2,3,4) plot(x,x) This leads to nothing. I've found a way

Re: [R] Graphics after invoking R from the command line

2007-11-28 Thread Mithun Jacob
me that it's not possible to run graphics in batch mode? Thanks for all the help! Regards, Mithun On 11/28/07, Prof Brian Ripley <[EMAIL PROTECTED]> wrote: > On Wed, 28 Nov 2007, Mithun Jacob wrote: > > > I've tried running graphics commands like plot by invoking R

[R] ggplot2: specifying legend titles

2009-03-20 Thread Etches Jacob
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 CRAN

[R] xYplot, error bars, log scale

2008-02-05 Thread Jacob Etches
I'm new to Hmisc and trying to get the following to work, but if I un- comment the y-scale list (in order to get a log-scale for the hazard ratio), the error bars become strangely large. The dataframe is simply ODS output from TPHREG in SAS. Can someone point me towards what I'm sure is a

[R] ggplot2: can one have separate ylim for each facet?

2008-11-17 Thread Etches Jacob
lation="free"))) will give you a scale in each subplot with a range equal to the range of y within each subplot. Is this possible using ggplot2? qplot(x,y,data=w) + facet_grid(j~k) + ylim(-2,2) produces a plot with the same range in each subplot. Can the latt

[R] mixed exponential distribution

2008-12-10 Thread Fazekas, Jacob
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

[R] R in Windows 11

2024-10-09 Thread Jacob Williams
Hi, What is the minimum version of R supported on Windows 11? Thanks, Jacob __ 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 https://www.R-project.org

Re: [R] Help with dates

2017-12-28 Thread Simmering, Jacob E
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 the

Re: [R] Problems with "predict" function

2018-01-31 Thread Simmering, Jacob E
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 is

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

2017-04-10 Thread Simmering, Jacob E
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 check for as to wh

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

2017-04-12 Thread Simmering, Jacob E
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 mailto:bogaso.christo...@gmail.com>> wrote: Hi again, I could not load the 'rgl' package with below Erro

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

2017-05-30 Thread Simmering, Jacob E
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 i

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

2016-08-08 Thread Strunk, Jacob (DNR)
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] interaction between clusterMap(), read.csv() and try() - try does not catch error

2016-08-08 Thread Strunk, Jacob (DNR)
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] package.skeleton fails

2016-08-24 Thread Strunk, Jacob (DNR)
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

Re: [R] package.skeleton fails

2016-08-29 Thread Strunk, Jacob (DNR)
figuration 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) ; r-help@r-project.org Subject: Re: [R] package.skeleton fails Your code works for

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

2009-08-12 Thread Jacob Nabe-Nielsen
-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 --

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

2009-08-12 Thread Jacob Nabe-Nielsen
() 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?

  1   2   >