HI,
Possibly R FAQ: 7.31
data1New-data1
data1New$A- round(data1New$A,2)
data2New- data2
data2New$A- round(data2New$A,2)
merge(data1New,data2New,by=A)
# A B C
#1 0.0 0.9 10.0
#2 1.1 0.6 11.1
#3 1.4 0.7 11.4
#4 3.1 0.4 13.1
#5 4.4 0.8 14.4
A.K.
Hello,
Lets say we have these two
Greetings,
I have just downloaded R onto a 64bit PC running Microsoft 7 Home Edition
via Rgui. I have quite a bit of programming experience, though not as a
professional programmer. I am a Measurement and Assessment professional
(standardized testing). I would like to be able to write R
Hello Guys,
I am trying to build R using the Intel compilers: ICC, ICPC. The first time I
ran the built I got the error below:
mlutils.c(107): error: floating-point operation result is out of range
double NA_REAL = ML_NAN;
^
mlutils.c(108): error: floating-point operation
On Tue, 28 May 2013, Marcus Nunes marcus.nu...@gmail.com writes:
Hello all
I'm trying to use simulated annealing to optimize a function. I know I can
use ?optim with method=SANN to do it.
However, I'd like to make optim to search for the best solution among a set
of possible points in my
hallo,i have a problem with running some code in R2Winbugs am getting the
following error.Error in FUN(X[[2L]], ...)
What might be the problem.with kind regards.Alex
[[alternative HTML version deleted]]
__
R-help@r-project.org mailing list
Hallo.
Is there any package / code snippet to test the distribution
assumption, heteroskedasticity, omitted variables, and linearity with
the models estimated by maximum likelyhood? I especially need it for
three type of models:
* binary choice (probit and probit with non-normal distribution) --
Hi,
I do not use R for Windows. But I would say that you have to run
'Rscript.exe' in a CMD prompt, if I am not mistaken. Not in 'Rgui'.
In 'Rgui', use 'source'.
Hope this helps,
Pascal
On 05/29/2013 03:07 AM, Mark Russell wrote:
Greetings,
I have just downloaded R onto a 64bit PC
On May 28, 2013, at 19:26 , Mark Breman wrote:
Thanks everyone for the suggestions.
If the problem is caused by a restriction in one of the (software)
components i'm using (OS, terminal etc), should the remaining of the pasted
code simply be missing after the point where the error occurs?
Hello,
I would say:
1. Which R version?
2. Which OS?
3. Where is the requested commented, minimal, self-contained,
reproducible code?
Regards,
Pascal
On 05/29/2013 10:05 AM, alex mutuku wrote:
hallo,i have a problem with running some code in R2Winbugs am getting the
following error.Error
On May 29, 2013, at 08:48 , Pascal Oettli wrote:
Hello,
I would say:
1. Which R version?
2. Which OS?
3. Where is the requested commented, minimal, self-contained, reproducible
code?
4. Use a sensible Subject: line, instead of contributing to the (no subject)
thread,
which on my
On 28-05-2013, at 20:07, Mark Russell gibsons...@cox.net wrote:
Greetings,
I have just downloaded R onto a 64bit PC running Microsoft 7 Home Edition
via Rgui. I have quite a bit of programming experience, though not as a
professional programmer. I am a Measurement and Assessment
Dear all,
I have a probably silly question on how to specify model for a
split-plot design and for a within-subject design (1way). I have two
factors: factor A with 5 levels and factor B with 2 levels, for each
level of factor A I have 4 subjects that experiment both levels of
factor B, so this
Hallo.
Is there any package / code snippet to estimate quantile regression
for a binary choice model (like probit) and selection model (like
heckit)? I found that quantreg package can estimate tobit-like model,
but I can't figure out how to tweak it for probit / heckit.
Best wishes,
Michal
Hi,
I would like to estimate a VAR(5) model, but including lags T-1, T-7, T-14,
T-21 and T-28 instead of the usual T-1, T-2, T-3, T-4, T-5.
But it seems I cannot accomplish this by using the below function.
VAR(y, p = 1, type = c(const, trend, both, none),
season = NULL, exogen = NULL, lag.max =
Dear all,
I have been experiencing difficulties in running the x12 package. I'm
encountering the following error message which may or may not be R-related:
CMD.EXE was started with the above path as the current directory.
UNC paths are not supported. Defaulting to Windows directory.
I don't
Hi,
It is a Windows problem.
CMD.exe does not support having a UNC path as the current directory.
Regards,
Pascal
On 05/29/2013 07:39 PM, Jose Iparraguirre wrote:
Dear all,
I have been experiencing difficulties in running the x12 package. I'm
encountering the following error message which
On 05/29/2013 02:02 AM, jcrosbie wrote:
Hi, I'm trying to download data from:
http://www.ngx.com/settlehistory.html
Is it possible to fetch the data with R?
Hi jcrosbie,
The simplest way seems to be to highlight the desired spreadsheet (less
the title row), copy (Ctrl-C) and paste (Ctrl-V)
Hi,
The combination read.table (and its arguments) + stdin also can be
used, directly in R.
read.table(stdin(), ...)
Regards,
Pascal
On 29/05/2013 21:35, Jim Lemon wrote:
On 05/29/2013 02:02 AM, jcrosbie wrote:
Hi, I'm trying to download data from:
http://www.ngx.com/settlehistory.html
I tried reading in the data using the XML package, but I can't figure out
how to read either ALL of the tables or a particular table. The code below
just reads the first table. Maybe someone else will know how.
Jean
library(XML)
look - readHTMLTable(http://www.ngx.com/settlehistory.html;)
Does anybody know whether R has a package for supermodular function
optimization?
Regards
Arjan Ruijs
[[alternative HTML version deleted]]
__
R-help@r-project.org mailing list
https://stat.ethz.ch/mailman/listinfo/r-help
PLEASE do read the
Dear all
I try to assign class attributes to strings. Depending on the approach I use,
it works (including method dispatch) or fails.
Let me clarify with an example.
Let:
AONtptmodelist - c(FOOT,BICY)
When I assign class attributes directly to these strings, it fails:
class(BICY) -
Hi,
You might need to check library(boot). I have never used that before. So, I
can't comment much. It is better to post on R-help list. I had seen your
postings on Nabble in the past. Unfortunately those postings were not accepted
in R-help. You have to directly post at
You probably want from windows GUI:
source(test.R)
To read help to learn about this command type:
?source
at the command prompt. (Similar pattern get help on other commands too -
its an important R skill/habit).
If your are going to do a lot of script development, I would highly
recommend
Second the RStudio but also suggest Tinn-R for a beginner as the more extenseve
code highlighting can be very useful.
John Kane
Kingston ON Canada
-Original Message-
From: rb...@atsu.edu
Sent: Wed, 29 May 2013 08:17:47 -0500
To: gibsons...@cox.net
Subject: Re: [R] Rank Amateur
No attachment. The R-help list tendst to strip out many type of attached files
though pdf and txt , among others get through.
It is better to supply the example in the email itself if possible. Have a look
at https://github.com/hadley/devtools/wiki/Reproducibility for suggestions.
John Kane
Inline:
On Tue, May 28, 2013 at 9:22 PM, Berend Hasselman b...@xs4all.nl wrote:
On 28-05-2013, at 19:31, bcrombie bcrom...@utk.edu wrote:
How do I assign a variable to R output that is generated by default
functions such as:
stat.desc(mydata, basic=TRUE, desc=TRUE, norm=FALSE, p=0.90)
library(boot) is wholly unnecessary.
Hint:
1000 * 5 = 5000
?sample
-- Bert
On Wed, May 29, 2013 at 6:15 AM, arun smartpink...@yahoo.com wrote:
Hi,
You might need to check library(boot). I have never used that before. So, I
can't comment much. It is better to post on R-help list. I had
There are two separate issues that seem to be unclear.
Concerning the class assignment: You cannot assign a class to a
character string. You can assign a class to an object that contains a
character string:
a - b
# ok
class(a) - AONmode
# not ok
class(b) - AONmode
Error in class(b) - AONmode
AIC is a different story. To do hypothesis tests on terms, use anova()
or dropterm() (as done in the book):
library(MASS)
example(polr)
dropterm(house.plr, test = Chisq)
Single term deletions
Model:
Sat ~ Infl + Type + Cont
DfAIC LRT Pr(Chi)
none3495.1
Infl2 3599.4
This is really odd, and probably 100% local to me, but I'm at a loss as to a
next step.
After narrowing things down, here's how to reproduce:
library(rgl)
showSomething - function() {
open3d()
points3d(rnorm(10),rnorm(10),rnorm(10))
axes3d()
This is a bit like asking how should I tweak my sailboat so I can explore the
ocean floor.
url:www.econ.uiuc.edu/~rogerRoger Koenker
emailrkoen...@uiuc.eduDepartment of Economics
vox: 217-333-4558University of Illinois
fax: 217-244-6678
a) You cannot assign attributes to literal values such as BICY. You must
assign them to variables, as in
x - BICY
class(x) - AONmode
b) You cannot give the separate elements of a vector their own separate
attributes unless the vector is of mode list. Ordinarily, attributes apply to
the whole
Dear R users.
We are currently developing a R package (GTrendsR) that allows to retrieve
data from google trends. To do so, Im using the RCurl library. At this
point everything works perfectly (i.e. the data obtained from R is identical
to the data obtained directly from the web site).
Dear all,
I want to make a baysian optimisation of the parameters values of a model,
written in FORTRAN. I thus need a parallel modeling procedure.
I started to use the R package snow but I went through trouble and I
hope that you will bring me light on this problem.
Briefly, I first make a
Sorry forthat and thanks for the hint.
At marked the appended text as rawtext, maybe that is the reason why
it was cut.
Here we go again.
##
x=runif(100, 1, 2)
y=runif(100, 2, 4)
z=runif(100, 1, 4)
### 2d example
Thank you, I will have to wait until I get home form work to test the XML. But
I'm looking to do something more along the lines of an automated process over
copy and pasting.
James
From: Adams, Jean jvad...@usgs.gov
To: jcrosbie ja...@crosb.ie
Cc: R help
Hi,
Don't know if this helps:
set.seed(244)
data1- generate(10)
data1
x1 x2 x3 x4 x5 x6
1 -0.89558948 -0.5609722069 -1.0782688 -1.2461548 -0.58857050 1.5708187
2 -0.82379187 -0.9056306209 -1.5700755 -0.6045090 -0.19683863 -0.8969354
3
This is most relevant to R-sig-mac. There are two different rgl devices
on OS X, depending how you are running this. One based on X11 and one
on Apple's GL. In particular, are you using command-line R or R.app?
I seem to be able to reproduce it using R.app, in which case it is most
Thank you Brian. Yes, problem is on R.app I will send to R-sig-mac. Thanks,
Bryan
On May 29, 2013, at 11:07 AM, Prof Brian Ripley rip...@stats.ox.ac.uk wrote:
This is most relevant to R-sig-mac. There are two different rgl devices on
OS X, depending how you are running this. One based
Dear R-users,
in case of linear model,
lm(Y ~ X)
is equivalent to
glm(Y ~ X, family=gaussian(link=identity))
In case of the exponential model
lm(log(Y) ~ X)
why
glm(Y ~ X, family=gaussian(link=log))
is not equivalent? Is there an equivalent glm to lm(log(Y) ~ X)?
Thank you for your help
Stefano
Dear Nello
Thanks for the clarification.
It was of course an elementary mistake to assume that the for loop would modify
the original AONtptmodelist.
I also understand now that the issues was that I cannot attribute a class to a
literal.
Your proposed alternative seems to work now, including the
Hi,
I have to do some Radon modelling and I have categorical and non
categorical datasets. I have been considering Artificial Neural Networks to
do this. I was wondering has anybody done anything like this before and
have you any advice before I start and where there might be some good
tutorials
On Wed, 29 May 2013, Shane Carey wrote:
I have to do some Radon modelling and I have categorical and non
categorical datasets. I have been considering Artificial Neural Networks
to do this. I was wondering has anybody done anything like this before and
have you any advice before I start and
On Wed, May 29, 2013 at 7:57 AM, Stefan Lüdtke slued...@gfz-potsdam.dewrote:
x=runif(100, 1, 2)
y=runif(100, 2, 4)
z=runif(100, 1, 4)
data_xyz=as.data.frame(cbind(x, y, z, a=rep(c(1:10), 10), b=rep(c(1:2),
each=50)))
custom.panel = function(x, y, z, subscripts, ...)
{
R-helpers:
I'm trying to retrieve the contents of a directory from an ftp site
(ideally, the file/folder names as a character vector):
ftp://e4ftl01.cr.usgs.gov/MOTA/MCD12C1.005/;
# (MODIS data)
Where I get the following error via RCurl:
require(RCurl)
url -
On May 29, 2013, at 17:23 , Stefano Sofia wrote:
Dear R-users,
in case of linear model,
lm(Y ~ X)
is equivalent to
glm(Y ~ X, family=gaussian(link=identity))
In case of the exponential model
lm(log(Y) ~ X)
why
glm(Y ~ X, family=gaussian(link=log))
is not equivalent?
Y is assumed
Dear R users,
There has been discussion on how to collect R usage statistics before and the
discussion has led to some good results such as
http://neolab.stat.ucla.edu/cranstats/ and http://crantastic.org/. Recently, I
tried to put together a package that allows R users to rate, comment, and
data1$A[8]
#[1] 1.4
data2$A[15]
#[1] 1.4
data2$A[15]==data1$A[8]
#[1] FALSE
You can check these links:
http://r.789695.n4.nabble.com/Comparing-decimal-numbers-td3251437.html
http://rwiki.sciviews.org/doku.php?id=misc:r_accuracy
A.K.
The FAQ says Other numbers have to be rounded to (typically)
Hi,
I have a problem that appears to be quite common - warnings that 'x' must
be numeric. How can I get rid of that problem?
specifically, I am trying to run an mrpp with a species-site data set plus
a second file containing factors. Both files are saved as .csv and comma
separated. I have tried
a) I've recently come across the global Goodness-of-ï¬t tests for complex
survey logistic regression. Has it been implemented in R?
Paper
http://med.stanford.edu/medicine/ArcherLemeshowHosmer.pdf
Implementation in Stata
R list members,
I am trying to replicate a command script in STATA and find
equivalency in R to the following Stata command:
egen stratida = concat(year byregn2 group) if sample==1
egen stratidb = concat(year byregn2 group) if sample==2
egen stratidc = concat(year byregn2 group) if sample==3
Dear Jonathan,in the MODIS package I use the following function to list files
within a http or a ftp folder. I'm not really practical in XML stuff but I got
it working somehow.
The LP DAAC has changed from FTP to HTTP, I'm not sure if it is good idea to
use the ftp protocol anymore. I also had
That looks good, thanks a lot!
Iwill give it a try in my real problemsoon.
Cheers,
stefan
On 05/29/2013 06:02 PM, ilai wrote:
On Wed, May 29, 2013 at 7:57 AM, Stefan Lüdtke slued...@gfz-potsdam.dewrote:
x=runif(100, 1, 2)
y=runif(100, 2, 4)
z=runif(100, 1, 4)
Dear users,
I am new to spatial data analysis in R. I have downloaded
species occurrence point from the gbif database but when i'm trying to give
coordinates to these points an error message comes up saying that the rows in
the dataframe are not the same with spatial points. Can
I am getting the error Error in names(data)[ix] : only 0's may be mixed with
negative subscripts and cannot quite figure out how to resolve this. Any
ideas?
attempt - read.csv(attempt.csv)
str(attempt)
'data.frame': 9424 obs. of 9 variables:
$ cgPlaId : int 647 647 647 647 647 647 647
What do you get when you type in these commands?
lapply(mzbtaxa, class)
lapply(mzbsites, class)
Do you get the same error if you use the as.is=TRUE argument in the
read.csv() function?
mzbtaxa-read.csv(MZBTaxa.csv, header=TRUE, dec=,, sep=;, as.is=TRUE)
mzbsites-read.csv(MZBSites.csv,
You are limiting yourself to getting help from only those people that
understand the STATA language. You may improve your chances for a helpful
response if you tell us in English what the STATA commands mean.
Jean
On Wed, May 29, 2013 at 10:37 AM, Daniel Tucker dtuck...@u.rochester.eduwrote:
Hi,
I am not sure if this is what you wanted.
set.seed(125)
simdata- replicate(S,generate(500))
set.seed(125)
data1- generate(500)
data2-
replicate(S,as.data.frame(t(sample(as.data.frame(t(data1)),100,replace=TRUE
dim(data2)
#[1] 12 1000
data2[1:5,1:5]
# [,1] [,2]
Hey all!
I have a time series dataset like this:
DateTime Var
112
1 14
1 1 5
1 2 8
1 2 8
1 2 9
213
21 4
214
I created
May be this helps:
mymatrix- matrix (1:5, nrow=5, ncol=columnas, byrow=TRUE,
dimnames=(list(myrownames, c(1,2,3,4,5
#Error in matrix(1:5, nrow = 5, ncol = columnas, byrow = TRUE, dimnames =
(list(myrownames, :
# object 'columnas' not found
myrownames[which.max(nchar(myrownames))]
#[1]
Hi,
May be I misunderstood your question:
dat- read.table(text=
Date Time Var
1 1 2
1 1 4
1 1 5
1 2 8
1 2 8
1 2 9
2 1 3
2 1 4
2 1 4
The actual date col is in -mm-dd format, and when I apply this code
to my actual data, it mess up the order
On Wed, May 29, 2013 at 11:37 AM, arun smartpink...@yahoo.com wrote:
Hi,
May be I misunderstood your question:
dat- read.table(text=
DateTime Var
112
Hi,
You're most likely to get an informative answer if you tell us exactly what
you did. A reproducible example was best. You might also be better off
posting to the r-sig-geo list rather than to R-help.
Sarah
On Wednesday, May 29, 2013, jickngea alexand wrote:
Dear users,
I am
dat1- read.table(text=
Date Time Var
11/1/2012 1 3
11/1/2012 1 1
11/1/2012 1 1
11/1/2012 2 3
11/1/2012 2 1
11/1/2012 2 1
11/2/2012 1 1
11/2/2012 1 3
11/2/2012 1 1
Interesting idea but would you be willing to change the name. Twitter,
Stackoverflow, and the like use rstats to refer to R itself and confusion seems
probable.
https://mobile.twitter.com/search?q=%23rstats
Michael
On May 29, 2013, at 19:21, Johnny Zhang johnny...@yahoo.com wrote:
Dear R
On Thu, May 30, 2013 at 3:29 AM, Marco Pomati pt...@bristol.ac.uk wrote:
a) I've recently come across the global Goodness-of-ï¬t tests for complex
survey logistic regression. Has it been implemented in R?
It's quite hard to definitively say that something hasn't been implemented
in R, but I
Hello,
I am trying to run glmm to test the effect of the three fixed effects [AGE
(weaned vs. unweaned claf), LOCATION (zoo vs. park), MOTher's social status
(matriarch vs. nonmatriarch)] and one random effect [ID (12 different
calves of whom I have multiple but unbalanced observations)] on the a
Hello,
I have a whole bunch of data to two decimal places. I've done some
arithmetic with them, so floating point becomes an issue.
x - c(1, 0.15,(0.1+.05),0.4)
I want to do something like this:
x[x==0.15]
But you'll notice that is troublesome with the well known floating point
issue. So
Let me explain what I am planning on doing; some background is probably
necessary and in hindsight the first email doesn't accurately describe what
I am looking to do. I'm trying to run a conditional logit model, but I want
to have a different strata depending on the value of another variable
if
Hi Stefan
Not a perfect solution but may get you some of the way.
I can't spend any more time on it
I think your groups argument is not valid so made a groups argument
Try this
data_xyz$gp - sapply(data_xyz$a == 1, pmatch, c(TRUE,FALSE))
cloud(z ~ x*y|b, data=data_xyz,
groups= gp,
On May 29, 2013, at 3:10 PM, Daniel Tucker wrote:
Let me explain what I am planning on doing; some background is
probably
necessary and in hindsight the first email doesn't accurately
describe what
I am looking to do. I'm trying to run a conditional logit model, but
I want
to have a
On May 29, 2013, at 6:27 PM, Peter Lomas wrote:
Hello,
I have a whole bunch of data to two decimal places. I've done some
arithmetic with them, so floating point becomes an issue.
x - c(1, 0.15,(0.1+.05),0.4)
I want to do something like this:
x[x==0.15]
But you'll notice that is
72 matches
Mail list logo
