Benilton Carvalho schrieb:
Well, AFAIK, the definition of a p-value is the probability of
observing something at least as extreme as the observed data.
If you observed z, and Z follows a std-normal
p-value = P( Z -abs(z) ) + P( Z abs(z) )
= 2*P ( Z abs(z) )
= 2*pnorm(z,
MHHS == Martin Henry H Stevens [EMAIL PROTECTED]
on Mon, 11 Jun 2007 17:23:46 -0400 writes:
MHHS Hi Folks, I was wondering what everyone thought about
MHHS adding a sentence to each package update announcement
MHHS that described what the package did. R extensions are
MHHS
Dear all R users
:
I am a IT student with few statistical background and new R user for only
have two month exprience. I have a data named medcost, import by
read.table() as follow for example (real dataset has 500 cases), the
heander id means case id, member means members in a family and
Tim Holland wrote:
Is there a way in R to select certain characters from a line of text? I
have some data that is presently in a large number of text files, and I
would like to be able to select elements of each text file (elements are
always on the same line, in the same position) and
-Original Message-
From: [EMAIL PROTECTED] [mailto:[EMAIL PROTECTED]
On Behalf Of Mary Royerr
Sent: Monday, June 11, 2007 8:14 PM
To: r-help@stat.math.ethz.ch
Subject: Re: [R] Textpad help
Thx for the response John and Jim.
I had heard before that textpad is a good editor for R
If you want to test whether ' member' has an effect on 'cost' (a continuous
numerical variable), I do not recommend using a chi.square test,
but rather a simple linear regression or a one-way analysis of variance.
Chi.square are for categorical variables and unless you have a good reason,
there
spime sabya23 at gmail.com writes:
Hi all,
I need some solution in the following problem. The following error appears
when i use mgcv package for implementing GAM. But the same formula works
fine in gam package.
model.gam - gam(formula = RES ~
Dear Charlie
dat - data.frame(id = 1:10, member = c(4,2,3,5,6,2,4,6,3,4),
cost = c(320,150,420,330,540,310,169,647,347,567))
dat[,costF] - cut(dat[,cost], breaks = seq(100, 700, by=100))
table(dat[,costF], dat[,member])
This should create the table you like.
Best regards,
Is there a standardized way in R to perform post-hoc comparisons for main
adn interaction effects of between- and within-subject factors?
For instance, I have a data set of performance of adults of different age
groups (20-30, 60-70,70-80) performing a WM task (n-back, with n=1,2,3,4) in
two
G'day Hadley,
On Wed, 30 May 2007 10:57:54 +0200
hadley wickham [EMAIL PROTECTED] wrote:
We're pleased to announce a one day course covering static and dynamic
graphics using R, ggplot and GGobi. The course will be held just
before the JSM, on Saturday, 28 July 2007, in Salt Lake City. The
Dear Mr. Oksanen,
First of all thanks for your reply. I have solved this problem in this way.
My data consists of some categorical(CAT..) predictors and also some
numerical variables(NUM..) have only {0,1} 0r {0,1,2,3} values. For applying
GAM i just didnot consider their splines. I had came to
Hello,
I work under windows with two versions of R : 2.4.1 and 2.2.1
Ghostscript tool is also installed (gswin32c Version 8.54).
I have a small sample program:
bitmap(a.bmp)
m=(1)
plot(m)
dev.off()
With R 2.2.1, it works well, but with R 2.4.1 it fails when gswin32c is launched
with
Hi everyone,
I want to make use of the save function but it did not work.
I'm using vista and R 2.5.0, winzip is installed too.
Here's the code (from example ?save):
x - runif(20)
y - list(a = 1, b = TRUE, c = oops)
save(x, y, file = xy.Rdata)
Fehler in gzfile(file, wb) : kann
Thanks again. Both layouts look very usefull and
certainly a lot better than I was getting on my own.
--- Deepayan Sarkar [EMAIL PROTECTED] wrote:
On 6/11/07, John Kane [EMAIL PROTECTED] wrote:
--- [EMAIL PROTECTED] wrote:
On 6/9/07, John Kane [EMAIL PROTECTED] wrote:
Is it
I would like to compute a direct hedonic prize-index, calld
Time-Dummy-Variable-Method. Is this method implmented in R? and further
indirect methods?
Best regards
[[alternative HTML version deleted]]
__
R-help@stat.math.ethz.ch mailing list
So you are running R somewhere in your file system where you do not have
permission to write.
Did you create a shortcut with a working directory set as the rw-FAQ
advised you to?
On Tue, 12 Jun 2007, Maja Schröter wrote:
Hi everyone,
I want to make use of the save function but it did not
Hi Martin,
Maybe we should appeal to each others' vanity (Just think how many
hundreds---even thousands---of people will use your package once they
know about it) Hmmm. Maybe that would have the opposite of the
intended effect.
Cheers,
Hank
On Jun 12, 2007, at 3:16 AM, Martin Maechler
The following gives two functions for producing distribution graphs:
distribution-graph
produces a single graph, and
multiple.distribution.graph
produces a number of graphs side by side.
Regards,
Tore Wentzel-Larsen
statistician
Centre for Clinical research
Armauer Hansen house
Haukeland
Hello,
I work under windows xp and I am trying to build a R package with a subroutine
written in fortran 90.
I have installed all the updated tools and I am working with R-2.4.0 or
R-2.5.0.
When I check a package with a subroutine in fortran 77 (and extension f)
everything is ok.
When I
Dear Xiaohui,
Thanks a lot for your help! It worked for me.
Best regards!
JP-
--On Monday, June 11, 2007 9:25 AM -0700 Xiaohui [EMAIL PROTECTED]
wrote:
see ?try
Jianping Jin wrote:
Dear group:
I wrote a code to iterate a non-linear fit with a set of data. The
entire process didn't
The tools do not include an F9x compiler: you need to edit MkRules to
specify the path to one (after installing it) (and you did last March).
If I do that, I can install and check packages such as oc and wnominate
from CRAN that make use of F9x.
On Tue, 12 Jun 2007, Cinzia Viroli wrote:
Hello,
I have a question about the lars package. I am using this package to get the
coefficients at a specific LASSO parameter s.
data(diabetes)
attach(diabetes)
object - lars(x,y,type=lasso)
cvres-cv.lars(x,y,K=10,fraction = seq(from = 0, to = 1, length = 100))
fits - predict.lars(object,
1) I have a data that includes some bad data. I want to make a
trellis plot where each panel contains
a) A scatter plot of the good data
b) A scatter plot of the bad data in a different color
c) A best fit line of all the data, and
d) A best fit line of the good data.
I tried using xyplot and
I apologize for the off-topic post, but my Google search did not turn
up much and I thought people on this list my have knowledge of this.
I am looking for examples of data brushing (i.e. dynmaic linked
plots) either on a web site, or in a web-based application, such as
an AJAX app.
I am having trouble figuring out the right form for the nlme arguments. I do
have examples in Modern and Applied Statistics with S and from other sources,
but I still can't figure it out.
I am trying to estimate species richness (sr) in streams across minnesota. My
predictor variables are
On 6/12/07, Roy Mendelssohn [EMAIL PROTECTED] wrote:
I apologize for the off-topic post, but my Google search did not turn
up much and I thought people on this list my have knowledge of this.
I am looking for examples of data brushing (i.e. dynmaic linked
plots) either on a web site, or in a
I need to process some datasets where the configuration information was
stored in .INI-like files, i.e., text files with sections like this:
[Section1]
var1=value1
var2=value2
[Section2]
A=value3
B=value4
...
From Google and other searches I haven't found any package, or function
within a
I would like to fit a Pareto Distribution and I am using the following codes.
First, I thought the fitted (fit1) should be the fitted value for the data,
is it correct? As the result of the fitted turns out to be the same value.
fit=vglm(ycf1 ~ 1, pareto1(location=alpha), trace=TRUE, crit=c)
Earl:
.Ini files are, for lack of a better description, ancient.
There are old windows functions such as GetProfileString.
However you will have to make reference to load these from the windows
Kernel.dll.
Probably not worth the effort to code really old things as .ini files.
From what I see of
Hi,
This is probably trivial to most people out there, but I'm struggling with this.
I have a data set which contains the closing prices (properly adjusted for
dividends and splits) for several hundred securities and the closing prices for
a general stock market index (SP 500). I have no
Hi,
I would like to get hold of the R version 2.4.0 windows installer, it
doesn't seem to be available on CRAN (except the source, which needs
compiling). Does anyone know if it's still available anywhere?
Thanks,
jeremy
--
Jeremy Miles
Learning statistics blog:
I am trying to access a dll with dyn.load, but I get an error message box
titled R Console: Rgui.exe - Unable to Locate Component. The error message
itself states this application has failed to start because libifcoremdd.dll
was not found. Re-installing the application may fix this problem. I
Hi all,
I have written a script in R that simulates genetically informative
data - it is posted on my website and available to the public. This is
my first time to write a script for use by others and am learning that
it isn't as easy as it seems.
To the issue. My script runs fine on my machine
Dear users,
In my psychometric test i have applied logistic regression on my data. My
data consists of 50 predictors (22 continuous and 28 categorical) plus a
binary response.
Using glm(), stepAIC() i didn't get satisfactory result as misclassification
rate is too high. I think categorical
Here is some code. It replaces [ and ] with = sign and reads the result
into a data frame, DF. DF2 is similar except the section is now in V3.
DF3 is like like DF2 except sections are carried forward and finally
we remove the rows which only had sections.
Lines.raw - [Section1]
var1=value1
R 2.5 under windows XP
R 2.5 under Linux FC 6
I currently use R under windows and when I do so, I get to work in a nice
windows environment complete with a very useful menu bar that includes a FILE,
EDIT, MISC, PACKAGE, WINDOWS, and HELP menus. I have begun to use R under
LINUX. When I run
Hi,
First of all, I apologize for sending out this message, as I'm sure
the answer is on the archives, but I just can't find it.
Not long ago, there was a discussion about building Windows packages
from the source code and someone posted a link to a website to which
we could submit the
Dear list,
First apologize that this is trivial and just betrays my slothfulness at the
keyboard. I'm sick of having to type a long name just to get a glimpse of
something. For example, if my data frame is named
'AuroraStochasticRunsJune1.df and I want to see what the middle looks like, I
Jeremy Miles wrote:
Hi,
I would like to get hold of the R version 2.4.0 windows installer, it
doesn't seem to be available on CRAN (except the source, which needs
compiling). Does anyone know if it's still available anywhere?
Thanks,
jeremy
Hi Jeremy:
The windows binary for 2.4.0
Dear Latticer,
I want to give individual colors to all elements in a simple stacked
barchart. I know why the example below does not work (and it is a excellent
default), but is there any workaround for this?
Dieter
# This only colors red and green, but I want blue and gray for Peatland.
At 01:45 PM 6/12/2007, Tirtha wrote:
Dear users,
In my psychometric test i have applied logistic regression on my data. My
data consists of 50 predictors (22 continuous and 28 categorical) plus a
binary response.
Using glm(), stepAIC() i didn't get satisfactory result as misclassification
rate is
There is, two candidates are ESS-emacs and JGR as a gui. Have a look at:
http://www.r-project.org/GUI
, on Tuesday 12 June 2007 19:42:33 John Sorkin wrote:
JS R 2.5 under windows XP
JS R 2.5 under Linux FC 6
JS
JS I currently use R under windows and when I do so, I get to work in a
nice
Ian McCarthy wrote:
I am trying to access a dll with dyn.load, but I get an error message box
titled R Console: Rgui.exe - Unable to Locate Component. The error
message itself states this application has failed to start because
libifcoremdd.dll was not found. Re-installing the application
On 6/10/07, R.H. Koning [EMAIL PROTECTED] wrote:
Hello, I have some problems in using chron, Hmisc, and lattice. First,
using both chron and Hmisc, I get an error message when describing data:
df$Date - chron(df$Date,format=c(d/m/y))
ll - latex(describe(df),file=..//text//df.tex)
Error in
Hy,
no need to find the best s value. CV does it for you:
cvres-cv.lars(X,Y,K=10,type='lasso')
sAtBest-cvres$fraction[which.min(cvres$cv)]
fits - predict.lars(object, type=coefficients, s=sAtBest, mode=fraction)
...
Ciao
Bruno
Hello,
I have a question about the lars package. I am using this
On 6/12/07, Dieter Menne [EMAIL PROTECTED] wrote:
Dear Latticer,
I want to give individual colors to all elements in a simple stacked
barchart. I know why the example below does not work (and it is a excellent
default), but is there any workaround for this?
Dieter
# This only colors red
On 6/12/07, Dieter Menne [EMAIL PROTECTED] wrote:
Dear Latticer,
I want to give individual colors to all elements in a simple stacked
barchart. I know why the example below does not work (and it is a excellent
default), but is there any workaround for this?
Dieter
# This only colors red
On Tue, 12 Jun 2007, Benilton Carvalho wrote:
Hi,
First of all, I apologize for sending out this message, as I'm sure
the answer is on the archives, but I just can't find it.
Not long ago, there was a discussion about building Windows packages
from the source code and someone posted a link
Hi Patrick,
Here is a simple R code for locating ALL the local maximum of a bivariate
function, which is known on a rectangular grid. I have illustrated it with
a function called the Branin function, which is commonly used as a test
function in the global optimization literature. It has 6 local
hadley wickham h.wickham at gmail.com writes:
On 6/12/07, Dieter Menne dieter.menne at menne-biomed.de wrote:
Dear Latticer,
I want to give individual colors to all elements in a simple stacked
barchart. I know why the example below does not work (and it is a excellent
default), but
Patrick,
Here is an example closer to what you are interested in - locating bumps in
kernel density estimator. I am using the example from package KernSmooth,
using the function bkde2D().
# Another example for locating maxima in kernel density estimation
data(geyser, package=MASS)
x -
Thanks for the information. These commands work perfectly fine and the
?panel.superpose help was most informative.
If you don't mind, I will certainly come back to you as there will
certainly be additionnal complexities in my datasets that I will be
unable to handle (e.g. data in more than one
I need to create artificial datasets with specific correlation
coefficients (i.e. a dataset that returns r = 0.30, etc.) as examples
for a lab I am teaching this summer. Is there a way to do that in R?
Thanks.
Jim Milks
Graduate Student
Environmental Sciences Ph.D. Program
136 Biological
On 6/12/07, Seb [EMAIL PROTECTED] wrote:
Hello
I apologize in advance if this question has already be posted on the
list, although I could not find a relevant thread in the archives.
I would like to overlay xyplots using different datasets for each plot.
I typically work on the following
see mvrnorm in MASS and especially the empirical argument
James Milks james.milks at wright.edu writes:
I need to create artificial datasets with specific correlation
coefficients (i.e. a dataset that returns r = 0.30, etc.) as examples
for a lab I am teaching this summer. Is there a
Dear Listers:
I have a couple of data frames to report and each corresponds to
different condtions, e.g. conditions=c(10, 15, 20, 25). In this
examples, four data frames need to be exported in a pretty report.
I knew Perl has some module for exporting data to Excel and after
googling, I found R
Hello,
It looks to me as if the ! logical operator cannot be called when
subsetting a data.frame. In the example below, the value column has two
factor levels (but my typical datasets have more), and what I am trying
to do is to exclude all lines for which the value is different from
A. I
Just realized append=T might help for xtable, but I think it still not
very pretty :)
On 6/12/07, Weiwei Shi [EMAIL PROTECTED] wrote:
Dear Listers:
I have a couple of data frames to report and each corresponds to
different condtions, e.g. conditions=c(10, 15, 20, 25). In this
examples, four
Hello everyone,
Hope that my question could have interest for more than one people here.
I know that in order to compare mixed models with different fixed effect
I need to use ML method.
But ... what about comparing models with the same fixed effects and
different random effects ?
Thanks in
Sébastien pomchip at free.fr writes:
Can you please point to me my syntax mistake or indicate a method to get
this type of data.frame subset ?
Thank you in advance
ID value
1 1 1.2
2 2 1.2
3 3 1.2
4 4 1.2
5 5 A
6 6 A
7 7 A
8 8 A
I am new to using R and would appreciate some advice on
which books to start with to get up to speed on using R.
My Background:
1-C# programmer.
2-Programmed directly using IMSL (Now Visual Numerics).
3- Used in past SPSS and Statistica.
I put together a list but would like to pick the best of
Hi,
I have all CRAN packages installed on my Linux. Now I have problems with JGR.
When I make a plot and close a device, the device dont work anymore, I nedd to
use before javaGD() and after plot(). When I try do close JGR and save a
session, It return a erro and dont close. Look:
Exception
I'm trying to proess date/time fields from files that were given to me to
analyze.
Any clues what I'm doing wrong with strptime? This seems to fail the same
way under Linux or Windows.
For ?strptime would it make sense to explain %OS3 somewhere besides the
Examples?
# Why does %OS3 work
Generating Excel or reports seems to be two different questions:
1. Reports. You can use xtable together with Sweave. See figure 1 in this
link for an example:
http://www.ci.tuwien.ac.at/~leisch/Sweave/Sweave-Rnews-2002-3.pdf
2. Regarding writing Excel files this can be done on Windows using
[EMAIL PROTECTED] a écrit :
I am new to using R and would appreciate some advice on
which books to start with to get up to speed on using R.
My Background:
1-C# programmer.
2-Programmed directly using IMSL (Now Visual Numerics).
3- Used in past SPSS and Statistica.
I put together a list
On 6/12/07, Alan S Barnett [EMAIL PROTECTED] wrote:
1) I have a data that includes some bad data. I want to make a
trellis plot where each panel contains
a) A scatter plot of the good data
b) A scatter plot of the bad data in a different color
c) A best fit line of all the data, and
d) A
Alain,
Can you tell us what you plan to use R for?
Regards,
-Cody
[EMAIL PROTECTED] a écrit :
I am new to using R and would appreciate some advice on
which books to start with to get up to speed on using R.
My Background:
1-C# programmer.
2-Programmed directly using IMSL (Now Visual
hi
i just learn how to use R in my research
i used extreame value package to get the return level of rainfall.
the output i got as the graph plot between return level and period.
i wonder how i can get the value from the graph as a table format.
please help me
thanks
jam
[[alternative
On 12 June 2007 at 16:54, Earl F. Glynn wrote:
| Any clues what I'm doing wrong with strptime? This seems to fail the same
| way under Linux or Windows.
|
| For ?strptime would it make sense to explain %OS3 somewhere besides the
| Examples?
|
| # Why does %OS3 work here?
|
Hi,
I have created a dll in Fortran and used the Visual Studio 2005 Compiler. I've
read that certain problems can arise based on the compiler used and that these
problems can sometimes cause R not to be able to access everything
appropriately, but I've not found how to fix any of these
In thinking about this a bit more here is an even shorter solution where
Lines.raw is as before:
# Lines - readLines(myfile.ini)
Lines - readLines(textConnection(Lines.raw))
Lines2 - chartr([], ==, Lines)
DF - read.table(textConnection(Lines2), as.is = TRUE, sep = =, fill = TRUE)
L - DF$V1 ==
On 12-Jun-07 17:45:44, Tirthadeep wrote:
Dear users,
In my psychometric test i have applied logistic regression
on my data. My data consists of 50 predictors (22 continuous
and 28 categorical) plus a binary response.
Using glm(), stepAIC() i didn't get satisfactory result as
On 12-Jun-07 20:54:05, Ken Knoblauch wrote:
see mvrnorm in MASS and especially the empirical argument
James Milks james.milks at wright.edu writes:
I need to create artificial datasets with specific correlation
coefficients (i.e. a dataset that returns r = 0.30, etc.) as examples
for
Hi
I'm trying to cluster a continuous dataset with a varying number of clusters
and with a restriction that each cluster must have more than 'x' number of
observations.
I have tried the clara function, using silhouette to give me the neighbouring
cluster mediod of each observation, then
I don't understand, what's wrong with using prcomp in this situation?
On Sunday 10 June 2007 12:50 pm, Ranga Chandra Gudivada wrote:
Hi,
I was wondering whether there is any package implementing Principal
Component Analysis for Binary data
At 5:01 PM -0400 6/12/07, Weiwei Shi wrote:
Dear Listers:
I have a couple of data frames to report and each corresponds to
different condtions, e.g. conditions=c(10, 15, 20, 25). In this
examples, four data frames need to be exported in a pretty report.
I knew Perl has some module for exporting
In clinical trial data preparation and many other data situations, the
statistical programmer needs to merge and re-merge multiple input
files countless times. A syntax for merging files that is clear and
concise is very important for the statistical programmer's
productivity.
Here is how Vilno
Hi,I know only mle function but it seems that in mle one can only specify the
bound of the unknowns forming the likelihood function. But I would like to
specify something like, a = 2b or a = 2b where 'a' and 'b' could be my
parameters in the likelihood function. Any help would be really
Here is yet another simplification. This one uses na.locf from the zoo
package to shorten it further and also make it easier to understand.
Below we have one line to read in the .ini file, one line to transform the
characters [ and ] to = and =, the read.table line parses the result and
the next
At 09:13 PM 6/12/2007, Don wrote:
At 5:01 PM -0400 6/12/07, Weiwei Shi wrote:
Dear Listers:
I have a couple of data frames to report and each corresponds to
different condtions, e.g. conditions=c(10, 15, 20, 25). In this
examples, four data frames need to be exported in a pretty report.
I
Thanks to all of the people who responded. What I was trying to do is to
turn my matrix or frame containing index level returns and stock returns
into a matrix of betas. I don't really need to worry about risk-free
interest rates. I just need to be able to come up with a number that shows
the
The problem with applying prcomp to binary data is that it's not
clear what problem you are solving.
The standard principal components and factor analysis models
assume that the observations are linear combinations of unobserved
common factors (shared variability), normally
Jam,
If you use fpot the GPD of package evd to calculate your
return periods, you can generate the table you want by
using the following commands:
RP - c(10,100,500,1000) #return periods
rf - your clean vector of rainfall data
thres - the threshold for your data
nopy - number of observations
I'm confused:
Have you considered the 'examples' in the 'polr' help file? The
first example ends summary(house.plr). The print of this summary
includes standard errors. If you want those numbers for subsequent
computations, you can try str(summary(house.plr)) or
You might try (detrended) correspondence analysis, which is designed for
count data, if it makes sense to treat your binary data that way.
I've used ade4 and also vegan, and they are both good packages for these
types of ordinations. You could also look at non-metric multidimensional
scaling.
Dear Listers
I have tried to install Frank Harrell's two libaries: Hmisc
and Design.
I found that Hmisc was listed in the list of packages from
the Install Packages command on the Packages menu, but
Design was not. I installed Hmisc from this list, and when I
issued the library(Hmisc)
On 6/13/07, Robert Wilkins [EMAIL PROTECTED] wrote:
The point is : there are lots of data preparation scenarios where
large numbers of merges need to be done. This is an example where
Vilno and SAS are easier to use than the competition. I'm sure an Awk
programmer can come up with something,
87 matches
Mail list logo