Dear R users,
I have a data frame (test) including five columns of upper (numeric), lower
(numeric), observed (numeric), best_sim (numeric) and stname (factor with 80
levels, each level with different length). Now I would like to write a short
program to draw one graph as follow for each level
Hello,
See
?Devices
?dev.print
e.g.
X11() # opens a new screen device
plot(1:10) # plots points on the device
dev.print(jpeg, [args]) # prints from the X11 device to a jpeg file
dev.off() # close X11 device
or
jpeg([args]) # opens a jpeg device (file) for plotting
plot(1:10) # plots points
Creating more than one graphic windows is, as far as I know, not
possible in R. But it's no problem to run a script which create multiply
jpeg's. See ?jpeg on how to create jpeg's.
Cheers,
Thierry
ir. Thierry
Hi all,
I am doing a two-sided DW test:
H0: rho = 0
H1: rho =/= 0
My understanding is that most test statistics tables are one-sided. It's the
way they created the table.
So from online, by doing Googling, I found a bunch of DW tables for
Confidence Level 5%.
Those tables can answer my
Hello sir;
I use the function heatmap.2 to draw a heatmap of microarray data,which
consists of logratios.
Q1
But the lengend shows the Z score and the corresponding color.But I want the
legend to show the logratios and the corresponding color.How can I do it?
Q1
How can I control that cluster
On Wed, 2007-03-07 at 09:39 +0100, Faramarzi Monireh wrote:
Dear R users,
I have a data frame (test) including five columns of upper (numeric),
lower (numeric), observed (numeric), best_sim (numeric) and stname
(factor with 80 levels, each level with different length). Now I would
like to
Hi all,
I have some residuals from regression, and i suspect they have correlations
in them...
I am willing to cast the correlation into a ARMA(p, q) framework,
what's the best way to identify the most suitable p, and q, and fit ARMA(p,
q) model and then correct for the correlations in
On Wed, 2007-03-07 at 15:33 +1000, Peter Dunn wrote:
Hi all
I love Sweave; use it all the time.
But I recently received a new computer, and ever since I
have had a problem I've never seen before.
For example, I place the following in my Snw file:
Try this in the preamble of your Snw
Buongiorno a tutti
mi chiamo Diego Pettena e sono uno studente universitario iscritto alla
facoltà di economia. Ho da poco iniziato ad utilizzare il sitema R per
l'analisi statistica, tuttavia con la presente non intendo chiederVi
nulla riguardo al linguaggio o all'utilizzo del programma. Per
Michael:
I am doing a two-sided DW test:
H0: rho = 0
H1: rho =/= 0
My understanding is that most test statistics tables are one-sided. It's the
way they created the table.
...because rho 0 is the alternative of interest in most applications.
So from online, by doing Googling, I found a
Dear Alberto,
please note that special characters (eg a space character) in the Excel
sheet names mess up the simple way of querying provided by sqlFetch.
If you have a regular case of all sheets like Sheet1:
plan1 - sqlFetch(channel,Sheet1) # should work
But if you have Sheet 1 (
Michael wrote:
Hi all,
I am doing a two-sided DW test:
H0: rho = 0
H1: rho =/= 0
My understanding is that most test statistics tables are one-sided. It's the
way they created the table.
So from online, by doing Googling, I found a bunch of DW tables for
Confidence Level 5%.
On 07-Mar-07 ONKELINX, Thierry wrote:
Creating more than one graphic windows is, as far as I know, not
possible in R.
It is if you are running R on Linux (in which the X Windowing
System -- or X -- is the standard graphical system). As it
says in ?X11:
This can only be done on machines that
The down side to R's factor solution:
The numerical values of factors are always 1 to number of levels. Thus, it
can be tough and requires great care to work with studies that have both
numerical values different from this and value labels. This situation is
currently not well-supported by R.
Matthias Kohl wrote:
in distr you can do:
library(distr)
N - Norm(mean = 1, sd = 2)
p(N)(0.5)
r(N)(100)
!!! not: p(N, 0.5) or r(N, 100) !!!
A detailed description of package distr is given in package distrDoc.
library(distrDoc)
vignette(distr)
Thanks!!! This is almost perfect.
Gavin Simpson [EMAIL PROTECTED] writes:
On Wed, 2007-03-07 at 15:33 +1000, Peter Dunn wrote:
But I recently received a new computer, and ever since I
have had a problem I've never seen before.
For example, I place the following in my Snw file:
Try this in the preamble of your Snw file:
Use the iplots package. The interactive querying works effectively
for largish datasets (over 50,000 cases) where sunflower plots may not.
Antony Unwin
Professor of Computer-Oriented Statistics and Data Analysis,
University of Augsburg,
86135 Augsburg, Germany
Tel: + 49 821 5982218
Ted Harding wrote:
Creating more than one graphic windows is, as far as I know, not
possible in R.
But, as to whether/to what extent X or equivalent is available for
MS Windows, that is another question on which I have no expertise.
X11() seems to work for Windows XP.
Alberto Monteiro
Hi R-users,
when carrying out a multiple regression, say lm(y~x1+x2), we can use an
anova of the regression with summary.aov(lm(y~x1+x2)), and afterwards
evaluate the relative contribution of each variable using the global Sum of
Sq of the regression and the Sum of Sq of the simple regression
Hello all,
I'm a new user of R, experienced with Octave/MATLAB and therefore
struggling a bit with the new syntax.
One of the easy things in Octave or MATLAB is to plot multiple lines or
sets of points by using a matrix where either the columns or the rows
contain the y-values to be plotted.
Hi
see matplot, matlines.
or use forbidden for cycle.
for (i in 1:n) lines(x,y[,i], col=i)
or if you want to use more colours use built in rainbow, topo.colors
or generate your own set.
Regards
Petr
On 7 Mar 2007 at 12:30, Joseph Wakeling wrote:
Date sent: Wed, 07 Mar 2007
On Wed, 2007-03-07 at 12:30 +, Joseph Wakeling wrote:
Hello all,
I'm a new user of R, experienced with Octave/MATLAB and therefore
struggling a bit with the new syntax.
One of the easy things in Octave or MATLAB is to plot multiple lines or
sets of points by using a matrix where
Dear AIX useR's.
Will you test it in various environment?
There is not xlf in the environment where I can approach.
cf.
32bit, 64bit, --enable-R-shlib, --enable-BLAS-shlib ...various cases!
http://prs.ism.ac.jp/%7enakama/AIX/
--
EI-JI Nakama [EMAIL PROTECTED]
\u4e2d\u9593\u6804\u6cbb
Dear Alberto,
It is better to assign a name to an area of data and not to use
the name of the sheet, because this can have graphs and other data.
If you this interested can see:
http://tarwi.lamolina.edu.pe/~fmendiburu/Rsolutions.htm
I hope that this also helps,
Felipe.
-Original
Hi,
My name is Cristina. I'm interested in studying which continuos predictor
variables (such as grooming received, rank, etc.) affect grooming given, as
well a continuos variable. I'm having problems finding an appropriate family
distribution to fit the GLMM I'm doing. The response variable,
David,
I wouldn't give up on windows so fast. Many people have gotten the 3Gb
switch to work. One used to have to modify the header of the Rgui.exe
program to use the switch, but now the binary comes ready for that, so
its really quite easy. I would like to hear more about why its not
working
Try na.locf from the zoo package and then use merge with specified suffixes:
library(zoo)
f - function(x) {
rownames(x) - NULL
merge(x, na.locf(x[-1], na.rm = FALSE), by = 0, suffixes = c(, .by))[-1]
}
do.call(rbind, by(x, x$id, f))
On 3/7/07, Jon Olav Vik [EMAIL PROTECTED] wrote:
Dear
Ulrike Grömping wrote:
The down side to R's factor solution:
The numerical values of factors are always 1 to number of levels. Thus, it
can be tough and requires great care to work with studies that have both
numerical values different from this and value labels. This situation is
Alberto Monteiro wrote:
Ted Harding wrote:
Creating more than one graphic windows is, as far as I know, not
possible in R.
But, as to whether/to what extent X or equivalent is available for
MS Windows, that is another question on which I have no expertise.
X11() seems to work for Windows
Jason Barnhart wrote:
- Original Message -
From: John Kane [EMAIL PROTECTED]
To: lamack lamack [EMAIL PROTECTED]; R-help@stat.math.ethz.ch
Sent: Tuesday, March 06, 2007 2:13 PM
Subject: Re: [R] R and SAS proc format
--- lamack lamack [EMAIL PROTECTED] wrote:
Dear all,
Hello,
I was wondering if anybody could help me with this?
I have plotted an acf function for a time series and am very happy with it.
Now I am interested in calculating for myself the two values for the confidence
intervals that are plotted on the graph of the acf.
The confidence intervals do
Dear people,
I've a problem in doing a power calculation. In Fryer and Nicholson
(1993), ICES J. mar. Sci. 50: 161-168 page 164 an example is given with
the following characteristics
T=5, points in time
R=5, replicates
Var.within=0.1
q=10, a 10% increase per year
The degrees of freedom for the
Hi,
I have a little problem with the installation of a new packages. The
installation of R software is correct, but my server required authentication
for use it, and for load a new package directly from R it is not possible.
Is there a code or process for server authentication (put my login and
I have used RODBC to get the database i can view the tables in RGUI-2.4.1
how can i query the records in R
i tried with sqlQuery
need some help
JJ
--
Lecturer J. Joshua Thomas
KDU College Penang Campus
Research Student,
University Sains Malaysia
--
Lecturer J. Joshua Thomas
KDU College
Dear R users,
I thought this might be useful for users of the R package HOPACH
(Windows version). Using R 2.4.1 in Windows XP, I found that the
functions: makeoutput, boot2fuzzy and hopach2tree is not recognized by
HOPACH versions 1.4.3 and 1.8.0. However, it works perfect with HOPACH
version
library(RGtk2)
library(rattle)
rattle()
click the ODBC option it as the DSN i am a bit confused with this i already
put my *.mdb file in C:drive
i try put the DSN name as Microsoft Access driver, in the appropriate text
box but i couldnt locate the table
i tried the other way round open- locate
Hmm, this does not seem to be over-documented :-)
But try
?plot.acf
and
getAnywhere(plot.acf)
Then you can find in the code how the values are actually calculated.
Petr
Gladwin, Philip [CIB-FI] napsal(a):
Hello,
I was wondering if anybody could help me with this?
I have plotted an acf
On 3/7/07, Peter Dalgaard [EMAIL PROTECTED] wrote:
Jason Barnhart wrote:
- Original Message -
From: John Kane [EMAIL PROTECTED]
To: lamack lamack [EMAIL PROTECTED]; R-help@stat.math.ethz.ch
Sent: Tuesday, March 06, 2007 2:13 PM
Subject: Re: [R] R and SAS proc format
---
Dear R-help community,
I have had trouble in the past installing the latest version of R: we got the
errors shown below (the computer specifications and version of R are below
that). Does anybody have tips for compiling the latest version of R so that I
can avoid these errors?
configure
make
Hi,
On 3/7/07, j.joshua thomas [EMAIL PROTECTED] wrote:
I have used RODBC to get the database i can view the tables in RGUI-2.4.1
how can i query the records in R
i tried with sqlQuery
need some help
A bit of code showing what you did would help others to track down what
may have gone
Dear all, Why the transform function does not accept two statistics
functions?
a = data.frame(matrix(rnorm(20),ncol=2))
transform(a,M.1=mean(X1),M.2=mean(X2)) # does not works
#while:
transform(a,M.1=mean(X1),M2=log(abs(X2))) #works
Best regards
JL
Hi,
I'm trying to aggregate date values using the aggregate function. For example:
aggregate(data,by=list(weekdays(LM),months(LM)),FUN=length)
I would also like to aggregate by year but there seems to be no
years() function.
Should there be one? Is there any alternative choice?
Also, a hours()
Gavin Simpson wrote:
You want maplot here. See ?matplot but here is an example:
Great! Thanks to you and Petr for pointing this out, it's exactly what
I wanted. Petr's other suggestions look interesting and I'll explore
them at length later.
Note the changed axis range in the right-hand
Meesters, Erik wrote:
Dear people,
I've a problem in doing a power calculation. In Fryer and Nicholson
(1993), ICES J. mar. Sci. 50: 161-168 page 164 an example is given with
the following characteristics
T=5, points in time
R=5, replicates
Var.within=0.1
q=10, a 10% increase per year
The
Hi,
Is there a smart way in the R graphs to create a line that is broken in
intervals based on the indicator given below.
following is a small test graph
Location,indicator,otherinfo
1.2,1,2.2
2.5,1,2.5
3.7,1,2.3
20.1,2,4.3
22.5,2,5.2
25.0,2,3.4
27.3,2,2.2
35.1,3,3.4
37.0,3,7.2
38.0,3,6.1
On Wed, 2007-03-07 at 15:11 +, Joseph Wakeling wrote:
Gavin Simpson wrote:
You want maplot here. See ?matplot but here is an example:
Great! Thanks to you and Petr for pointing this out, it's exactly what
I wanted. Petr's other suggestions look interesting and I'll explore
them at
Dear Monireh,
try using lattice:
library(lattice)
set.seed(1234)
dat - data.frame(months=rep(1:10,80),upper = rnorm(800)+1,
lower = rnorm(800)-1,
observed = rnorm(800), best.sim = rnorm(800),
stname = factor(gl(80, 10)))
jpeg(filename =
Hello all,
I was wondering if anyone could help me formulate a Two-way ANOVA for
unbalanced multiple sample data?
We have a new study method aimed to help students to study for tests
using computers. (I am a computer scientists, hence my
soon-to-be-apparent lack of statistical knowledge).
To
The variances of the random effects and the residual variances are given
by the summary function. Maybe VarCorr or varcomp gives you the answer
you are looking for:
library(nlme)
library(ape)
?VarCorr
?ape
JR
El mié, 07-03-2007 a las 13:09 +0100, Berta escribió:
Hi R-users,
when carrying out
Jenny Barnes wrote:
Dear R-help community,
I have had trouble in the past installing the latest version of R: we got the
errors shown below (the computer specifications and version of R are below
that). Does anybody have tips for compiling the latest version of R so that I
can avoid
Thanks José Rafael, I will try with library(ape) (at the moment I cannot
load it).
VarCorr gives the variance estimates for the random effect and the error
terms. However, what I am looking for is a measure of the explained
proportion of variance, such as it is R2 in regression models, and
Thanks for the tips, Roger.
fyi: When I added /3GB to the boot.ini, the resulting desktop was incomplete
and locked - no chance to try starting R. Searching the web lead me to
believe that this was possibly a dead-end, so I abandoned this effort. Any
hints on getting this to work, anyone?
If you insert an NA (or row of NA's) into the data at each place you
want a break (after indicator increases), then the regular plot with
type='l' will break the line for you.
Is this what you want?
--
Gregory (Greg) L. Snow Ph.D.
Statistical Data Center
Intermountain Healthcare
[EMAIL
Hi
you shall probably cooperate with segments, so you need to extract
start and end points for your lines e.g.
x-c(1:6, 10:15,20:25)
y-rep(c(1,2,3), each=6)
plot(x,y, type=l)
plot(x,y)
segments(sapply(split(x,y), min),1:3, sapply(split(x,y),max),1:3)
Regards
Petr
On 7 Mar 2007 at
I was finally able to get R to 'configure', 'make', and 'install' on Mandriva
2007. Itried to install gnomeGUI and I received an error. See below. At what
step do I make R a shared library? Where did I go wrong?
Ross
==
downloaded 74Kb
*
Hi,
New to both R and SEM, so this may be a very simple question. I am
trying to run a very simple path analysis using the sem package.
There are 2 exogenous (FARSCH, LOCUS10) and 2 endogenous (T_ATTENT,
RMTEST) observed variables in the model. The idea is that T_ATTENT
mediates the
Hello,
Has anyone successfully fit empirical data to a tCopula using the fitCopula
function? If so, are there ways to pick intelligent starting values to
avoid the errors such as a minor matrix not being positive definite and the
initial value of 'vmmin' not being finite?
I've been able to
Dear r-helpers,
I can run the examples on the mcsamp help page. For example:
M1 - lmer (y1 ~ x + (1|group))
(M1.sim - mcsamp (M1))
fit using lmer,
3 chains, each with 1000 iterations (first 500 discarded)
n.sims = 1500 iterations saved
lamack lamack wrote:
Dear all, Why the transform function does not accept two statistics
functions?
a = data.frame(matrix(rnorm(20),ncol=2))
transform(a,M.1=mean(X1),M.2=mean(X2)) # does not works
#while:
transform(a,M.1=mean(X1),M2=log(abs(X2))) #works
It's a variation of this
It's not the correlation as such that is the problem; it's because you
only have 10 degrees of freedom available with four observed
variables, and you are estimating 10 parameters, which is why you get
a chi square of zero. When you remove any one free parameter (such as
the correlation), the
Hi Greg,
Thank you for your response and a previous posting about Macros in R.
Thank you also to Ken Knoblouch (Ken had the same idea as Greg's and
Peter Pikal (who proposed the use of segments function).
There is only a technical specific that when applying max function to
find the limit of y
Dear David and Ista,
I haven't looked at this model carefully, but the fact that the df are
0 suggests that the model is just-identified and therefore necessarily
perfectly reproduces the covariances among the observed variables.
Removing a parameter would over-identify the model, making possible
Hi Everyone,
I'm so confused. I've been trying to save data to a table but I keep
getting an error that says the table
does not exist and at other times saying that it does. So here are some
statements:
sqlQuery(channel, select top 1 * from
TestDB.[SILICON\\holouis1].clep_tier_shift)
State
Hi
I have two dataframes
names(DF1) = c(id, val1, val2);
names(DF2) = c(id2);
Ids in DF2 are a complete subset of those in DF1
How can I extract entries from DF1 where id NOT IN
DF2.
I tried setdiff(DF1, DF2); setdiff(DF1$id, DF2$id),
etc.
Although the latter eliminates the ids as required, I
Hi,
I have a little problem with the installation of a new packages. The
installation of R software is correct, but my server required authentication
for use it, and for load a new package directly from R it is not possible.
Is there a code or process for server authentication (put my login and
David,
Here is what my boot.ini file looks like:
[boot loader]
timeout=5
default=multi(0)disk(0)rdisk(0)partition(1)\WINDOWS
[operating systems]
multi(0)disk(0)rdisk(0)partition(1)\WINDOWS=Microsoft Windows XP
Professional /noexecute=optin /fastdetect /3gb
The easiest way to edit the boot.ini
try something along these lines (untested):
DF1[DF1$id %in% DF2$id2, c(val1, val2)]
DF1[!DF1$id %in% DF2$id2, c(val1, val2)]
I hope it helps.
Best,
Dimitris
Dimitris Rizopoulos
Ph.D. Student
Biostatistical Centre
School of Public Health
Catholic University of Leuven
Address:
Hi Aldi,
Yet another way to do what you want. 'd' is your data frame. You cannot easily
modify the attributes of the lines however...
# just to set up coordinates
plot( d$Location, d$otherinfo, type=n)
# plot the lines
invisible(by(d, d$indicator, function(dd) lines(dd$Location,
I`m doing some functions on C that gives me the x and y coordinates.
I`d like to now how I can get these coordinates (both are a vector of
number) on R to that I can make a graphic.
I`ve already made a package with my functions, so I just wanna how
about how to get the coordinates.
Thanks,
Luis Garavito [EMAIL PROTECTED] writes:
Hi,
I have a little problem with the installation of a new packages. The
installation of R software is correct, but my server required authentication
for use it, and for load a new package directly from R it is not possible.
Is there a code or process
Dear all,
I wanted to use ATLAS to increase the computational speed of my
program. I have installed R in a straightfoward way with .exe file
(not building from source).
However, after getting the P4 ATLAS Rblas.dll file from the CRAN
directory bin/windows/contrib/ATLAS, and replacing the
Luis,
posting a message several times usually won't help you to get answers...
You did not tell us what kind of error you obtain, what packages did you
try to install and how, what operating system are you using etc. To be
honest, I have no idea what kind of server can be involved - R does not
Dear r-help users,
I have the following simple problem: Reading data from a file. The
file is a .txt file exported (save as...) from Excel (see below for
an example). The Excel file consists of two header rows (first row
consists of ticker symbols of stocks, the second row consists of
Dear R list,
I have a question in R, it could be very simple, but I don't know how to do it?
for example:
I assign 6 to x in beginning of of my R script code
x-6
..
After many line code, I forget using x variable before, I use x
again, and do assignment like this
x-45
x
[1] 45
then
More problems. If I run
sim(fm1 - lmer(Reaction ~ Days + (Days|Subject), sleepstudy))
from the lmer() help page.
I get the error
Error in mvrnorm(n.sims, bhat[j, ], V.beta) :
'Sigma' is not positive definite
On Mar 7, 2007, at 1:30 PM, Michael Kubovy wrote:
Dear r-helpers,
I can run
Sérgio Nunes snunes at gmail.com writes:
Hi,
I'm trying to aggregate date values using the aggregate function. For example:
aggregate(data,by=list(weekdays(LM),months(LM)),FUN=length)
I would also like to aggregate by year but there seems to be no
years() function.
Should there be
Hi Heloise,
there's a manual online to help you with this - see here
http://cran.r-project.org/doc/manuals/R-exts.html
Cheers
Andrew
On Wed, Mar 07, 2007 at 05:29:10PM -0300, Heloise Mattos wrote:
I`m doing some functions on C that gives me the x and y coordinates.
I`d like to now how I can
From the help of weekdays:
Note:
Other components such as the day of the month or the year are very
easy to compute: just use 'as.POSIXlt' and extract the relevant
component.
Yet another option:
help(package=chron)
JR
El mié, 07-03-2007 a las 15:35 +, Sérgio Nunes
dear all,
how can i get a vector that shows the number of the column of matrix
that contains the maximum of the row ??
can´t believe in need a loop for this...
i have a 100 x 3 matrix and want to get a 100 x 1 vector with values
1,2,3 .
there must be a simple solution. i just cannot find
dear all,
how can i get a vector that shows the number of the column of
matrix that contains the maximum of the row ??
can´t believe in need a loop for this...
i have a 100 x 3 matrix and want to get a 100 x 1 vector with
values 1,2,3 .
there must be a simple solution. i just
?apply
?which.max
m - matrix(rnorm(12),nrow=4)
m
apply(m,1,which.max)
Petr
bunny , lautloscrew.com napsal(a):
dear all,
how can i get a vector that shows the number of the column of matrix
that contains the maximum of the row ??
can´t believe in need a loop for this...
i have a
apply(yourMatrix,1,which.max)
Bert Gunter
Nonclinical Statistics
7-7374
-Original Message-
From: [EMAIL PROTECTED]
[mailto:[EMAIL PROTECTED] On Behalf Of bunny ,
lautloscrew.com
Sent: Wednesday, March 07, 2007 2:12 PM
To: r-help@stat.math.ethz.ch
Subject: [R] hwo can i get a vector
can´t believe in need a loop for this...
No, you don't ;)
want to get a 100 x 1 vector
Has each row unique values? If yes:
mat - matrix(rnorm(300),nr=100)
vet - apply(mat,1,function(x) {return(which(x==max(x)))})
scionforbai
__
Andrew Robinson has gently chided me for not including more
information. So here goes:
R version 2.4.1 (2006-12-18)
powerpc-apple-darwin8.8.0
locale:
C
attached base packages:
[1] stats graphics grDevices utils methods base
other attached packages:
foreign car
Check out which.max
Peter Alspach
-Original Message-
From: [EMAIL PROTECTED]
[mailto:[EMAIL PROTECTED] On Behalf Of bunny
, lautloscrew.com
Sent: Thursday, 8 March 2007 11:20 a.m.
To: R-help@stat.math.ethz.ch
Subject: Re: [R] hwo can i get a vector that...
dear all,
Michael wrote:
Hi all,
I have some residuals from regression, and i suspect they have correlations
in them...
I am willing to cast the correlation into a ARMA(p, q) framework,
what's the best way to identify the most suitable p, and q, and fit ARMA(p,
q) model and then correct for the
Decompose your code into small understandable functions.
On 3/7/07, Aimin Yan [EMAIL PROTECTED] wrote:
Dear R list,
I have a question in R, it could be very simple, but I don't know how to do
it?
for example:
I assign 6 to x in beginning of of my R script code
x-6
..
After many
___
max.col sounds like what you're after.
--
Hong Ooi
Senior Research Analyst, IAG Limited
388 George St, Sydney NSW 2000
+61 (2) 9292 1566
-Original Message-
From: [EMAIL PROTECTED] [mailto:[EMAIL
Hello R users,
Problem...I do not understand how to use aggregate,by, or the
appropriate apply to perform a function on data with more than one
factor on unbalanced data...
I have a data frame in the long format that does not contain balanced
data. The ID is a unique identifier corresponding
I would have thought that a csv file written out by Excel would have looked
like this for your data:
ADS GY Equity, ,ALV GY Equity,
Date, Px Last,Date , Px Last
07/02/04, 41,395 ,07/01/31 , 130,234
07/02/05, 42,134
Here is one way of doing it:
# create the rows for each unique combination
x.split - split(seq(nrow(mydata)), list(mydata$time, mydata$treatment,
+ mydata$expREP, mydata$techREP), drop=TRUE)
# now go through the list of indices and add the median
mydata$Y50 - 0 # add the dummy median
Hi,
I am really new with R, so I don't know anything about it. I have
written a script (attached) which tries to do really basic stuff (such
as computing basic statistics and basic plots). When I try to plot a
histogram and pairs, for example, I get the following message:
source(project.R)
XLSolutions Corporation is proud to announce our April 2007 R/S:
Fundamentals and Programming Techniques - @ 6 USA locations:
http://www.xlsolutions-corp.com/Rfund.htm
Ann Arbor, MI *** April 23-24, 2007
Raleigh, NC *** April 19-20, 2007
Princeton, NJ *** April 26-27 , 2007
Hartford, CT ***
I have the same question that Eusebio had:
Is there a function similar to sink that redirect also R code to a
file that
is:
sink(R001)
x - c(2,-6,-4,8,5,4,1,3,4,-9,0,1)
A - matrix(x, ncol=3)
A
A.prima - t(A)
A.prima
dim(A)
dim(A.prima)
sink()
create a file R001 with contents:
Thanks to Gavin Simpson and Patrick Drechsler:
The solution was simple. (No idea how I would have
discovered these answers without this R mailing list!)
Many thanks.
P.
Gavin Simpson [EMAIL PROTECTED] writes:
On Wed, 2007-03-07 at 15:33 +1000, Peter Dunn wrote:
But I recently received a new
Dear Alan,
I think that podria to be of utility the function tapply.stat () of the package
agricolae.
see ?tapply.stat
Regards,
Felipe.
for example:
library(agricolae)
attach(mydata)
set1-tapply.stat(mydata[,2:5],Y,median)
set2-tapply.stat(time,Y,function(x) median(x))
Hi,
I'm hoping someone has a suggestion for handling a simple problem. A
client gave me a comma separated value file (call it x.csv) that has
an id and name and address for about 25,000 people (25,000 records).
I used read.table to read it, but then discovered that there are
#!/bin/sh
echo 'a=${1}; b=${2}; source(myRcodeUsing_a_and_b.R); ' | R
--vanilla --quiet
and you can run this from shell like:
./callR valueOfa valueOfb
Best, Oleg
akintayo holder wrote:
Hi,
Does any one know if it is possible to create an R script that can use
command line parameters. I
How do you define a carriage return in the middle of a line if a carriage
return is also used to delimit a line? One of the things you can do is to
use 'count.fields' to determine the number of fields in each line. For
those lines that are not the right length, you could combine them together
Hello:
I am sure this question was dealt with several years ago. Is the function
vmmin() available from Rmath Standalone? If not is it possible to call
optim() or nlm() from Rmath in C. Thank you.
Mervyn
__
R-help@stat.math.ethz.ch mailing
1 - 100 of 112 matches
Mail list logo