#
par(mar = c(4.5, 4.5, 1, 1) + 0.1)
###
plot(D_mean,y,type="l",yaxt="n",lty=2,lwd=2,col="black",
ylab = list(expression(paste(dN/dlogD[agg]," ["*cm^-3*"]"))),
xlab = expression(paste(D[agg]," [nm]"
pe="l",yaxt="n",lty=2,lwd=2,col="black",
ylab = list(expression(paste(dN/dlogD[agg]," ["*cm^-3*"]"))), xlab =
expression(paste(D[agg]," [nm]")),
cex.lab=1.2
)
axis(2, mgp=c(0, 0.2, -2))
dev.off()
Charles Annis, P.E.
[EMAIL PROTECT
ery other major
bookseller, including the nearest university bookstore.
Charles Annis, P.E.
[EMAIL PROTECTED]
phone: 561-352-9699
eFax: 614-455-3265
http://www.StatisticalEngineering.com
-Original Message-
From: [EMAIL PROTECTED]
[mailto:[EMAIL PROTECTED] On Behalf Of Donatas G.
Sent: Wednesday,
Jessie:
How many pixels would you need to allocate for each of these 1000 parts? Is
that feasible?
Charles Annis, P.E.
[EMAIL PROTECTED]
phone: 561-352-9699
eFax: 614-455-3265
http://www.StatisticalEngineering.com
-Original Message-
From: [EMAIL PROTECTED]
[mailto:[EMAIL PROTECTED
e", box.lty=0)
You can look at the code for legend (?legend) and modify it, or use parts of
it for your own function.
Charles Annis, P.E.
[EMAIL PROTECTED]
phone: 561-352-9699
eFax: 614-455-3265
http://www.StatisticalEngineering.com
-Original Message-
From: [EMAIL PROTECTED]
Look at the definition for the transform. For example in the "car" package,
?box.cox
Then do the simple algebraic manipulations yourself.
Charles Annis, P.E.
[EMAIL PROTECTED]
phone: 561-352-9699
eFax: 614-455-3265
http://www.StatisticalEngineering.com
-Original Message
plot( x=rnorm(25, 0.5, 0.3), y=rnorm(25, 4, 1), xlim=c(0,1), ylim=c(2,7))
# ^^
for example
Charles Annis, P.E.
[EMAIL PROTECTED]
phone: 561-352-9699
eFax: 614-455-3265
http://www.StatisticalEngineering.com
-Original
You might consider having R do everything: R can read the Excel sheet, do
what needs to be done, and write the results to an Excel sheet.
Charles Annis, P.E.
[EMAIL PROTECTED]
phone: 561-352-9699
eFax: 614-455-3265
http://www.StatisticalEngineering.com
-Original Message-
From
's "insert" to point to
the new graphic.
Charles Annis, P.E.
[EMAIL PROTECTED]
phone: 561-352-9699
eFax: 614-455-3265
http://www.StatisticalEngineering.com
-Original Message-
From: [EMAIL PROTECTED]
[mailto:[EMAIL PROTECTED] On Behalf Of LL
Sent: Saturday, May 26, 2007 1:
Perhaps format is what you are looking for:
?format
Charles Annis, P.E.
[EMAIL PROTECTED]
phone: 561-352-9699
eFax: 614-455-3265
http://www.StatisticalEngineering.com
-Original Message-
From: [EMAIL PROTECTED]
[mailto:[EMAIL PROTECTED] On Behalf Of Pietrzykowski,
Matthew (GE
n the icon, as usual.) To do that, right-click on the program's icon and
choose "Run as administrator." Perhaps this will solve your RWinEdt
problems as well.
Best wishes!
Charles Annis, P.E.
[EMAIL PROTECTED]
phone: 561-352-9699
eFax: 614-455-3265
http://www.StatisticalEngineering
Monte Carlo Statistical
Methods, Springer, 1999
15. Venables and Ripley, Modern Applied Statistics with S, 4th ed.,
Springer, 2002
Charles Annis, P.E.
[EMAIL PROTECTED]
phone: 561-352-9699
eFax: 614-455-3265
http://www.StatisticalEngineering.com
-Original Message-
From: [EMAIL PROTEC
Sure.
rpois(n, lambda)
... will do it. But you should tell us something about how you want your
numbers to be distributed, since rpois() produces integers having a Poisson
distribution.
Charles Annis, P.E.
[EMAIL PROTECTED]
phone: 561-352-9699
eFax: 614-455-3265
http
Or you can right-click on the R icon and choose Run as administrator. That
way you won't alter the security settings and forget to re-set them. After
the packages are installed R will load in the usual way by clicking on the
icon.
Charles Annis, P.E.
[EMAIL PROTECTED]
phone: 561-352-9699
y as a mechanism for
implementing the statistical thinking you've acquired.
Charles Annis, P.E.
[EMAIL PROTECTED]
phone: 561-352-9699
eFax: 614-455-3265
http://www.StatisticalEngineering.com
-Original Message-
From: [EMAIL PROTECTED]
[mailto:[EMAIL PROTECTED] On Behalf Of Paul Lynch
Se
ty.
In my compiled html files I have my e-mail address and homepage, however
they do not have hyperlinks.
How can I include the hyperlinks? (Or how did I mess up and have R CMD
build not do that?)
Thanks.
Charles Annis, P.E.
[EMAIL PROTECTED]
phone: 561-352-9699
eFax: 614-455
Thank you Dan!
I dunno if I would have ever found that.
Thanks!
Charles Annis, P.E.
[EMAIL PROTECTED]
phone: 561-352-9699
eFax: 614-455-3265
http://www.StatisticalEngineering.com
-Original Message-
From: [EMAIL PROTECTED]
[mailto:[EMAIL PROTECTED] On Behalf Of [EMAIL PROTECTED
ensions<-matrix(c(strwidth(expression(theta),cex=5),strheight(expression(
theta),
cex=5)),nrow=1)
symbols(0.5,0.5
,rectangle=dimensions,bg='white',fg='white',add=TRUE,inches=FALSE)
text(0.5,0.5,expression(theta),cex=5)
~~~~~~~~~~~~~
Charles Annis, P.E.
[E
o difficult to figure
out), then install the libraries, then re-set the security levels.
Charles Annis, P.E.
[EMAIL PROTECTED]
phone: 561-352-9699
eFax: 614-455-3265
http://www.StatisticalEngineering.com
-Original Message-
From: [EMAIL PROTECTED]
[mailto:[EMAIL PROTECTED] On Behalf Of S
It might be helpful to those not familiar with Matlab to tell us what
function "persistent" does.
Charles Annis, P.E.
[EMAIL PROTECTED]
phone: 561-352-9699
eFax: 614-455-3265
http://www.StatisticalEngineering.com
-Original Message-
From: [EMAIL PROTECTED]
[mailto:[EMAIL PRO
, loaded the new package and there were my
new menus.
Thanks again.
Sheepishly,
Charles Annis, P.E.
[EMAIL PROTECTED]
phone: 561-352-9699
eFax: 614-455-3265
http://www.StatisticalEngineering.com
-Original Message-
From: [EMAIL PROTECTED]
[mailto:[EMAIL PROTECTED] On Behalf Of Charles Annis
build
--binary.
Should I just declare a namespace and use onLoad, or is there a mechanism to
use my own .First.lib?
Charles Annis, P.E.
[EMAIL PROTECTED]
phone: 561-352-9699
eFax: 614-455-3265
http://www.StatisticalEngineering.com
-Original Message-
From: Prof Brian Ripley [mailto:[
place.
Is there a way that I can use R CMD build --binary, with its zip file
output, but substitute my .First.lib function for the default?
Thanks for any help.
Charles Annis, P.E.
[EMAIL PROTECTED]
phone: 561-352-9699
eFax: 614-455-3265
http://www.StatisticalEngineering.com
.
Charles Annis, P.E.
[EMAIL PROTECTED]
phone: 561-352-9699
eFax: 614-455-3265
http://www.StatisticalEngineering.com
-Original Message-
From: [EMAIL PROTECTED]
[mailto:[EMAIL PROTECTED] On Behalf Of Benilton Carvalho
Sent: Thursday, November 16, 2006 12:38 AM
To: [EMAIL PROTECTED]
Cc: r-help
event.
Still struggling.
Charles Annis, P.E.
[EMAIL PROTECTED]
phone: 561-352-9699
eFax: 614-455-3265
http://www.StatisticalEngineering.com
-Original Message-
From: [EMAIL PROTECTED]
[mailto:[EMAIL PROTECTED] On Behalf Of Prof Brian Ripley
Sent: Thursday, November 16, 2006 4:44
: cannot open compressed file 'mypkg/DESCRIPTION'
But when I click on that file in the tarball it opens and shows me what I
expected.
I had hoped that I had weathered the hard part - building the package - but
I still need some help:
1) How do I get a zipped file, rather than a tarba
Thank you, Peter. The code provides a great example of packing the widget
nicely too.
Charles Annis, P.E.
[EMAIL PROTECTED]
phone: 561-352-9699
eFax: 614-455-3265
http://www.StatisticalEngineering.com
-Original Message-
From: [EMAIL PROTECTED] [mailto:[EMAIL PROTECTED] On Behalf Of
Jones _Practical Programming in Tcl
and Tk_ but would appreciate any help from fellow R users.
Many thanks.
Charles Annis, P.E.
[EMAIL PROTECTED]
phone: 561-352-9699
eFax: 614-455-3265
http://www.StatisticalEngineering.com
__
R-help@stat.math.ethz.ch ma
em too. Since the number of
available colors and available symbols differ, the color/symbol combinations
will be unique for the recycled usage.
Thanks for your counsel.
Charles Annis, P.E.
[EMAIL PROTECTED]
phone: 561-352-9699
eFax: 614-455-3265
http://www.Statisti
owth Class
1 1 4808 0.86212 0.0066964 S
3 3 2084 0.93579 0.4407900 S
Charles Annis, P.E.
[EMAIL PROTECTED]
phone: 561-352-9699
eFax: 614-455-3265
http://www.StatisticalEngineering.com
-Original Message-
From: [EMAIL PROTECTED]
[mailto:[EMAIL PROTECTED] On Behalf Of Chris Fon
D'oh!
You are right of course. How very obvious these things are on the other
side of the solution. I spent the entire weekend and completely missed it.
Thanks!
Charles Annis, P.E.
[EMAIL PROTECTED]
phone: 561-352-9699
eFax: 614-455-3265
http://www.StatisticalEngineerin
TRUE)
tcl.tk.example(PROBLEM = FALSE)
to turn the problem on and off.
I have also attached a tcl.tk.example.txt file since these long inclusions
within an R-help note often have the line feeds missing.
Thank you for your insights.
Charles Annis, P.E.
###
tcl.tk.example <-
2:
tkconfigure(answer.button1, variable=buttonValue, value="TRUE")
tkgrid(tklabel(frame.2, text="Question?yes "), answer.button1,
sticky="e")
^^^
I know this can not be as hard as I am making it.
Thanks
dows XP
machine, so the restarts are not onerous.
Many thanks to you all!
Charles Annis, P.E.
[EMAIL PROTECTED]
phone: 561-352-9699
eFax: 614-455-3265
http://www.StatisticalEngineering.com
-Original Message-
From: [EMAIL PROTECTED]
[mailto:[EMAIL PROTECTED] On Behalf Of Clint Bowm
es/R/R-2.4.0/library/file6fc97ac2/RColorBrewer/chtml/RColorBrewer.chm'
>
Any help would be appreciated. Thanks.
Charles Annis, P.E.
[EMAIL PROTECTED]
phone: 561-352-9699
eFax: 614-455-3265
http://www.StatisticalEngineering.com
__
R-help@stat.math.
sp1,Insp2,Insp3,Insp4,Insp5,Insp6,Insp7,Insp8,Insp9"
I can do that with a loop, but isn't there a more elegant way?
> result <- repeated.measures.FACTOR.names[[1]]
> for(i in 2:length(repeated.measures.FACTOR.names)) {
result <- paste(result, repeated.measures.FACTOR.names[[i]],
D'oh!
I've been using cat() but somehow never got the bigger picture.
Thanks!!!
Charles Annis, P.E.
[EMAIL PROTECTED]
phone: 561-352-9699
eFax: 614-455-3265
http://www.StatisticalEngineering.com
-Original Message-
From: Marc Schwartz [mailto:[EMAIL PROTECTED]
Sent:
quot;, "\"", "default.FACTOR.labels <- c('Probe1', 'Probe2', 'Probe3'))")
[1] "default.FACTOR.labels <- c(\"Probe1\", \"Probe2\", \"Probe3\"))"
^ ^ ^ ^ ^
ws XP on a Dell with 2Meg.
Thanks.
Charles Annis, P.E.
PS - Yes, I have upgraded to R version 2.4.0 (2006-10-03) but do not want to
change in the middle of a project.
[EMAIL PROTECTED]
phone: 561-352-9699
eFax: 614-455-3265
http://www.StatisticalEngineering.com
od surface.)
You may be looking for the resulting confidence bounds on the glm fit for
which I also have code that iteratively interrogates the loglikelihood
surface without plotting it.
If any of this is interesting, please send me a note so we won't clog the
bandwidth.
Charles Anni
}
Perfect!
Thanks!
Charles Annis, P.E.
[EMAIL PROTECTED]
phone: 561-352-9699
eFax: 614-455-3265
http://www.StatisticalEngineering.com
-Original Message-
From: [EMAIL PROTECTED]
[mailto:[EMAIL PROTECTED] On Behalf Of Charles Annis, P.E.
Sent: Thursday, September 28, 2006 10:13 PM
To:
your counsel.
Charles Annis, P.E.
[EMAIL PROTECTED]
phone: 561-352-9699
eFax: 614-455-3265
http://www.StatisticalEngineering.com
__
R-help@stat.math.ethz.ch mailing list
https://stat.ethz.ch/mailman/listinfo/r-help
PLEASE do read the posting guide
n(P(abs(phantom(i)*Z*phantom(i))>1.96,
H[1])==0.85,
P(abs(phantom(i)*Z*phantom(i))>1.96,H[0])==0.05), bty="n")
text(0, .2, quote(H[0]:~~mu[1]==mu[2]))
text(3, .2, quote(H[1]:~~mu[1]==mu[2]+delta))
#
Charles Annis, P.E.
[EMAIL PR
with R_ (Springer 2002). The
writing is relaxed and succinct, not condescending as some texts might
appear to a newcomer. It's just a good book.
Charles Annis, P.E.
[EMAIL PROTECTED]
phone: 561-352-9699
eFax: 614-455-3265
http://www.StatisticalEngineering.com
-Original Message
How close do you think it should be, given finite resolution with digital
computing?
> acos(0.5) - pi/3
[1] 2.220446e-16
Charles Annis, P.E.
[EMAIL PROTECTED]
phone: 561-352-9699
eFax: 614-455-3265
http://www.StatisticalEngineering.com
-Original Message-
From: [EMAIL PROTEC
Intercept) exped
5.2394254 0.9294818
Charles Annis, P.E.
[EMAIL PROTECTED]
phone: 561-352-9699
eFax: 614-455-3265
http://www.StatisticalEngineering.com
-Original Message-
From: [EMAIL PROTECTED]
[mailto:[EMAIL PROTECTED] On Behalf Of Dimitrios Rizopoulos
Sent: Sunday, September
Under Windows mine is located here
C:\Program Files\R\R-2.3.1\library\base\R
The file name, however is not .Rprofile, but rather Rprofile
Charles Annis, P.E.
[EMAIL PROTECTED]
phone: 561-352-9699
eFax: 614-455-3265
http://www.StatisticalEngineering.com
-Original Message-
From
Romain:
> a <- tempfile()
> cat(' document.write(screen.width) ;
', file=a)
> browseURL(a)
>
The object "a" was created, but no browser opened.
> ls()
[1] "a"
> a
[1] "C:\\DOCUME~1\\CHARLE~1\\LOCALS~1\\Temp\\RtmpRgWrqb\\file678418be&quo
that the same R code on this machine generates a
crowded jpeg on the older machine. My hope was to check the capability of
the current machine and make some modest adjustments to my R code to produce
an acceptable jpg on any machine running it.
Thanks.
Charles Annis, P.E.
[EMAIL PROTECTED]
phone:
Using sub rather than gsub appears to have no effect. The DOS window opens,
then closes having had no noticeable effect.
Charles Annis, P.E.
[EMAIL PROTECTED]
phone: 561-352-9699
eFax: 614-455-3265
http://www.StatisticalEngineering.com
-Original Message-
From: Gabor Grothendieck
splayed; nothing created. But the DOS window DID open and
something happened, but I don't know what.
Agin thanks for your patience.
Charles Annis, P.E.
[EMAIL PROTECTED]
phone: 561-352-9699
eFax: 614-455-3265
http://www.StatisticalEngineering.com
-Original Message-
From: Gabor Grothend
.vbs', intern = TRUE),
value = TRUE)))
The DOS window opened, some magic occurred in the blink of an eye, and the
DOS window closed. I haven't the foggiest idea what to do next since I can
see no evidence of having done anything.
Thanks for your patience.
Charles Annis, P.E.
[EMAIL
Charles Annis, P.E.
[EMAIL PROTECTED]
phone: 561-352-9699
eFax: 614-455-3265
http://www.StatisticalEngineering.com
-Original Message-
From: [EMAIL PROTECTED]
[mailto:[EMAIL PROTECTED] On Behalf Of bogdan romocea
Sent: Monday, August 28, 2006 10:55 AM
To: [EMAIL PROTECTED]
Cc: r-help
idn't anticipate it.
My code will be used on machines with varying graphics (and memory)
capacity. Is there a way I can check the native resolution of the machine
so that I can make adjustments to my code for the possible limitations of
the machine running it?
Thanks.
Charles Annis, P.E.
[EMAIL
Please visit the R site http://www.r-project.org/ and search the mailing
list for "paste expression" We discussed the topic recently. The "trick"
is that you don't paste expressions, you make an expression containing
paste.
Charles Annis, P.E.
[EMAIL PROTECTED]
phone
aster
than the competition with some system.time trials.
____________
Many thanks to all.
Charles Annis, P.E.
[EMAIL PROTECTED]
phone: 561-352-9699
eFax: 614-455-3265
http://www.StatisticalEngineering.com
-Original Message-
From: [EMAIL PROTECTED]
[mailto:
rst person to attempt this, and 2) I've made
this WAY harder than it is.
Would some kind soul please instruct me (and perhaps subsequent searchers)
how to convert the elements of a string into numbers?
Thank you.
Charles Annis, P.E.
[EMAIL PROTECTED]
phon
This'll work.
theta <- 2.1
plot(NA, xlim=c(0,1), ylim=c(0,1), xlab=bquote(theta == .(theta)),
ylab=bquote(theta == .(theta)), main=bquote(paste("Results for ",theta ==
.(theta
Charles Annis, P.E.
[EMAIL PROTECTED]
phone: 561-352-9699
eFax:
Wow! Exactly what I needed.
Thank you, Professor Ripley!
Charles Annis, P.E.
[EMAIL PROTECTED]
phone: 561-352-9699
eFax: 614-455-3265
http://www.StatisticalEngineering.com
-Original Message-
From: [EMAIL PROTECTED]
[mailto:[EMAIL PROTECTED] On Behalf Of Prof Brian Ripley
Sent
files,
but closing the file with dev.off() to produce the first graphic means that
I can't add to that and produce the subsequent graphic and file.
Is there a way to do this without generating the entire plot from the
beginning for each graphic?
Thanks
Charles Annis, P.E.
[EMAIL PROTEC
niversal way to find the name of a file at the end of a string
with "\\" separators.
Can I get there from here? (I've looked through previous R-help listing of
similar problems but if it's there, I missed it.)
Thanks.
Charles Annis, P.E.
[EMAIL PROTECTED]
phone:
has. I must wait for the calling routine to
end before I can use the new info, which is correctly place in the
workspace, in subsequent R routines.
Is there a way I can use the updated values in the same routine that created
the widget?
Thanks for your advice and patience.
Charles Annis, P.E
t is incomplete and should
also contain something like
Themes[i], sep=""
I am using R Version 2.3.1 (2006-06-01) on a DELL WindowsXP system with 2
gig RAM.
Many thanks for your counsel.
Charles Annis, P.E.
[EMAIL PROTECTED]
phone: 561-352-9699
eFax: 614-455-3265
http://www.Stati
osite.expression,lwd=3,type="l",cex.lab=1.3)
Charles Annis, P.E.
[EMAIL PROTECTED]
phone: 561-352-9699
eFax: 614-455-3265
http://www.StatisticalEngineering.com
-Original Message-
From: [EMAIL PROTECTED]
[mailto:[EMAIL PROTECTED] On Behalf Of Brian Phillip Weaver
Sent: Monday, June 1
Why create an enormous matrix? Why not read each company's info and
immediately write it to the file using write.csv( ... append = TRUE ...)?
^^^^^
Charles Annis, P.E.
[EMAIL PROTECTED]
phone: 561-352-9699
eFax: 614-455-3265
Department of Statistics at the University of South Carolina, but
as of May 2006 with the Department of Statistics at Texas A & M University,
http://www.stat.sc.edu/~west/javahtml/Histogram.html and
http://www.stat.tamu.edu/~west/
Charles Annis, P.E.
[EMAIL PROTECTED]
phone: 561-352-9699
eFax:
u're
happy with the defaults.) And the line end, lend=2, makes the ends square.
Best wishes.
Charles Annis, P.E.
[EMAIL PROTECTED]
phone: 561-352-9699
eFax: 614-455-3265
http://www.StatisticalEngineering.com
-Original Message-
From: [EMAIL PROTECTED]
[mailto:[EMAIL PROTECTED] On B
Many thanks to Duncan Murdoch, Thomas Lumley, Patrick Burns, and Seth
Falcon, for the illuminating advice, which will be found in the R-help
archives.
Charles Annis, P.E.
[EMAIL PROTECTED]
phone: 561-352-9699
eFax: 614-455-3265
http://www.StatisticalEngineering.com
-Original Message
d.
Thanks.
Charles Annis, P.E.
[EMAIL PROTECTED]
phone: 561-352-9699
eFax: 614-455-3265
http://www.StatisticalEngineering.com
__
R-help@stat.math.ethz.ch mailing list
https://stat.ethz.ch/mailman/listinfo/r-help
PLEASE do read the posting guide! ht
cex=0.8 (or whatever you like) as in
legend(0.3,0.4,c("name1","name2",...),col=1:20,lty=1:20, cex=0.8)
Charles Annis, P.E.
[EMAIL PROTECTED]
phone: 561-352-9699
eFax: 614-455-3265
http://www.StatisticalEngineering.com
-Original Message-
From: [EMAIL PROTE
Please give us an EXAMPLE of the loop you have in mind. (It's likely that
you can use simpler methods than a loop, but without an example we'd be
guessing.)
Charles Annis, P.E.
[EMAIL PROTECTED]
phone: 561-352-9699
eFax: 614-455-3265
http://www.StatisticalEngineering.com
---
="")==.(avar1)), adj=0)
#
Note that these look much the same but how they're constructed is rather
different.
Also, have a look at
?substitute
and
?bquote
Best Wishes,
Charles Annis, P.E.
[EMAIL PROTECTED]
phone: 561-352-9699
eFax: 614-455-3265
http://ww
Vielen Dank, Uwe!
I dunno how I missed it. Looking only at the demo(plotmath), I guess.
Thanks again!
Charles Annis, P.E.
[EMAIL PROTECTED]
phone: 561-352-9699
eFax: 614-455-3265
http://www.StatisticalEngineering.com
-Original Message-
From: [EMAIL PROTECTED]
[mailto:[EMAIL
(0,1.6), xlab="", ylab="", tck=-0.015)
text(1.6, 1, expression(paste("slope = ", frac(paste(" f( ", hat(theta),
" )"), hat(theta))), sep=""), adj=0)
text(2.25, 1.03, "\\pd", vfont=c("sans serif","plain")
hings up. Otherwise the file
will be there, but it won't be useable. So your final statement should be
dev.off()
Charles Annis, P.E.
[EMAIL PROTECTED]
phone: 561-352-9699
eFax: 614-455-3265
http://www.StatisticalEngineering.com
-Original Message-
From: [EMAIL PROTECTED]
[mai
7;time', 2=left censored, 3=interval censored.
Although unusual, the event indicator can be omitted, in
which case all subjects are assumed to have an event.
Charles Annis, P.E.
[EMAIL PROTECTED]
phone: 561-352-9699
eFax: 614-455-3265
http://www.StatisticalEngineering.
Go to the R site and look. http://www.r-project.org/
There has been some recent traffic on this topic.
Click on search: "Tobit"
Charles Annis, P.E.
[EMAIL PROTECTED]
phone: 561-352-9699
eFax: 614-455-3265
http://www.StatisticalEngineering.com
-Original Message-
Fr
You could substitute
lines(kernelgraf)
for your last line:
points(kernelgraf, xlab="Probability", xlim=c(0,100), ylim=c(0,.1),
col=rgb(0,0,1), main="")
See
?lines
Charles Annis, P.E.
[EMAIL PROTECTED]
phone: 561-352-9699
eFax: 614-455-3265
http://www.Statis
from your website.
Thanks.
Charles Annis, P.E.
[EMAIL PROTECTED]
phone: 561-352-9699
eFax: 614-455-3265
http://www.StatisticalEngineering.com
-Original Message-
From: [EMAIL PROTECTED]
[mailto:[EMAIL PROTECTED] On Behalf Of Kort, Eric
Sent: Thursday, January 26, 2006 2:39 PM
To: Thomas
Please don't forget engineering! (e.g. fatigue and reliability - censored
regression and survival; quantitative nondestructive evaluation - GLM)
Charles Annis, P.E.
[EMAIL PROTECTED]
phone: 561-352-9699
eFax: 614-455-3265
http://www.StatisticalEngineering.com
-Original Me
You can do something similar with Microsoft's browser but it isn't quite as
easy as Foxfire:
Right-click on the frame and choose Properties. Then highlight and copy the
URL and paste into the address window and click Go.
Then save the page.
Charles Annis, P.E.
[EMAIL PROTECTED]
t;- runif(100)
n.points <- length(my.data)
expected.cdf <- ((1:n.points)-0.5)/(n.points)
qqplot(my.data, expected.cdf, las=1)
# Use the "interocular trauma test" for goodness-of-fit:
my.lm <- lm(expected.cdf ~ sort(my.data))
abline(coef=coef(my.lm), lty=2)
Charles Annis, P.E.
[EMAI
Your problem sounds like it could be modeled with logistic regression
whereby the propensity for one result or another is "linked" to the factors
that control it. Logistic regressions are a special case of generalized
linear models. Look at ?glm
Charles Annis, P.E.
[EMAIL PROTEC
lot(data.df$x, data.df$y, log="xy")
This plot suggests that while Experiments B and C might have a similar
relationship between x and y, Experiment A differs.
Since I know nothing of the physical meaning of these observations I am
unqualified to comment further.
Best wishes.
Charl
nput.
Thanks to all!
Charles Annis, P.E.
[EMAIL PROTECTED]
phone: 561-352-9699
eFax: 614-455-3265
http://www.StatisticalEngineering.com
-Original Message-
From: Sundar Dorai-Raj [mailto:[EMAIL PROTECTED]
Sent: Thursday, November 10, 2005 5:00 PM
To: [EMAIL PROTECTED]
Cc: [EMAIL PR
cult as I have found it to be. Can anyone help?
Copious Thanks.
Charles Annis, P.E.
[EMAIL PROTECTED]
phone: 561-352-9699
eFax: 614-455-3265
http://www.StatisticalEngineering.com
__
R-help@stat.math.ethz.ch mailing list
https://stat.ethz.ch/mail
LL WinXP machine.
Charles Annis, P.E.
[EMAIL PROTECTED]
phone: 561-352-9699
eFax: 614-455-3265
http://www.StatisticalEngineering.com
-Original Message-
From: [EMAIL PROTECTED]
[mailto:[EMAIL PROTECTED] On Behalf Of jon butchar
Sent: Thursday, October 27, 2005 8:02 PM
To: Srinivas Iyy
Works well with both IE and Firefox on my 2 year old DELL WinXP machine.
Charles Annis, P.E.
[EMAIL PROTECTED]
phone: 561-352-9699
eFax: 614-455-3265
http://www.StatisticalEngineering.com
-Original Message-
From: [EMAIL PROTECTED]
[mailto:[EMAIL PROTECTED] On Behalf Of Gabor
Try
plot(Day, gene1)
lines(Day, gene2)
see
?lines
for more details.
Charles Annis, P.E.
[EMAIL PROTECTED]
phone: 561-352-9699
eFax: 614-455-3265
http://www.StatisticalEngineering.com
-Original Message-
From: [EMAIL PROTECTED]
[mailto:[EMAIL PROTECTED] On Behalf Of IAIN
" "formula" "terms"
"data"
[26] "offset""control" "method""contrasts"
"xlevels"
There are lots of ways to retrieve the parameter
ltitude, depend = tree, logi.mod = 1, type =
"dit", boxp = TRUE, rug = TRUE, las.h = 1)
which does not use the histograms but instead uses "dit plots" to provide a
helpful, visceral feel for the behavior of the observations.
Charles Annis, P.E.
[EMAIL PROTECTED]
phone: 56
with S_,
Fourth Edition. Springer, 2002. ISBN 0-387-95457-0.
http://www.stats.ox.ac.uk/pub/MASS4/. This book is more demanding and
covers a broad spectrum of contemporary statistical practice.
Charles Annis, P.E.
[EMAIL PROTECTED]
phone: 561-352-9699
eFax: 614-455-3265
http://www.Statistical
_Monte Carlo Statistical Methods_
by Christian P. Robert, George Casella
Springer, 2nd ed 2005
This book (I have edition 1) is a dandy. It will be rough sledding unless
you have a reasonable background in math stats but I think it is just what
you are looking for.
Charles Annis, P.E.
[EMAIL
?qqplot
Charles Annis, P.E.
[EMAIL PROTECTED]
phone: 561-352-9699
eFax: 614-455-3265
http://www.StatisticalEngineering.com
-Original Message-
From: [EMAIL PROTECTED]
[mailto:[EMAIL PROTECTED] On Behalf Of [EMAIL PROTECTED]
Sent: Friday, July 15, 2005 9:03 PM
To: r-help
x <- rep(NA, 3)
for (i in 1:length(x)){
x[i] <- ...
}
will do the job, but you may be able to take advantage of R's vectorization
and do what you want with no loop at all.
Charles Annis, P.E.
[EMAIL PROTECTED]
phone: 561-352-9699
eFax: 614-455-3265
http://www.StatisticalEngi
Thank you, Professor Ripley.
cbind(log(pr1$fit) - 1.96*pr1$se.fit/pr1$fit, log(pr1$fit) +
1.96*pr1$se.fit/pr1$fit)
... is precisely what had eluded me, self-evident though it appears after
you have illuminated the way.
Again, thank you.
Charles Annis, P.E.
[EMAIL PROTECTED]
phone: 561-352
I have tried everything I could
imagine I cannot reconcile the differences.
I believe that the confidence bounds for both models should agree. After
all, both calls to survreg() produce identical parameter estimates.
So I have missed something. Would some kind soul please point out my error?
Thank
.htm
I've also included an explanation of what is happening at each iteration.
Charles Annis, P.E.
[EMAIL PROTECTED]
phone: 561-352-9699
eFax: 614-455-3265
http://www.StatisticalEngineering.com
-Original Message-
From: [EMAIL PROTECTED]
[mailto:[EMAIL PROTECTED] On Behalf Of Paul
Vito:
Please plot your data:
y <- c(1, 1, 1, 1, 1, 1, 1, 0, 0, 0, 0, 0, 0)
x <- c(37, 35, 33, 40, 45, 41, 42, 20, 21, 25, 27, 29, 18)
plot(x, y)
You will see that ANY step function between 29 < x < 33 will describe these
observations perfectly.
Charles Annis, P.E.
[EMAIL PROT
ty statements for them.
Best wishes.
Charles Annis, P.E.
[EMAIL PROTECTED]
phone: 561-352-9699
eFax: 614-455-3265
http://www.StatisticalEngineering.com
-Original Message-
From: [EMAIL PROTECTED]
[mailto:[EMAIL PROTECTED] On Behalf Of
[EMAIL PROTECTED]
Sent: Thursday, January 20, 200
1 - 100 of 125 matches
Mail list logo
