-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
" 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
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
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
://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-help
>Subject: Re: [R] Codes to conduct network
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 probab
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
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
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
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
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
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
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
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
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
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
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
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
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
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
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
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
. 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
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
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
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
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
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
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
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
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 interaction
ld appear that f(y) = dlnorm(y,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().
>
> c
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
param
h a variable binwidth? I.e. if I set my
binwidth 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 (OC
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
:
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
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
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
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
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.
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.
___
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
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
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"
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
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)
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),
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 mor
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
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
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
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
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
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
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
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
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
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
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
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
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
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
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
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 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
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
)
, scales=list(
ppp=list( alternating=3)
, y=list(
relation="sliced"
, alternating=3
, rot=0
, log=T
)
)
( alternating=3)
, y=list(
relation="sliced"
, alternating=3
, rot=0
, log=T
)
)
)
# In comparing these plots, we also
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
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
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
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.
(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
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
(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
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 o
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)
=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:
>
$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
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|
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
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
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 so
() 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?
-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
--
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
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
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://
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"
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)
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.)
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"
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
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
[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
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
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
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
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
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 - 100 of 109 matches
Mail list logo