I think what you are asking for is not a tilde. That is a raised symbol
(an accent), and not as in TeX's $\sim$ . It is character 126 in the
Adobe Symbol encoding (Adobe's name is '\similar'), so one way is
expression(X*symbol(\176)*N(mu, sigma^2))
There are others, and in most
Hi!
Does anyone know a R package for multivariate (not univariate!) shape
preserving interpolating splines. For example spline must preserve
monotonocity in both directions, etc...
Tine Mlač
__
R-help@r-project.org mailing list
Hello!
I have an maybe a to a simple question but I cant get it right.
I have a dataframe with in one column an id-variable and in another a
year-variable. One id-number can occur several years.
I have sorted the dataframe on id then on year so the same id-number is sorted
by year with the
Hi Ullrich,
The model is
RT.aov - aov(RT~Cond + Error(Subj/Cond), WMU3C)
I understand that TukeyHSD only works with an aov object, but that
RT.aov is an aovlist object.
You want to use lme() in package nlme, then glht() in the multcomp package.
This will give you multiplicity adjusted
If df is your dataframe, try
aggregate(df$Year,list(df$Id),min)
--- On Fri, 30/5/08, Jojje Andersson [EMAIL PROTECTED] wrote:
From: Jojje Andersson [EMAIL PROTECTED]
Subject: [R] Aggregate by First case
To: r-help@r-project.org r-help@r-project.org
Received: Friday, 30 May, 2008, 5:13 PM
Jojje Andersson:
I have a dataframe with in one column an id-variable and in another a
year-variable. One id-number can occur several years. I have sorted the
dataframe on id then on year so the same id-number is sorted by year with
the first occurens at top.
Now I want to make a subset of
Dimitris Rizopoulos:
dat[!duplicated(dat$Id), ]
Very elegant. I didn't know about 'duplicated()' before now.
Seems like a very useful function.
--
Karl Ove Hufthammer
__
R-help@r-project.org mailing list
https://stat.ethz.ch/mailman/listinfo/r-help
On Thu, 29-May-2008 at 11:25PM -0400, Charilaos Skiadas wrote:
On May 29, 2008, at 11:02 PM, [EMAIL PROTECTED] wrote:
Suppose I have a plot
plot(1:10, pch = )
And I want some text to indicate a Normal distrubition. I could do
this:
text(5, 6, substitute(XN(mu, sigma^2)), adj = 0)
Dear list,
I need to apply a function on each column of each matrix contained in
a list. Consider the following code,
x - 1:3
my.data - list(matrix(c(1,2,3,4,5,6),ncol=2),
matrix(c(4,5,6,7,8,9),ncol=2))
par(mfrow=c(2,2))
results - sapply(1:length(my.data),
function(ii)
On Fri, 30-May-2008 at 07:29AM +0100, Prof Brian Ripley wrote:
I think what you are asking for is not a tilde. That is a raised symbol
(an accent), and not as in TeX's $\sim$ . It is character 126 in the Adobe
Symbol encoding (Adobe's name is '\similar'), so one way is
I need to apply a function on each column of each matrix contained in
a list. Consider the following code,
x - 1:3
my.data - list(matrix(c(1,2,3,4,5,6),ncol=2),
matrix(c(4,5,6,7,8,9),ncol=2))
par(mfrow=c(2,2))
results - sapply(1:length(my.data),
function(ii)
Hi there everybody,
I would like to be able to modify a string like 'Param_i = Val_i-1' in 'Param_i
= Val_i-2' in a text file
Indeed I have two arrays :
-One for the name of my parameters
-One for the new values for each parameter
I would like to write a loop like this one :
...
SCAN -
Thank you for the suggestions (off-list as well). I think the best
option may eventually be an explicit for loop to make things clearer.
To clarify a bit, I've used the plot function in the example where in
fact it is a numerical integration (which is why I need to pass an
additional
See R News 6/2.
On Thu, May 29, 2008 at 8:04 AM, Birgit Lemcke
[EMAIL PROTECTED] wrote:
Thanks that was a huge help.
Now I am using this:
pdf(InfLengMaxVarHomogeneity.pdf)
plot(inflorescence_length_Max~Sex,xaxt=n #unterdrückt normale x-Achse
, ylab=inflorescence length max, main=Bartletts
Great, thanks a lot, guys!
Only thing is I now have two working versions, that yield *slightly*
different results.
ACCaov - aov(Acc ~ Cond + Error(Subj/Cond), WMU3C)
ACClme - lme(Acc ~ Cond, random = ~1 | Subj/Cond, WMU3C) # what does
'~1 | Subj/Cond' mean?
summary(glht(ACClme,
Thanks for all your replies and sorry for a negligence in my examples. They
are very simplified to reflect very roughly the structure of the case I deal
with, which is too complicated to be quoted here.
What I deal with is:
1) 'theta' which is vector with length 2, being known to me (eg.
Hi
This code should explain what I'm trying to do
strptime(30-Jan-08, %d-%b-%y)
[1] 2008-01-30
format(strptime(30-Jan-08, %d-%b-%y) , %b-%y)
[1] Jan-08
strptime(format(strptime(30-Jan-08, %d-%b-%y) , %b-%y) ,
%b-%y)
[1] NA
I have a string of the form Jan-08 but I can't seem to get
On 30 May 2008, at 11:25, Romain wrote:
...
SCAN - scan(File.txt,sep=\n, what=raw,blank.lines.skip=F)
For (i in 1:Nb_param)
{
sub('Param[i] = Value_i-1','Param[i] = Value_i-2', SCAN)
}
...
I Know how to modify a string with sub when it is a fixed string :
sub((K =)([0-9]*),paste(\\1,
As yearmon in the zoo package can convert a month/year to
an object of class yearmon and that can be converted to other
classes, e.g.
library(zoo)
x - Jan-08
as.Date(as.yearmon(x, %b-%y))
[1] 2008-01-01
On Fri, May 30, 2008 at 7:58 AM, Creighton, Sean
[EMAIL PROTECTED] wrote:
Hi
This code
On May 30, 2008, at 5:37 AM, baptiste Auguié wrote:
Thank you for the suggestions (off-list as well). I think the best
option may eventually be an explicit for loop to make things
clearer. To clarify a bit, I've used the plot function in the
example where in fact it is a numerical
Carlos,
There are many sources of real datasets (in R itself, on the web), you
just need to look a little. For teaching purposes, I think it is
always better to use real datasets than to use simulated ones.
One thing bothers me, though. You imply that in all the examples you
have the
Dear All,
I suspect I'm becoming a bit anoying but... following the question
about the low r2 vs signific regression (thru ANCOVA), the reviewer
suggested (and I quote here):
'If you used RNA/DNA ratios as the dependent variable, but gut content
as the independent variable, and fish
2008/5/29 Frank E Harrell Jr [EMAIL PROTECTED]:
Patrizio Frederic wrote:
dear Harrell,
thank you for quick reply and suggestions. I still have the problem:
library(Design)
x = rnorm(100)
y = runif(100)(exp(x)/(1+exp(x)))
y = 0*y+1*y
d = datadist(x,y)
options(datadist=d)
fit =
On 5/30/2008 8:54 AM, Doran, Harold wrote:
I have a very large tab delimited file (~ 1.97 GB) that I need to read
in to R. The data contain 10 columns and there are millions of rows.
I need all rows of the data, but I only need the first column in the
data. I was looking at the ?read.delim and
Very nice, i really like this one! May confuse non-R users, but
that's not a concern here.
Thanks a lot,
baptiste
On 30 May 2008, at 13:19, Charilaos Skiadas wrote:
I think this is where the beauty of ... comes in, the following
should be doing just what you want:
sapply(my.data,
Thanks,
read.delim('cleanitem.txt', header = TRUE, nrows = 10, colClasses =
rep(c(factor, NULL), c(1,6)))
Is exactly what I needed.
-Original Message-
From: Duncan Murdoch [mailto:[EMAIL PROTECTED]
Sent: Friday, May 30, 2008 9:25 AM
To: Doran, Harold
Cc: R-help@r-project.org
Hi Ullrich,
# what does '~1 | Subj/Cond' mean?
It is equivalent to your aov() error structure [ ... +Error(Subj/Cond) ].
It gives you a set of random intercepts, one for each level of your nesting
structure.
## To get some idea of what's being done it helps to have a continuous
covariate in
Dear All,
I would like to automate the creation of new environments.
What I currently do to create environments is the following:
en1 - new.env(parent = .GlobalEnv)
en2 - new.env(parent = .GlobalEnv)
en3 - new.env(parent = .GlobalEnv)
...
I would like to be able to do something like
Dear R community,
Using nlme library I have developed a nonlinear mixed model. Incorporating
an autoregressive model gives me an error that I can't allocate vector of
size X. The problem is that my computer does not have enough physical memory
most probably due to a large number of observations
Powell, Jeff:
I would like to be able to do something like the following:
nms - c(en1, en2, en3, ...)
for(i in nms) i - new.env(parent = .GlobalEnv)
Try
for(i in nms) assign(i, new.env(parent = .GlobalEnv))
--
Karl Ove Hufthammer
__
Dear R-ers,
My aggregation saga continues.
Using the following sequence, I can calculate any statistic for row
groups and merge the result back to all associated rows ...
WM = by( D60, D60[ KeyProfA], FUN=function(x) weighted.mean( x$IAC, x$Wt))
D60$IAC.WM = as.numeric( WM[ D60$KeyProfA])
Hi Carlos,
Carlos López wrote:
I´m trying to find datasets that will give me residuals, after applying
the lm function, with no normality, non linearity, and heteroscedacity
so I can try to exemplify
those cases in the linear regression model. Can you give any advice on
what datasets would be
Hello,
We have R working with Rmpi/openmpi, but I'm a little worried. Specifically,
(a) the -np flag doesn't seem to override the hostfile (it works fine with
fortran hello world) and (b) I appear to have twice as many processes running
as I think I should.
Rmpi version 0.5.5
Openmpi
Dear Roland and Carlos,
There are some examples in my R and S-PLUS Companion to Applied Regression;
the data sets are in the car package.
I hope this helps,
John
--
John Fox, Professor
Department of Sociology
McMaster University
Hamilton, Ontario, Canada
web:
Hi all
Four questions regarding Unicode.
Three Windows questions. I am using
- a PC with Windows XP (Build 20600.xpsp080413-2111 (Service Pack 3);
- the following R version:
R.version
platform i386-pc-mingw32
arch i386
os mingw32
system i386, mingw32
status
Dear R-help.
When plotting 3D wireframe plots with the lattice package, there is often
a lot of space between the actual plot (i.e., cube) and panel borders. For
a single wireframe, this is not a problem, but when plotting multiple
wireframes, the wasted space decreases the size of the plots
Yes! Again, thank you ALL very, very much.
Even simply lurking on the list generates many gems worth collecting.
DaveT.
From: [EMAIL PROTECTED] [mailto:[EMAIL PROTECTED]
I totally agree both of you. This is a super place to mature the R.
I learn a lot from this R heaven!
Chunhao
Quoting Esmail
I'm looking for something in R to fit a multivariate linear model
robustly, using
an M-estimator or any of the myriad of other robust methods for linear
models
implemented in robustbase or methods based on MCD or MVE covariance
estimation (package rrcov).
E.g., one can fit an mlm for the iris
I was trying to plot some data in R. I used the following code to draw a loess
fit and got the output as
?lines(lowess(log(abs(t(res))), log(abs(t(synthesised, col=red)
Error in lowess(log(abs(t(res))), log(abs(t(synthesised :?? NA/NaN/Inf in
foreign function call (arg 1)
Then I
Dear R community,
Using nlme library I have developed a nonlinear mixed model. Incorporating an
autoregressive model gives me an error that I cant allocate vector of size X.
The problem is that my computer does not have enough physical memory most
probably due to a large number of
On 5/30/08, Karl Ove Hufthammer [EMAIL PROTECTED] wrote:
Dear R-help.
When plotting 3D wireframe plots with the lattice package, there is often
a lot of space between the actual plot (i.e., cube) and panel borders. For
a single wireframe, this is not a problem, but when plotting multiple
Patrizio Frederic wrote:
2008/5/29 Frank E Harrell Jr [EMAIL PROTECTED]:
Patrizio Frederic wrote:
dear Harrell,
thank you for quick reply and suggestions. I still have the problem:
library(Design)
x = rnorm(100)
y = runif(100)(exp(x)/(1+exp(x)))
y = 0*y+1*y
d = datadist(x,y)
Hi,
to put it simply. Windows cannot handle utf-8 data. There is no utf-8
locale available.
If your corpus only contains Russian data, maybe English glosses etc.
you can try to set lang of Rgui.exe to Russian.
Then at least you can use grep, strsplit because they are depending
on the
I was trying to plot some data in R. I used the following code to draw a loess
fit and got the output as
?lines(lowess(log(abs(t(res))), log(abs(t(synthesised, col=red)
Error in lowess(log(abs(t(res))), log(abs(t(synthesised :?? NA/NaN/Inf in
foreign function call (arg 1)
Then I thought
On Fri, May 30, 2008 at 12:04 PM, threshold [EMAIL PROTECTED] wrote:
Thanks for all your replies and sorry for a negligence in my examples. They
are very simplified to reflect very roughly the structure of the case I deal
with, which is too complicated to be quoted here.
What I deal with is:
On 5/30/2008 12:58 PM, Hans-Jörg Bibiko wrote:
Hi,
to put it simply. Windows cannot handle utf-8 data. There is no utf-8
locale available.
Code page 65001 is utf-8. Most text editors (including Notepad) include
an option to save in the UTF-8 encoding.
Some programs don't fully support
1. Please don't double post.
2. You probably have some synthesized values that are 0 and so get -Inf
for their logs, which is NaN, thereby producing the error messages.
Failing that, your data may be messed up in some way (character values,
stray non-characters, ...)
-- Bert Gunter
I just wonder that if you know theta and x. Why can't you specify
these in your functions. such as
f1-function(theta1){1+2+theta1[1]} (Let me assume theta[1]=1, theta[2]=2)
f2-function(theta1){f1(theta1)*exp(3)*theta1[2]} (again, I assume x=3)
f3-function(theta1){f1(theta1)-f2(theta1)}
You could do something like this:
mydata=c(1,2,1,1,6,7,-1,-1,5,-1)
color= as.numeric(mydata== -1) +1
plot(mydata,col=color)
This will give you a plot where the -1's are in red (color = 2) and the
other numbers in black (color=1).
Julian
uv wrote:
Hi. I am plotting graphs for values
Hello all,
I used pushViewport to create a viewport with a 9*8 layout. In each cell, I
used lattice to draw plots. In the final output, if I used the
trellis.device(),i.e., the default device, the 9*8 plots I intend to draw
seem ok in terms of the letter size and everything in each cell. But if
I have a general question about how to interpret the plotcp() graph.
The cross-validation xerror value typically follows a decreasing pattern,
from approximately 1.0 at the root node, then it crosses the 1SE boundary,
reaches a plateau, and decreases further when the tree gets very complex
What happens if plotcp() shows an V-shape profile?
It seems that the predictors (I have a few) do not help the splits.
Would it be reasonable to prune it at size 6?
Or perhaps rpart() is not suitable for this analysis?
What I meant on the last sentence was perhaps I made a
I noticed that the ksmooth documentation mentions the function is real slow and
that better smoothing functions are available. Does this mean that the method
itself is a poor method for smoothing data or just a slow function? Thanks.
[[alternative HTML version deleted]]
I am extremely puzzled by this behavior in R. I have a data frame called
Trials in which I have results from an experiment. I am trying to do a
subjects analysis, but getting weird results. Each row has 1 trial in it,
which includes a column for the subject number I get the list of subject
Quoting Duncan Murdoch [EMAIL PROTECTED]:
On 5/30/2008 12:58 PM, Hans-Jörg Bibiko wrote:
to put it simply. Windows cannot handle utf-8 data. There is no
utf-8 locale available.
Code page 65001 is utf-8. Most text editors (including Notepad)
include an option to save in the UTF-8
On 5/30/2008 4:12 PM, Hans-Joerg Bibiko wrote:
Quoting Duncan Murdoch [EMAIL PROTECTED]:
On 5/30/2008 12:58 PM, Hans-Jörg Bibiko wrote:
to put it simply. Windows cannot handle utf-8 data. There is no
utf-8 locale available.
Code page 65001 is utf-8. Most text editors (including Notepad)
Hi,
I often run predictive models on large datasets with multiple
combination of parameter space. I am wondering if there is any way to
quickly check the execution plan, like how much time does it take to
run a model. Here is a more specific example:
I have a set of datasets, S, of size
Thanks for the solution. But I wish someone would tell us why the solution
makes sense. If you have changed your R environment in any way from the
defaults please let us know.
Frank
--
Frank E Harrell Jr Professor and Chair School of Medicine
Department
Mark,
Thanks for the reply.
hi: subset doesn't know what i is but I don't have enough knowledge
about scope to know what it's actually doing in that case. my
point is that i wouldn't put i inside a subset command and expect it to
know the value. scope is very complex in R so
doing
Patrizio Frederic wrote:
Thanks for the solution. But I wish someone would tell us why the solution
makes sense. If you have changed your R environment in any way from the
defaults please let us know.
Frank
--
Frank E Harrell Jr Professor and Chair School of Medicine
On Fri, 30 May 2008, Patrizio Frederic wrote:
Thanks for the solution. But I wish someone would tell us why the solution
makes sense. If you have changed your R environment in any way from the
defaults please let us know.
Frank
--
Frank E Harrell Jr Professor and Chair School of
On 30/05/2008 5:50 PM, Robert Felty wrote:
Mark,
Thanks for the reply.
hi: subset doesn't know what i is but I don't have enough knowledge
about scope to know what it's actually doing in that case. my
point is that i wouldn't put i inside a subset command and expect it to
know the value.
On Fri, 30 May 2008, Duncan Murdoch wrote:
On 5/30/2008 4:12 PM, Hans-Joerg Bibiko wrote:
Quoting Duncan Murdoch [EMAIL PROTECTED]:
On 5/30/2008 12:58 PM, Hans-Jörg Bibiko wrote:
to put it simply. Windows cannot handle utf-8 data. There is no utf-8
locale available.
Code page 65001 is
Prof Brian Ripley wrote:
On Fri, 30 May 2008, Patrizio Frederic wrote:
Thanks for the solution. But I wish someone would tell us why the
solution
makes sense. If you have changed your R environment in any way from the
defaults please let us know.
Frank
--
Frank E Harrell Jr Professor
Duncan, Mark, et al.,
It's a scoping problem. To be helpful, subset looks up variables within
the dataframe first. Since it can find Subj there, that's the one it
uses in the test.
The easiest solution is just to make sure that Subj is named something
that isn't a column of Trials, e.g.
Hi,
I would like to get all X iterations in optim output in matrix form.
I know about the follow approach:
sink(reportOptim)
optim( ..., control=list( trace=6,..) )
sink()
all_iterOptim - readLines(reportOptim)
unlink(reportOptim)
all_iterOptim - all_iterOptim[ grep( '^X',
Hi R Gurus:
I have used the following in the past:
g77 -O2 -c test.f
g77 -shared -o test.dll test.f
with success.
Now I'm getting
ld: dllcrt2.o no such file
I have the latest version of Rtools28.exe
Any suggestions would be much appreciated.
Sincerely,
Edna Bell
Hi list,
I have a question about using *read.table()* to read in a txt file.
Basically, it
consists of 16346 rows, 6 columns (no header). The code I used is:
exprSet - read.table('process_all4_GSA2.txt', row.names = 1,header =FALSE)
and I got an error message:
exprSet -
Hi Steve,
You can use write.table:
write.table(x, file=/Users/Desktop/Data.txt,
sep=,row.names=F,col.names=F)
Cheers,
Julian
Stropharia wrote:
Dear R users,
I've had no joy finding a solution to this online or in any of my R books.
Many thanks in advance for any help you can give.
I'm
On May 30, 2008, at 7:56 PM, ss wrote:
and I got an error message:
exprSet - read.table('process_all4_GSA2.txt', row.names = 1,header
=FALSE)
Error in read.table(process_all4_GSA2.txt, row.names = 1, header
= FALSE)
:
duplicate 'row.names' are not allowed
I would say that's pretty
On Fri, 30 May 2008, Frank E Harrell Jr wrote:
Prof Brian Ripley wrote:
On Fri, 30 May 2008, Patrizio Frederic wrote:
Thanks for the solution. But I wish someone would tell us why the
solution
makes sense. If you have changed your R environment in any way from the
defaults please let us
Rtools28 does not contain g77 to my knowledge (that is, it did not in the
version I last downloaded). But it does contain dllcrt2.o.
The Fortran compiler is gfortran. I think you have another GCC toolset in
your path.
On Fri, 30 May 2008, Edna Bell wrote:
Hi R Gurus:
I have used the
72 matches
Mail list logo
