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 foll
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 kn
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 e
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
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 <- read.table('process_
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,
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', al
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.
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 and
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
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. sc
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
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
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 t
> 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
> Depa
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 4X700
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
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 encodin
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
numb
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]]
__
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 poor
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 [e.
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
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 ranging
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)}
optim
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
-Orig
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 ut
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 wit
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 thoug
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 locale
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)
options
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 mult
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 observatio
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 thou
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 d
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 Es
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 very
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 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: socserv
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 version
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
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$KeyProf
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 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
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
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
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
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, apply
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 a
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)
>> optio
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 am trying to see if it is
possible to tell thi
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 leng
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 data
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 integrat
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:
> H
On Fri, May 30, 2008 at 12:58:27PM +0100, Creighton, Sean wrote:
> 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-Ja
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
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"
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.
theta=
Dear All (particularly Peter, Steve, Bert, Ben F, Ben B, Christian and Gerrit)
After clicking the "send msg" button of my previous msg I felt that
probably it's better to leave the "citing issue" out of it all!
Nonetheless, what I feel now (as a - fish - biologist aspiring to be a
reasonable u
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",
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, linfct=mcp(Co
I found "The R Book" by Michael J. Crawley quite satisfying for the
purpose of learning R.
On May 28, 10:25 pm, "Neil Gupta" <[EMAIL PROTECTED]> wrote:
> Hi I am still fairly new to R but picking it up quickly. I have some
> problems manipulating data in tables. I was wondering if anyone new any g
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 variab
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 <- scan("
> 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),
> >
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
>
> expression(X*s
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 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, sig
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
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 subs
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"
> Received: Friday, 30 May, 2008, 5:13 PM
> Hello!
>
try this:
dat <- read.table(textConnection(
"Id Year
123456 2001
123456 2002
123456 2003
655432 2001
655432 2002
655432 2003"), header = TRUE)
closeAllConnections()
dat[!duplicated(dat$Id), ]
# or
dat[tapply(row.names(dat), dat$Id, head, 1), ]
I hope it helps.
Best,
Dimitris
Dim
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 adju
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 fi
75 matches
Mail list logo
