On Fri, 8 Jun 2007, M. Jankowski wrote:
This is my first post requesting help to this mailing list. I am new
to R. My apologies for any breach in posting etiquette.
For future reference, telling us your version of R and exact OS would have
helped here. The R posting guide suggests showing
On Thu, 7 Jun 2007, Zack Weinberg wrote:
Given
D = data.frame(o=gl(2,1,4))
this works as I expected:
evalq(o, D)
[1] 1 2 1 2
Levels: 1 2
but neither of these does:
f - function(x, dat) evalq(x, dat)
f(o, D)
Error in eval(expr, envir, enclos) : object o not found
g - function(x,
Hi R-experts,
I have a data frame (A) , and a subset (B) of this data frame. I am trying
to create a new data frame which gives me all the rows of B, plus the 5th
next row(occuring in A). I have used the below code, but it gives me all 5
rows after the matching row. I only want the 5th.
Hi
I have a timeSeries object (X) with monthly returns. I want to display the
returns with a barplot, which I can fix easily. But my problem is labaling
the x-axis, if I use the positions from the timeseries It gets very messy. I
have tried rotating and changing the font size but it doesn't do
Hi R-users,
I need your help in the following problem. Suppose we have a regression
problem containing 25 predictor variables of 1000 individuals. I want to
divide the data matrix ( 1000 x 25 ) into two partitions for training (70%)
and testing(30%). For this reason, i sample 70% of data into
Hi,
you could use the sample function:
sample-sample(1:1000)
m.training-m[sample[1:700],]
m.test-m[sample[701:1000],]
Matthias
spime wrote:
Hi R-users,
I need your help in the following problem. Suppose we have a regression
problem containing 25 predictor variables of 1000
Dear list,
I have a very short question,
Suggest a dataframe of four columns.
df - data.frame(w,x,y,z)
I want this ordered the following way:
first by :x, decreasing = FALSE
and
secondly by: z, decreasing =TRUE
How can this be done ?
Thanks
Gunther
probably the function sort.data.frame() posted in R-help some time ago
can be useful; check:
RSiteSearch(sort.data.frame)
I hope it helps.
Best,
Dimitris
Dimitris Rizopoulos
Ph.D. Student
Biostatistical Centre
School of Public Health
Catholic University of Leuven
Address:
Tom.O schrieb:
Hi
I have a timeSeries object (X) with monthly returns. I want to display the
returns with a barplot, which I can fix easily. But my problem is labaling
the x-axis, if I use the positions from the timeseries It gets very messy. I
have tried rotating and changing the font size
Hi thanks for the respone, but cant you be more specific with your example. I
cant see that this will do the trick. What Im looking for is a function that
remembers each position but only displays every n'th date.
For example
positionReturns Disply Date
2003-01-31 1
Sorry I forgot the around the dates
x - c(dates(01/31/03),dates(06/30/07))
But I think your problem is the plot area.
You must first define the plot area with type =n for no plotting,
afterwards you could fill in the data.
I did this with times() but I am afraid the displayed dates/times
Hi all,
I have written a R-script under Windows using choose.dir. Now, I have
seen that this function is missing at MacOS. Does anybody know an
alternative?
Antje
__
R-help@stat.math.ethz.ch mailing list
https://stat.ethz.ch/mailman/listinfo/r-help
Sorry for double posting - was wrong e-mail adress , thougt this one
will run into Spam filter
__
R-help@stat.math.ethz.ch mailing list
https://stat.ethz.ch/mailman/listinfo/r-help
PLEASE do read the posting guide
Hi,
Can you provide examples of data formats that are problematic to read and
clean with R ?
The only problematic cases I have encountered were cases with multiline
and/or varying length records (optional information). Then, it is sometimes
a good idea to preprocess the data to present in a
On 6/8/07, Tom.O [EMAIL PROTECTED] wrote:
Hi
I have a timeSeries object (X) with monthly returns. I want to display the
returns with a barplot, which I can fix easily. But my problem is labaling
the x-axis, if I use the positions from the timeseries It gets very messy. I
have tried rotating
Steve,
You can do this with the package gstat. Look for ?krige of
?predict.gstat
Post further question on this topic on the R-sig-geo list. You'll get
more response.
Cheers,
Thierry
ir. Thierry Onkelinx
Instituut
Thanks for pointing at this. But you know, the user is writable. R is
installing Packages in /Documents/R/win-library which works fine so I
find it absolutely naturally that update should work as well. Especially
since when I install the packages it gets the latest version, library
loads this
Moving from S-plus to R I encountered many great features and a much
more stable system.
Currently, I am left with 2 problems that are handled differently:
1) I did lots of overplots in S-Plus using
par(new=T,xaxs='d',yaxs='d') to fix the axes
-What is the workaround in R ?
2) In S-Plus I could
On 08-Jun-07 08:27:21, Christophe Pallier wrote:
Hi,
Can you provide examples of data formats that are problematic
to read and clean with R ?
The only problematic cases I have encountered were cases with
multiline and/or varying length records (optional information).
Then, it is
see sort_df() in the reshape package
Justin BEM
Elève Ingénieur Statisticien Economiste
BP 294 Yaoundé.
Tél (00237)9597295.
- Message d'origine
De : Gunther Höning [EMAIL PROTECTED]
À : r-help@stat.math.ethz.ch
Envoyé le : Vendredi, 8 Juin 2007, 7h58mn 53s
Objet : [R] Sorting dataframe
Dear R users
I would like to ask a question regarding to icc (intraclass correlation) or many
biologists refer it to as repeatability. It is very useful to get icc for many
reasons and it is easy to do so from linear mixed-effects models and many
packages like psy, psychometric, aod and irr have
also try
active.sample-sample(1:1000,n=700)
active.df-thedf[active.sample, ]
test.df-thedf[-active.sample, ]
Justin BEM
Elève Ingénieur Statisticien Economiste
BP 294 Yaoundé.
Tél (00237)9597295.
- Message d'origine
De : Matthias Kirchner [EMAIL PROTECTED]
À :
Hi there,
two questions.
1) Is there any possibility to look up the help pages within R for more
complex combinations of character strings, for example Bayesian AND
regression but not necessarily Bayesian regression?
2) Is there a package/command that does fully Bayesian linear regression
(if
Dear R-users,
When I installed rattle package with the command:
install.packages(rattle, dependencies=TRUE), I got
Warning message:
Dependency 'Design' is not available
Is this warning serious? How to avoid this warning?
Thanks
_
Dr.Ruixin ZHU
maybe this page could give you some hints:
http://www.ats.ucla.edu/STAT/r/faq/sort.htm
Regards Knut
__
R-help@stat.math.ethz.ch mailing list
https://stat.ethz.ch/mailman/listinfo/r-help
PLEASE do read the posting guide
I was pointed at that my message might be considered as impolite. It was
not intended so. I was just trying to formulate that there should be
some improvement since the solutions offered were either not optimal for
me (disabling security features) or where not working (FAQ).
I apologize for any
Hi,
Is it possible to make a world map matrix where land values are set to 0 and
sea values to 1?
Cheers,
Antonio
[[alternative HTML version deleted]]
__
R-help@stat.math.ethz.ch mailing list
https://stat.ethz.ch/mailman/listinfo/r-help
Ruixin ZHU wrote:
Dear R-users,
When I installed rattle package with the command:
install.packages(rattle, dependencies=TRUE), I got
Warning message:
Dependency 'Design' is not available
Version of R? OS? Please do read the posting guide!
If R-2.5.0 under Windows: Design did not pass
try:
FiveDaysLater - A[match(rownames(B),rownames(A))+5,]
On 6/7/07, Alfonso Sammassimo [EMAIL PROTECTED] wrote:
Hi R-experts,
I have a data frame (A) , and a subset (B) of this data frame. I am trying
to create a new data frame which gives me all the rows of B, plus the 5th
next
On 6/8/2007 6:52 AM, Antonio Rodríguez wrote:
Hi,
Is it possible to make a world map matrix where land values are set to 0 and
sea values to 1?
It's not hard to produce a bitmap of a world map with the maps package,
and then some image manipulation functions could convert it to 0's and
On 6/7/07, Robert Wilkins [EMAIL PROTECTED] wrote:
As noted on the R-project web site itself ( www.r-project.org -
Manuals - R Data Import/Export ), it can be cumbersome to prepare
messy and dirty data for analysis with the R tool itself. I've also
seen at least one S programming book (one of
Dear R-user,
Could anybody tell me of the key difference between data mining and text
mining?
Please make a list for packages about data/text mining.
And give me an example of text mining with R (any relating materials
will be highly appreciated), because a vignette written by Ingo Feinerer
Tim Bergsma said the following on 6/8/2007 5:57 AM:
Suppose I have a list of logicals, such as returned by lapply:
Theoph$Dose[1] - NA
Theoph$Time[2] - NA
Theoph$conc[3] - NA
lapply(Theoph,is.na)
Is there a direct way to execute logical or across all vectors? The
following gives the
Alan S Barnett wrote:
How do I add to a trellis plot the best fit line from a robust fit? I
can use panel.lm to add a least squares fit, but there is no panel.rlm
function.
How about using panel.abline() instead of panel.lmline()?
fit1 - coef(lm(stack.loss ~ Air.Flow, data = stackloss))
try the following:
as.logical(rowSums(is.na(Theoph)))
## or
!complete.cases(Theoph)
I hope it helps.
Best,
Dimitris
Dimitris Rizopoulos
Ph.D. Student
Biostatistical Centre
School of Public Health
Catholic University of Leuven
Address: Kapucijnenvoer 35, Leuven, Belgium
Tel:
A Ezhil wrote:
Hi All,
I have a vector of length 48, something like:
0 0 1 1 0 0 0 0 0 0 0 0 0 0 1 1 1 1 1 1 1 1 1 1 1 1 1
1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1
I would like to print (reformat) this vector as:
00110011
by simply removing the
Hi All,
I have a vector of length 48, something like:
0 0 1 1 0 0 0 0 0 0 0 0 0 0 1 1 1 1 1 1 1 1 1 1 1 1 1
1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1
I would like to print (reformat) this vector as:
00110011
by simply removing the spaces between them. I
Tim Bergsma wrote:
Suppose I have a list of logicals, such as returned by lapply:
Theoph$Dose[1] - NA
Theoph$Time[2] - NA
Theoph$conc[3] - NA
lapply(Theoph,is.na)
Is there a direct way to execute logical or across all vectors? The
following gives the desired result, but seems
Suppose I have a list of logicals, such as returned by lapply:
Theoph$Dose[1] - NA
Theoph$Time[2] - NA
Theoph$conc[3] - NA
lapply(Theoph,is.na)
Is there a direct way to execute logical or across all vectors? The
following gives the desired result, but seems unnecessarily complex.
How do I add to a trellis plot the best fit line from a robust fit? I
can use panel.lm to add a least squares fit, but there is no panel.rlm
function.
--
Alan S Barnett [EMAIL PROTECTED]
NIMH/CBDB
__
R-help@stat.math.ethz.ch mailing list
a little simplier:
apply(do.call(rbind,lapply(Theoph,is.na)),2,any)
or
!complete.cases(Theoph)
On 6/8/07, Tim Bergsma [EMAIL PROTECTED] wrote:
Suppose I have a list of logicals, such as returned by lapply:
Theoph$Dose[1] - NA
Theoph$Time[2] - NA
Theoph$conc[3] - NA
lapply(Theoph,is.na)
Dear All,
discussing with a statistician of a pharmaceutical company I received
this answer about the statistical package that I have planned to use:
As R is not a validated software package, we would like to ask if it
would rather be possible for you to use SAS, SPSS or another approved
On Fri, 2007-06-08 at 06:13 -0700, A Ezhil wrote:
Hi All,
I have a vector of length 48, something like:
0 0 1 1 0 0 0 0 0 0 0 0 0 0 1 1 1 1 1 1 1 1 1 1 1 1 1
1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1
I would like to print (reformat) this vector as:
On the 'lsoda' help page, I did not see any option to force some
or all parameters to be nonnegative.
Have you considered replacing the parameters that must be
nonnegative with their logarithms? This effective moves the 0 lower
limit to (-Inf) and seems to have worked well for me
--- Jason Barnhart [EMAIL PROTECTED] wrote:
Hi John,
a) The NA appears because '30/02/1995' is not a
valid date.
strptime('30/02/1995' , %d/%m/%Y)
[1] NA
I knew we should never have moved to the Gregorian
Calender!
Thanks. I accidently made up the date but this means
that
Dear all,
Is there a way to plot / calculate pointwise confidence bands or
interval values for a non parametric regression like sm.regression?
Thank you in advance.
Regards,
Martin
__
R-help@stat.math.ethz.ch mailing list
People, don't get angry at the pharma statistician, he is just trying to
abide by an FDA requirement that is designed to insure that test perform
reliably the same. There is no point in getting into which product is
better. As far as the FDA rules are concerned a validated system beats a
better
Thanks all for the many excellent suggestions!
!complete.cases(Theoph) is probably the most succinct form for the
current problem, while the examples with 'any' seem readily adaptable to
similar situations.
Kind regards,
Tim.
Dimitris Rizopoulos wrote:
try the following:
On 7 June 2007 at 17:22, Tim Keitt wrote:
| I'm just curious if anyone else has had problems with this
| configuration. I added the CRAN repository to apt and installed 2.5.0
| with apt-get. I then did an install.packages(Rmpi) on cluster nodes.
| Rmpi loads and lamhosts() shows the nodes, but
I had mentioned exactly the same thing to others and the feedback I got is -
'when you have a hammer, everything will look like a nail'
^_^.
On 6/7/07, Frank E Harrell Jr [EMAIL PROTECTED] wrote:
Robert Wilkins wrote:
As noted on the R-project web site itself ( www.r-project.org -
Manuals -
Is there an example available of this sort of problematic data that
requires this kind of data screening and filtering? For many of us,
this issue would be nice to learn about, and deal with within R. If a
package could be created, that would be optimal for some of us. I
would like to
I like to know the answer as well.
To be honest, I really have hard time to understand the mentality of
clinical trial guys and rather believe it is something related to job
security.
On 6/8/07, Giovanni Parrinello [EMAIL PROTECTED] wrote:
Dear All,
discussing with a statistician of a
Looks much better. I seldom use dates for much and
didn't think to look at the sort.POSIXlt function.
If I understand this correctly the sort.POSIXlt with
na.last = FALSE is dropping all the NAs. Very nice.
--- Gabor Grothendieck [EMAIL PROTECTED]
wrote:
Perhaps you want one of these:
On 6/7/07, Alan S Barnett [EMAIL PROTECTED] wrote:
How do I add to a trellis plot the best fit line from a robust fit? I
can use panel.lm to add a least squares fit, but there is no panel.rlm
function.
It's not trellis, but it's really easy to do this with ggplot2:
install.packages(ggplot2,
On the R wiki site there is a general-purpose function
(sort.data.frame) that allows you to do this:
sort(df, by=~ x-z)
See: http://wiki.r-project.org/rwiki/doku.php?id=tips:data-frames:sort
Regards,
Kevin
On 6/8/07, Gunther Höning [EMAIL PROTECTED] wrote:
Dear list,
I have a very short
--- Peter Lercher [EMAIL PROTECTED] wrote:
Moving from S-plus to R I encountered many great
features and a much
more stable system.
Currently, I am left with 2 problems that are
handled differently:
1) I did lots of overplots in S-Plus using
par(new=T,xaxs='d',yaxs='d') to fix the axes
Giovanni Parrinello wrote:
Dear All,
discussing with a statistician of a pharmaceutical company I received
this answer about the statistical package that I have planned to use:
As R is not a validated software package, we would like to ask if it
would rather be possible for you to use
hello,
I have just a question before the week end it's that I don't know how to do to
paste matrixs and these matrix they have one same column and I'd like to paste
its by this column
and I wanna paste its not below but just at right side hand
thanks good week end
Hi ComRades,
I want to make a matrix of frequencies from vectors of a continuous
variable spanning different values. For example this code
x-c(runif(100,10,40),runif(100,43,55))
y-c(runif(100,7,35),runif(100,37,50))
z-c(runif(100,10,42),runif(100,45,52))
a-table(ceiling(x))
b-table(ceiling(y))
On 7 June 2007 at 14:27, [EMAIL PROTECTED] wrote:
| how can I use R in a pipline like this
|
| $ ./generate-data | R --script-file=Script.R | ./further-analyse-data
result.dat
The 'r' in our 'littler' package can do that. One example we show on the
littler webpage is
$ ls -l /boot |
The code in my post uses Date class, not POSIX.
sort.POSIXlt is never invoked. Suggest you read the
help desk article in R News 4/1 for more.
On 6/8/07, John Kane [EMAIL PROTECTED] wrote:
Looks much better. I seldom use dates for much and
didn't think to look at the sort.POSIXlt function.
If
On 6/7/07, Prof Brian Ripley [EMAIL PROTECTED] wrote:
f - function(x, dat) evalq(x, dat)
f(o, D)
Error in eval(expr, envir, enclos) : object o not found
g - function(x, dat) eval(x, dat)
g(o, D)
Error in eval(x, dat) : object o not found
What am I doing wrong? This seems to be
I'm not at all certain I understand your question, but try
?cbind
Sarah
On 6/8/07, elyakhlifi mustapha [EMAIL PROTECTED] wrote:
hello,
I have just a question before the week end it's that I don't know how to do
to paste matrixs and these matrix they have one same column and I'd like to
Dear Ruixin:
Among others, text mining is dealing with non-structural data while
data mining mainly focuses on structural one. Many algorithms can be
shared b/w them; however, some necessary data preprocessing is
required for text mining. There are a lot of online-resource there.
As to packages
agree with Frank.
as far as I've known, FDA doesn't encourage or discourage the usage of
software.
On 6/8/07, Frank E Harrell Jr [EMAIL PROTECTED] wrote:
Giovanni Parrinello wrote:
Dear All,
discussing with a statistician of a pharmaceutical company I received
this answer about the
On 6/8/07, Dirk Eddelbuettel [EMAIL PROTECTED] wrote:
On 7 June 2007 at 17:22, Tim Keitt wrote:
| I'm just curious if anyone else has had problems with this
| configuration. I added the CRAN repository to apt and installed 2.5.0
| with apt-get. I then did an install.packages(Rmpi) on cluster
Peter Lercher wrote:
-Original Message-
From: [EMAIL PROTECTED]
[mailto:[EMAIL PROTECTED] On Behalf Of Peter Lercher
Sent: Friday, June 08, 2007 3:07 AM
To: r-help@stat.math.ethz.ch
Subject: [R] overplots - fixing scientific vs normal notation
in output
Moving from S-plus to
Sicotte, Hugues Ph.D. wrote:
People, don't get angry at the pharma statistician, he is just trying to
abide by an FDA requirement that is designed to insure that test perform
reliably the same. There is no point in getting into which product is
better. As far as the FDA rules are concerned a
I don't think the code below does what's requested, as it assumes a single
overall fit for all panels, and I think the requester wanted separate fits
by panel. This can be easily done, of course, by a minor modification:
xyplot( y ~ x | z,
panel = function(x,y,...){
Frank et. al:
I believe this is a bit too facile. 21 CFR Part 11 does necessitate a
software validation **process** -- but this process does not require any
particular software. Rather, it requires that those using whatever software
demonstrate to the FDA's satisfaction that the software does
On 6/7/07, Alan S Barnett [EMAIL PROTECTED] wrote:
How do I add to a trellis plot the best fit line from a robust fit? I
can use panel.lm to add a least squares fit, but there is no panel.rlm
function.
Well, panel.lmline (not panel.lm, BTW) is defined as:
panel.lmline
function (x, y, ...)
{
I may have overstated things a bit.
See section VIII
http://www.fda.gov/CDER/GUIDANCE/2396dft.htm
If you are analyzing data your statistical package does not necessarely
have to be validated. You may have to show that the statistical methods
are adequate/appropriate or that the results are
As I read 21 CFR 11, the regulation deals more with ensuring the security
of the electronic health record itself. Thus, it seemed to me that so long
as the software (SAS, R, Splus, etc.) could not alter the data base in any
way then you're fine (this may be naive, but that's how I understood it).
Hello,
maybe my question ist stupid, but I would like to calculate a new
variable for all cases in my dataset. Inspired by the dialog in Rcmdr
I tried
Datenmatrix$cohigha- with(Datenmatrix,mean (c(M2ORG, M5ORG, M8ORG,
M11ORG), na.rm = TRUE)
as output I got the same number for all my cases
I am working on a plot and would be like to click on a few points and then
have a line connect them. Could anyone help me with this or advise me in a
direction that would suit this. I know I would be using ievent.wait in iplot
but not sure about this.
thank you.
--
View this message in context:
Bert Gunter wrote:
Frank et. al:
I believe this is a bit too facile. 21 CFR Part 11 does necessitate a
software validation **process** -- but this process does not require any
For database software and for medical devices -
particular software. Rather, it requires that those using whatever
Matthias von Rad wrote:
Hello,
maybe my question ist stupid, but I would like to calculate a new
variable for all cases in my dataset. Inspired by the dialog in Rcmdr
I tried
Datenmatrix$cohigha- with(Datenmatrix,mean (c(M2ORG, M5ORG, M8ORG,
M11ORG), na.rm = TRUE)
as output I got the
Bert,
I just want to make sure what I said is not overstated to offend
statistician who use SAS. actually, i am using SAS daily and able to
use it pretty well. ^_^
What I meant are:
1) I don't understand the mentality
2) using SAS instead of R might be related to job-security.
which is very
Hi to all,
maybe the last question was not clear enough.
I did not found any hints how to decide whether it should use lower.tail
or not.
As it is an extra R-feature ( written in
http://finzi.psych.upenn.edu/R/Rhelp02a/archive/66250.html )
I do not find anything about it in any statistical books
Hi,
I am a complete newbe to R, so the following problem will probably be
trivial for most of you guys: I get an error message every time I try
to run a R file directly from the DOS shell.
My R file (test.R) is intended to create a basic graph and has a very
simple code:
x-rep(1:10,1)
You need to basically use table on factors with fixed pre-specified
levels. For example:
x - c(runif(100,10,40), runif(100,43,55))
y - c(runif(100,7,35), runif(100,37,50))
z - c(runif(100,10,42), runif(100,45,52))
xx - ceiling(x); yy - ceiling(y); zz - ceiling(z)
mylevels - min(
Martin Henry H. Stevens sent the following at 08/06/2007 15:11:
Is there an example available of this sort of problematic data that
requires this kind of data screening and filtering? For many of us,
this issue would be nice to learn about, and deal with within R. If a
package could be
Hi Martin,
Do please, at least, read the documentation for the package you are using!:
?sm.options ## sub: display
## Example
with(iris, sm.regression(Sepal.Length, Sepal.Width, display=se))
Regards,
Mark Difford.
M. P. Papadatos wrote:
Dear all,
Is there a way to plot /
R isn't in your path. Either change your path to include it or place
Rcmd.bat from batchfiles anywhere in your existing path:
http://code.google.com/p/batchfiles/
and then:
Rcmd BATCH ...whatever...
On 6/8/07, Sébastien Bihorel [EMAIL PROTECTED] wrote:
Hi,
I am a complete newbe to
For windows users, EpiData Entry http://www.epidata.dk/ is an
excellent (free) tool for data entry and documentation.--Dale
On 6/8/07, Chris Evans [EMAIL PROTECTED] wrote:
Martin Henry H. Stevens sent the following at 08/06/2007 15:11:
Is there an example available of this sort of
At 01:31 PM 6/8/2007, Carmen wrote:
Hi to all,
maybe the last question was not clear enough.
I did not found any hints how to decide whether it should use lower.tail
or not.
As it is an extra R-feature ( written in
http://finzi.psych.upenn.edu/R/Rhelp02a/archive/66250.html )
I do not find anything
Alternatively, use the full path in your call to R as I do below:
F:\Program Files\R\R-2.4.1pat\bin\R.exe CMD BATCH --vanilla --slave whatever.R
HTH,
Roger
-Original Message-
From: [EMAIL PROTECTED] [mailto:[EMAIL PROTECTED] On Behalf Of Gabor
Grothendieck
Sent: Friday, June 08,
I need to be able to run a generalized linear model with a log() link
and a Weibull family, or something similar to deal with an extreme
value distribution.
I actually have a large dataset where this is apparently necessary.
It has to do with recovery of forensic samples from surfaces, where
On 6/8/2007 11:33 AM, Zack Weinberg wrote:
On 6/7/07, Prof Brian Ripley [EMAIL PROTECTED] wrote:
f - function(x, dat) evalq(x, dat)
f(o, D)
Error in eval(expr, envir, enclos) : object o not found
g - function(x, dat) eval(x, dat)
g(o, D)
Error in eval(x, dat) : object o not found
Not to mention all the work that goes into PROC TEMPLATE and ANNOTATE to
make SAS graphs presentable! I suspect that a lot of companies don't use
SAS graphs or tables at all - they just export the data from SAS to Excel.
-Cody
Cody Hamilton, PhD
Edwards Lifesciences
What I would love to have
This is an expanded version of the question I tried to ask last night
- I thought I had it this morning, but it's still not working and I
just do not understand what is going wrong.
What I am trying to do is write a wrapper for lattice xyplot() that
passes a whole bunch of its secondary
On 6/8/07, Duncan Murdoch [EMAIL PROTECTED] wrote:
After a bit more experimentation I figured out that this does what I want:
h - function(x, d) eval(substitute(x), d, parent.frame())
but I don't understand why the substitute() helps, or indeed why it
has any effect at all...
Within
Hi,
Does anyone know how to count the number of modes in 2 dimensions using
kde2d function?
Thanks
Pat
__
R-help@stat.math.ethz.ch mailing list
https://stat.ethz.ch/mailman/listinfo/r-help
PLEASE do read the posting guide
On Fri, 2007-06-08 at 16:02 +0200, Giovanni Parrinello wrote:
Dear All,
discussing with a statistician of a pharmaceutical company I received
this answer about the statistical package that I have planned to use:
As R is not a validated software package, we would like to ask if it
would
On Fri, 8 Jun 2007, Robert A. LaBudde wrote:
I need to be able to run a generalized linear model with a log() link
and a Weibull family, or something similar to deal with an extreme
value distribution.
The Weibull with log link is not a GLM, but survreg() in package survival
can fit it, as
As you may not be surprised to hear, no sooner did I post the previous
message than I realized I had a really dumb mistake. I've now gotten
a bit farther but am still stuck. New code:
graph - function (x, data, groups, xlab) {
pg - function(x, y, group.number, ...) fnord
body(pg) -
Dale Steele wrote:
For windows users, EpiData Entry http://www.epidata.dk/ is an
excellent (free) tool for data entry and documentation.--Dale
Note that EpiData seems to work well under linux using wine.
Frank
__
R-help@stat.math.ethz.ch mailing
The fact that FDA statisticians are using R also assuages one of the main
concerns that I have heard voiced about using R for FDA submissions - that
there would be no statisticians available at FDA to review R code which
would seriously delay the review of a submission.
Mark also brings up a
Note that the number of modes (local maxima??) is a function of the
bandwidth, so I'm not sure your question is even meaningful.
Bert Gunter
Genentech Nonclinical Statistics
South San Francisco, CA 94404
650-467-7374
-Original Message-
From: [EMAIL PROTECTED]
[mailto:[EMAIL PROTECTED]
Thanks for the reply,
maybe I shall say bumps, I can use persp to show a density on a X Y
dimensions.
one peak is one mode I think. I try to find an automatic way to detect how
many peaks of the densities.
Pat
Note that the number of modes (local maxima??) is a function of the
bandwidth, so
1 - 100 of 112 matches
Mail list logo