On Apr 16, 2007, at 1:37 AM, Dong-hyun Oh wrote:
Dear expeRts,
I would simplify following code.
-
youtput - function(x1, x2){
n - length(x1)
y - vector(mode=numeric, length=n)
for(i in 1:n){
if(x1[i] =5 x1[i] = 10 x2[i] =5 x2[i]
Hi, this is Charlie and I am trying to embed R in my apache server. However, I
am having problem with installation. The R projet says that we can install
add-on package with command $R CMD (rapache.0.1.4.tar.gz). However, I try to
use it in command windows (windows xp) and an error shows.
Yuchen == Yuchen Luo [EMAIL PROTECTED]
on Sun, 15 Apr 2007 12:18:23 -0700 writes:
Yuchen Dear Friends.
Yuchen I tried to use nls.control() to change the 'minFactor' in nls( ),
but it
Yuchen does not seem to work.
yes, you do not seem to use it correctly.
No reason to jump to
Charlie wrote:
Hi, this is Charlie and I am trying to embed R in my apache server. However,
I am having problem with installation. The R projet says that we can install
add-on package with command $R CMD (rapache.0.1.4.tar.gz). However, I try to
use it in command windows (windows xp) and
bonjour,
je vous écris pour savoir si vous aviez trouver réponse pour l'écriture du
test de Newman Keuls sous R car j'en ai aussi besoin
merci encore
cordialement.
-
[[alternative HTML version deleted]]
Mandy Barron wrote:
Hello
I am having problems indexing a subset dataframe, which was created
as:
waspsNoGV-subset(wasps,site!=GV)
Fitting a linear model revealed some data points which had high
leverage, so I attempted to redo the regression without these data
points:
Which `gam' are you using? If it's mgcv::gam then all the smooth terms are
represented using a basis expansion, so it is possible to write down the
functional form used. The forms used are given in Chapter 4 of my book
(see ?gam for reference). It is also possible for `s' terms (but not `lo'
On Mon, Apr 16, 2007 at 09:03:27AM +0200, Martin Maechler wrote:
Yuchen == Yuchen Luo [EMAIL PROTECTED]
on Sun, 15 Apr 2007 12:18:23 -0700 writes:
Yuchen Dear Friends.
Yuchen I tried to use nls.control() to change the 'minFactor' in nls( ),
but it
Yuchen does not seem to
Mandy Barron wrote:
Hello
I am having problems indexing a subset dataframe, which was created
as:
waspsNoGV-subset(wasps,site!=GV)
Fitting a linear model revealed some data points which had high
leverage, so I attempted to redo the regression without these data
points:
On Sunday 15 April 2007 23:36, Jin Huang wrote:
Dear all,
I fitted a non-parametric model using GAM function in R. i.e.,
gam(y~s(x1)+s(x2)) #where s() is the smooth function
Then I obtained the coefficients(a and b) for the non-parametric terms.
i.e., y=a*s(x1)+b*s(x2)
-- do you mean
Dear R-users,
I have some trouble to perform a 2-samples KS test.
Apparently my 2 samples are numerical (see below) but R
complains that the y term is not
-
str(subset(mydata,Identified==NO,select=KD))
Joerg == Joerg van den Hoff [EMAIL PROTECTED]
on Mon, 16 Apr 2007 11:15:23 +0200 writes:
Joerg On Mon, Apr 16, 2007 at 09:03:27AM +0200, Martin Maechler wrote:
Yuchen == Yuchen Luo [EMAIL PROTECTED]
[]
[]
Joerg a final remark (off-topic?),
Greetings,
In the following one-way ANOVA I am attempting to calculate the means of
each treatment along with their 95% Tukey confidence intervals for the data
shown below using a routine from the HH package.
library(HH)
options(digits=10)
# load data
treat
voltage
1
130
1
Dear Partecipants to the list,
I am Enrico Foscolo, a student of the Faculty of Statistics, University of
Bologna,
and I am interested under consideration of the search of zeros of one function.
I would want to write a code using the software R.
I have already read the book Computing the Zeros of
On 4/16/2007 1:37 AM, cottrell wrote:
I'm looking for a way to 'reuse' existing rgl device windows. Right
now, every time I run my script I have to close the preexisting windows
and the new windows get assigned ever-increasing numbers.
I know how to do it for regular R plotting device
Dear Partecipants to the list,
I am Enrico Foscolo, a student of the Faculty of Statistics, University of
Bologna,
and I am interested under consideration of the search of zeros of one function.
I would want to write a code using the software R.
I have already read the book Computing the Zeros of
On Mon, 16 Apr 2007, Stéphane CRUVEILLER wrote:
Dear R-users,
I have some trouble to perform a 2-samples KS test.
Apparently my 2 samples are numerical (see below) but R
That says there are data frames, not numeric vectors. That would be easier
to see if you printed the value (rather than
There is no need to send this twice under different subject lines.
On Mon, 16 Apr 2007, enrico.foscolo2 wrote:
Dear Partecipants to the list,
I am Enrico Foscolo, a student of the Faculty of Statistics, University of
Bologna,
and I am interested under consideration of the search of zeros
Hello,
I just read this paper, and was surprised by:
This study performs a comparison of the latest versions of nine
different statistical software packages using StRD. These packages are
SAS (9.1), SPSS (12.0), Excel (2003), Minitab (14.0), Stata (8.1), Splus
(6.2), R (1.9.1), JMP (5.0), and
Dear Partecipants to the list,
I am Enrico Foscolo, a student of the Faculty of Statistics,
University of Bologna,
and I am interested under consideration of the search of zeros of
one function.
I would want to write a code using the software R.
I have already read the book Computing
Hi,
This is a very basic question, but apparently I am too stupid for it.
I have a large matrix A, and I need to avoid for loops. How could I
apply a function f(a,r,c) on each element of A, using the subscript (row
and column) of a as the other arguments?
Thanks in advance,
Markku Karhunen
Hello,
I have a problem with readLines.
I have a data file with many informations added with a different number of
measurments (example at the end).
I only want to read the measurments witch start with START OF HEIGHT DATA and
end with END OF HEIGHT DATA.
The difficulty is:
-I want to read
It would be helpful if you could be more specific of what exactly
you'd like to compute. Have a look also at the posting guide available
at:
http://www.R-project.org/posting-guide.html
Best,
Dimitris
Dimitris Rizopoulos
Ph.D. Student
Biostatistical Centre
School of Public Health
See library(multcomp) and
?glht
?contrMat
for several procedures for multiple comparisons.
The Newman Keuls test is not on the list. The related Tukey method
is on the list.
__
R-help@stat.math.ethz.ch mailing list
Hi Charlie,
Charlie wrote:
Hi, this is Charlie and I am trying to embed R in my apache server. However,
I am having problem with installation. The R projet says that we can install
add-on package with command $R CMD (rapache.0.1.4.tar.gz). However, I try to
use it in command windows
Dear R-helpers Friedrich Leisch,
I'm trying to fit a mixture of regression models on attached data set
by doing the following:
agl - read.table(agl.txt)
mod1 - flexmix(resp~gng+csc|Subject,data=agl,
model=FLXglm(family=binomial),k=2)
The result is a (varying) number of the following
Hi,
is there a platform independent way to log a complete R session (input +
output + warnings + errors + ???) into a text file? I can use sink(file,
split=T) and savehistory() but it is quite cumbersome to merge the files
and still there is some loss (eg. warnings).
I wondered that I only
The behavior you are seeing is part of the design of the
summary.glht function in the multcomp package.
There is some randomness in the adjustment for multiplicity.
See
library(multcomp)
?summary.glht
?adjusted
__
R-help@stat.math.ethz.ch mailing list
ESS is the answer for Windows.
NT emacs is available at
ftp://ftp.gnu.org/pub/gnu/emacs/windows
download emacs-21.3-fullbin-i386.tar.gz
ESS is available at
http://ess.r-project.org
Click source and then right-click ess-5.3.3.zip
__
--- Charilaos Skiadas [EMAIL PROTECTED] wrote:
On Apr 14, 2007, at 7:06 AM, John Kane wrote:
--- Adrian Dusa [EMAIL PROTECTED] wrote:
Charilaos Skiadas skiadas at hanover.edu
writes:
[...]
I save as csv format all the time, and it offers
me a choice to use
the labels instead
Will this work for you?
x.in - EXAM NUM:2
+ -
+ EXAM #1
+ ASTIG:-2.4D
+ AXIS:4.8
+ START OF HEIGHT DATA
+ S= 0 y=0.0 x=0.
+ S= 0 y=0.1 x=0.00055643
+ S= 9 y=4.9 x=1.67278117
+ S= 9 y=5.0 x=1.74873257
+ S=10 y=0.0 x=0.
+ S=10 y=0.1 x=0.00075557
+ S=99 y=5.3
[EMAIL PROTECTED] napsal dne 16.04.2007 14:52:55:
Hi,
This is a very basic question, but apparently I am too stupid for it.
I have a large matrix A, and I need to avoid for loops. How could I
apply a function f(a,r,c) on each element of A, using the subscript (row
and column) of a as
Ingmar,
I assume that resp is a vector of zeros and ones.
mod1 - flexmix(resp~gng+csc|Subject,data=agl,
model=FLXglm(family=binomial),k=2)
Please specify the model using the following command:
mod1 - flexmix(cbind(resp, 1-resp) ~ gng+csc | Subject, data=agl,
Dear R-helpers,
I am using xtable package to prepare a Latex code of some R tables.
Is this possible to have a greek symbols in xtable cells?
How can I get for example a string of : $\Delta$
$\Delta$
[1] $Delta$
And string: $\\Delta$
[1] $\\Delta$
Gives a latex aoutput
Folks,
I'm writing a gui wrapper around RODBC (having abandoned RMySQL) and some
custom analysis scripts. I'm hoping that I could get some advice.
All of my users will have access to a database server. I am going to limit
their access to certain databases--a sizeable warehouse in a SAN, and
Agreed. What Petr Pikal wrote, works exactly, so thank you all!
Best,
Markku
It would be helpful if you could be more specific of what exactly
you'd like to compute. Have a look also at the posting guide available
at:
http://www.R-project.org/posting-guide.html
Best,
Dimitris
On Mon, 16 Apr 2007, Pete Cap wrote:
Folks,
I'm writing a gui wrapper around RODBC (having abandoned RMySQL) and
some custom analysis scripts. I'm hoping that I could get some advice.
All of my users will have access to a database server. I am going to
limit their access to certain
Philippe Grosjean phgrosjean at sciviews.org writes:
Hello,
For a paper published in 2007, and submitted in April 2005, this is
still surprising. If I my calculation is correct, in 2004, they would
have used R 2.2.x, or something,... not 1.9.1?
Anyway, does someone know if there is a
On Mon, 16 Apr 2007, Philippe Grosjean wrote:
Hello,
I just read this paper, and was surprised by:
This study performs a comparison of the latest versions of nine
different statistical software packages using StRD. These packages are
SAS (9.1), SPSS (12.0), Excel (2003), Minitab (14.0),
Sean
Both Bill V and Peter D are right regarding traditional
repeated-measures ANOVA that assumes equality of the variance-covariance
matrices across groups and compound symmetry of the covariance matrix.
However, there is a multivariate approach to repeated measures that does
not require the
Depending on the assumptions you make about the error sturcture of your
model there are a couple of options:
1). Use transformations and lm (assumes lognormal multiplicitive
errors):
fit - lm( log(Y) ~ log(X) )
This finds log(bo) and b1
2). Use nonlinear least squares (assumes normal additive
Hi everyone.
I have to deal with a time series data and I'm looking for some
methodological help. The time series starts in January 1992 and until
April 1996 (52 observations) the values are ok (there is some trend
and seasonality in the data, but nothing special). From Mai 1996 to
June 1997 no
Unfortunately, it is not.
I am aware of Win Emacs, but it is not at all easy, sometimes almost
impossible to convince the admin(s) to install software. I was looking
(in the best case) for some R-built-in routine to do this.
Petr
Richard M. Heiberger napsal(a):
ESS is the answer for Windows.
emacs _installation_ does not require Windows Administrator privileges.
emacs can run from a CD, or USB, or a user area.
It is possible that your IT department enforces a policy by which
the emacs.exe won't _run_ if it's name is not on the policy list.
In that case, from the R Console, click on
Thomas Lumley tlumley at u.washington.edu writes:
On Mon, 16 Apr 2007, Philippe Grosjean wrote:
Hello,
[snip]
No -- there is a new x.y.0 twice per year. Checking the r-announce
archives shows that 2.0.0 came out in October 2004 and 1.9.1 in June 2004.
Describing 1.9.1 as the
I have a vector of character strings such as
mainnames-c(CAD,AUD) and another vector say
checknames-c(CAD.l1,AUD.l1,JPY.l1,EUR.l1,CAD.l2,AUD.l2,JPY
.l2,EUR.l2)
I want a new vector of character strings that is just
resultnames-c(JPY.l1,EUR.l1,JPY.l2,EUR.l2)
Because any time a name occurs in
Try this:
mainnames-c(CAD,AUD)
checknames-c(CAD.l1,AUD.l1,JPY.l1,EUR.l1,CAD.l2,AUD.l2,JPY.l2,EUR.l2)
x.match - grep(paste(mainnames, collapse=|), checknames)
x.match
[1] 1 2 5 6
checknames[-x.match]
[1] JPY.l1 EUR.l1 JPY.l2 EUR.l2
On 4/16/07, Leeds, Mark (IED) [EMAIL PROTECTED] wrote:
For one thing, get rid of the for() loop.
The following may or may not be faster, and may or may not use more
memory (an issue if x1 and x2 are long), but it's easier to read.
n - length(x1)
y - numeric(n)
tmp - x1=5 x1 = 10 x2 = 5 x2 = 10
y[tmp] - 0.631*x1[tmp]^0.55 * x2[tmp]^0.65
# and
Hi,
I'd like to simply add column-wise using Matrix objects (Csparse).
It looks like one can apply mosty any base function to these objects
(i.e., apply, colSums), but there is a nasty conversion to traditional
matrix objects if one does that.
Is there any workaround? I can see colSum listed
Hi, Simon,
Thank you very much for the helpful suggestion. But if the model is a
partially linear model, i.e.
y(x1)=a0+a1*x1+a2*x2+a3*s(x3)where x1 is the decision variable
I would like to evaluate y(x1) at a specific level of x2 and x3 (i.e. y is a
function of x1 only). In
Hello,
If it is not a problem for ODBC, you can use an SQL statemente to change the DB:
use NAME_OF_DB
Folks,
I'm writing a gui wrapper around RODBC (having abandoned RMySQL) and some
custom analysis scripts. I'm hoping that I could get some advice.
All of my users will have access to a
Lucke, Joseph F wrote:
Sean
Both Bill V and Peter D are right regarding traditional
repeated-measures ANOVA that assumes equality of the variance-covariance
matrices across groups and compound symmetry of the covariance matrix.
However, there is a multivariate approach to repeated measures
On Mon, 16 Apr 2007, Andris Jankevics wrote:
Dear R-helpers,
I am using xtable package to prepare a Latex code of some R tables.
Is this possible to have a greek symbols in xtable cells?
How can I get for example a string of : $\Delta$
$\Delta$
[1] $Delta$
And string: $\\Delta$
As Gabor said, ts series are not good for weekly data. That's why I created the
tis (Time Indexed Series) class, in the 'fame' package. You don't need the
Fame database to use them.
--
Jeff
__
R-help@stat.math.ethz.ch mailing list
I'm using R in a winXp machine. I have problems with the size of my database
(xyz), anyone know if is possible to use some kind of swap memory, in order
to expand my real RAM and use HD space for the processes?
I'm trying to make a variogram, and I don't want to sacrifice the precision
of the
I wish to set up a simple function using boxplot so that it will be
available to someone using R for Windows. (I myself use Linux.)
The way the data is organised makes it convenient to use the boxplot
function in a manner similar to this example given in the help.
mat - cbind(Uni05 =
Try this:
boxplot(data.frame(mat), names = colnames(mat))
On 4/16/07, Patrick Connolly [EMAIL PROTECTED] wrote:
I wish to set up a simple function using boxplot so that it will be
available to someone using R for Windows. (I myself use Linux.)
The way the data is organised makes it
Dear ListeRs,
I am trying to fit a heteroscedastic multilevel model using lmer{lme4-
package). Take, for instance, the (fictive) model below.
lmer(test.result ~ homework + Sex -1 + (1 | School))
Suppose that I suspect the error terms in the predicted values to
differ between men and women
Dear R-users,
I have 10 data files in gpr format (dat1.gpr,
dat10.gpr). I want to read in
these files one by one in R and then add one extra column (called log) to each
data file as below
On 4/16/2007 4:22 PM, Patrick Connolly wrote:
I wish to set up a simple function using boxplot so that it will be
available to someone using R for Windows. (I myself use Linux.)
The way the data is organised makes it convenient to use the boxplot
function in a manner similar to this example
On Mon, 16-Apr-2007 at 04:29PM -0400, Gabor Grothendieck wrote:
| Try this:
|
| boxplot(data.frame(mat), names = colnames(mat))
|
Great! Now why didn't I think of that?
--
~.~.~.~.~.~.~.~.~.~.~.~.~.~.~.~.~.~.~.~.~.~.~.~.~.~.~.~.~.~.~.~.~.~.~.~.
___Patrick Connolly
X2R is a bundle of three software libraries allowing the user to pass
structured data easily from Fortran, C/C++, or AD Model Builder to R.
An update to X2R has been sent to CRAN today and should be available at
mirrors shortly. This fixes a bug in the ADMB2R and C2R components in
which data
On Mon, 16 Apr 2007, Duncan Murdoch wrote:
On 4/16/2007 4:22 PM, Patrick Connolly wrote:
I wish to set up a simple function using boxplot so that it will be
available to someone using R for Windows. (I myself use Linux.)
The way the data is organised makes it convenient to use the boxplot
Dear List,
My first R function is a rip-off bagging algorithm from pg. 138 of
Everitt and Hothorn's Handbook of Statistical Analyses using R
(HSAUR). I'm using recursive partitioning to develop a set of useful
variables in diagnosing ADHD.
I'm running this in ESS in XEmacs 21.4.19, R 2.4.1 on
Hi R-helpers,
I wish to produce frequency histograms of all of the variables in my
dataframe (except some identifying variables).
I have tried
hist(dataframe[,3:20])
to produce histograms of the 3rd through 20th variables in my dataframe, but
R doesn't like that.
Could anyone provide a
try:
win.metafile(file=output.wmf)
lapply(dataframe[3:20], hist)
dev.off()
On 4/16/07, [EMAIL PROTECTED] [EMAIL PROTECTED] wrote:
Hi R-helpers,
I wish to produce frequency histograms of all of the variables in my
dataframe (except some identifying variables).
I have tried
Hi Jim, thanks for your help. That looks like it might work, but a couple of
things...
1) The resulting 73 (in my case) histograms will be named by their variable
number, not by the variable name contained in the first row. Any way to
include the variable name in the resulting histogram?
2) How
I forgot that win.metafile can only output one graph. This will write
separate files:
lapply(3:20, function(.ind){
win.metafile(filename=paste(file., .ind, .wmf, sep=''))
hist(dataframe[[.ind]], main=paste(dataframe[[, .ind, ]],
sep=''), xlab=)
dev.off()
})
On 4/16/07, [EMAIL
This will put the colnames on the histograms:
lapply(3:20, function(.ind){
win.metafile(filename=paste(file., .ind, .wmf, sep=''))
hist(dataframe[[.ind]], main=colnames(dataframe)[.ind], xlab=)
dev.off()
})
On 4/16/07, [EMAIL PROTECTED] [EMAIL PROTECTED] wrote:
Hi Jim, thanks for
__
R-help@stat.math.ethz.ch 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 commented, minimal, self-contained, reproducible code.
On Apr 16, 2007, at 10:41 AM, John Kane wrote:
--- Charilaos Skiadas [EMAIL PROTECTED] wrote:
It is not an export option, it is a save as
option. I don't have
a 14 to check, but on a 15 I just go to File - Save
As, and change
the Save as type field to Comma Delimited
(.*.csv). (I suppose
Hi, Joel:
I cannot reproduce the error w/o the data called probands but I assume
the line in bagger function
bootsamples - rmultinom (length(trees), n,
rep(1,n)/n)
might not succeed so that it complains bootsamples not found.
try to add one line like
Suppose you have a vector of data in x and response values in y. How
do you plot together both the points (x,y) and the curve that results
from the fitted model, if the model is not y ~ x, but a higher order
polynomial, e.g. y~poly(x,2)? (In other words, abline doesn't work
for this case.)
Thanks again Jim, I appreciate your time. I've been trying to debug the
code, but am running into a wall. I'm getting a syntax error after the line
containing the hist function. Here's the R session, any ideas?
Also, I'd like to be able to have each histogram use the same x-axis breaks
(0, 1,
On 4/16/07, Paul Lynch [EMAIL PROTECTED] wrote:
Suppose you have a vector of data in x and response values in y. How
do you plot together both the points (x,y) and the curve that results
from the fitted model, if the model is not y ~ x, but a higher order
polynomial, e.g. y~poly(x,2)? (In
On 4/16/07, Paul Lynch [EMAIL PROTECTED] wrote:
Suppose you have a vector of data in x and response values in y. How
do you plot together both the points (x,y) and the curve that results
from the fitted model, if the model is not y ~ x, but a higher order
polynomial, e.g. y~poly(x,2)? (In
Dear R-experts,
This should be a simple question, but I couldn't find a way yet.
I would like to create a function that accepts a vector, creates a
data.frame whose colname is that vector. For example, I want a function
temp.func() such that:
name = c( mike, john, steve)
result =
On 4/16/2007 8:19 PM, jiho.han wrote:
Dear R-experts,
This should be a simple question, but I couldn't find a way yet.
I would like to create a function that accepts a vector, creates a
data.frame whose colname is that vector. For example, I want a function
temp.func() such that:
name =
name - letters[1:4]
data.frame(name, stringsAsFactors=FALSE)
name
1a
2b
3c
4d
__
R-help@stat.math.ethz.ch mailing list
https://stat.ethz.ch/mailman/listinfo/r-help
PLEASE do read the posting guide
Dr. Paul Murrell will offer his Graphics in R course online at
statistics.com May 4 June 1. Participants can ask questions and exchange
comments with Dr. Murrell via a private discussion board throughout the period.
This course teaches you how to produce publication-quality statistical
Hi All R Experts,
Please help me. Query is given at the end of this R Session
+start of R Session+++
SCORES-read.csv(C:/SCORE.csv,header=TRUE,sep=,)
SCORES
SrNo Player Score
1 1 A12
2 2 B23
3 3 C34
4 4 A54
+
17
rownames is the command you want, I think.
rownames(A) - seq(1:8)
Try that.
Pat
On 4/16/07, [EMAIL PROTECTED] [EMAIL PROTECTED] wrote:
Hi All R Experts,
Please help me. Query is given at the end of this R Session
+start of R Session+++
82 matches
Mail list logo
