Dear Jan,
Thank you for your answers. They are very useful. I will try the LaF
package.
Cheers,
--
View this message in context:
http://r.789695.n4.nabble.com/Reading-big-files-in-chunks-ff-package-tp4502070p4537857.html
Sent from the R help mailing list archive at Nabble.com.
Hi,
I intend to use the BiodiversityR package to perfom community analyses. I
will do ordination analyses and diversity comparisons (e.g. accumulation
curves, etc.).
I have repeated measures of plots and want to control for it when performing
analyses. I have not decided yet which analysis I
Yes, I just tried the length() function and it worked beautifully.
Thank you.
On 4/5/2012 6:01 PM, William Dunlap wrote:
I think you are looking for the function called length(). I cannot recreate
your output, since I don't know what is in NZ_Conifers, but with the built-in
dataset mtcars I
The CPU time and elapsed time are essentially identical. (That is, the
system time is negligible.)
Using Rprof, I just ran the code twice. The first time, while
randomForest is doing its thing, there are 850 consecutive lines which
read:
.C randomForest.default randomForest randomForest.formula
Hello,
the R Installation and Administration handbook states in Section C 5.1:
For ‘amd64’ the builds have failed to complete in several different
ways, currently with relocation errors for libRblas.so.
To fix it: add '-shared' to the SHLIB_LDFLAGS, SHLIB_CXXLDFLAGS and
SHLIB_FCLDFLAGS
Hi
I have a dataset with a date variable.
There may be more than one Date Range. i've to create 1 for the dateRange
and 0 for the other.
Please can any one help in this
-
Thanks in Advance
Arun
--
View this message in context:
I completely agree that the development of full-information maximum
likelihood (FIML) estimation for use in packages like lm, lme, lmer, etc.
would make R much more attractive. FIML is better than other approaches to
missing data (e.g., multiple imputation; see Graham, Olchowski, Gilreath,
Hi,
The compute.es package is great for calculating estimated effect size.
However, it is limited to tests that have independent groups (i.e.,
independent samples t-test, ANOVA). Are there any other packages available
for download that can handle estimating effect sizes for repeated-measures
Hello,
Oops!
What happened to the function 'f'?
Forgot to copy and pasted only the rest, now complete.
f - function(x){
nr - nrow(x)
result - matrix(0, nrow=nr, ncol=ncol(x))
colnames(result) - colnames(x)
inp.ord - order(x)[1:nr] - 1 # Keep only one per row,
Hello,
I maybe missing something but this seems like an indexing problem
which doesn't require a loop at all.
Yes, but with 'order'.
# Original example
input - as.matrix(data.frame(a=c(5,1,3,7), b=c(2,6,4,8)))
(input)
desired.result - as.matrix(data.frame(a=c(100,0,100,0),
i am trying to replicate the following graph using xyplot :
attach(x)
plot ( jitter(type), mortality, pch=16, xlim = c(0.25, 3.75))
lines ( c(1-0.375,1.375) , c ( median(mortality[type==1]),
median(mortality[type==1])), lwd=5,col=2)
lines ( c(2-0.375,2.375) , c ( median(mortality[type==2]),
Dear All,
I would like to make a quadratic with a plateau model in R. Is there a
package in R doing this? The bentcableAR package seems won't work.
The link below describes what I am looking for in R exactly:
newbie to R, less than a week, and I ordered some books about R, but I learn
better by examples.. and thus far I cant find a good example of what I am
trying to do... which follows:
assuming one is using any instrument intra-day data... I want to..
open a file (lets name it signal) that will
On 06/04/2012 19:38, Michael Figiel wrote:
Hello,
the R Installation and Administration handbook states in Section C 5.1:
For ‘amd64’ the builds have failed to complete in several different
ways, currently with relocation errors for libRblas.so.
To fix it: add '-shared' to the SHLIB_LDFLAGS,
On Fri, Apr 06, 2012 at 10:25:15AM -0700, ens wrote:
a-sample(1:6,100,replace=T)
a
[1] 2 4 3 4 5 1 3 2 4 3 6 6 2 6 2 1 5 5 3 4 6 1 6 6 3 4 6 6 4 4 5 4 6 5 6
3 4 5 6 3 4 1 6 6 6 4 2 1 1 3 1 5 3 2 2 6 2 5
[59] 2 6 1 6 1 1 6 4 4 2 2 3 4 5 6 1 6 4 6 1 5 1 1 2 1 3 4 4 6 3 1 4 1 1 1
5 5 2 4 6
Dear list!
I have problem with buffer size (width) in package rgeos. I would like to
expand given geometry (points) to specified width based on the z value from
attribute table. Here is example:
point - data.frame(x=c(10,20), y=c(10, 10), z = c(2,7))
point_shp -
Hi there,
I want to estimate simultaneous equation model with panel data. The model looks
as follows
Y1=a0+a1*X1+a2*X2
Y2=b0+b1*X2+b2*X1
X1=Z1-(Y1/a1)
X2=Z2-(Y2/b1)
I
In this model Y1, Y2, X1 and X2 are endogenous variables; Z1, Z2 are exogenous
variables and a0, a1, a2, b0, b1 and b2 are
Hi Petr,
Thank you. Yes, normalized.matrix and normalized are the same, my mistakes.
Operators: like Cyclic loess normalization and quantile Normalization,
from limma or affy packages.
Would cv= sd/mean be a better replacement for sd?
On 4/5/12, Petr Savicky [via R]
On 04/ 7/12 08:13 AM, Prof Brian Ripley wrote:
On 06/04/2012 19:38, Michael Figiel wrote:
Hello,
the R Installation and Administration handbook states in Section C
5.1:
For ‘amd64’ the builds have failed to complete in several different
ways, currently with relocation errors for libRblas.so.
So you want us to deliver an entire backtesting architecture for you?
That's a pretty hefty request...
But being R, it's already been done. This has all basically been made
available in the quantstrat project (google it) -- the documentation
is online and you can see lots of worked examples in
Just use a loop and collect the results in a list. BTW, R-Sig-Geo is a
special interest group mailing list specifically for these sorts of
questions regarding sp and rgeos packages (which you did not mention).
library(sp)
library(rgeos)
point - data.frame(x=c(10,20), y=c(10, 10), z = c(2,7))
On Apr 6, 2012, at 9:41 PM, wcheckle wrote:
i am trying to replicate the following graph using xyplot :
attach(x)
plot ( jitter(type), mortality, pch=16, xlim = c(0.25, 3.75))
lines ( c(1-0.375,1.375) , c ( median(mortality[type==1]),
median(mortality[type==1])), lwd=5,col=2)
lines (
readline(press return to continue)
worked as a pause until I used it in 4 places in a different program where
none of them wait.
What problem do I look for?
--
View this message in context:
http://r.789695.n4.nabble.com/pause-code-in-R-code-tp4523209p4539399.html
Sent from the R help mailing
Are you sure you are executing that line of code?
Michael
On Sat, Apr 7, 2012 at 9:18 AM, Hurr hill0...@umn.edu wrote:
readline(press return to continue)
worked as a pause until I used it in 4 places in a different program where
none of them wait.
What problem do I look for?
--
View
One last try -- though if you cannot figure out how to take the name
of your object and put it in dput( NameOfObjectGoesHere ) after I gave
you a worked example of how to do so, it seems unlikely that this will
help:
# Attempt to recreate your data
library(zoo)
x - read.zoo(textConnection(monoMn
On Sat, Apr 7, 2012 at 6:18 AM, Hurr hill0...@umn.edu wrote:
readline(press return to continue)
worked as a pause until I used it in 4 places in a different program where
none of them wait.
What problem do I look for?
Your bad coding.
cf the posting guide -- especially reproducible examples
Dear Waseem
On 7 April 2012 09:20, waseem ahmad waseem...@yahoo.com wrote:
I want to estimate simultaneous equation model with panel data.
The model looks as follows
Y1=a0+a1*X1+a2*X2
Y2=b0+b1*X2+b2*X1
X1=Z1-(Y1/a1)
X2=Z2-(Y2/b1)
I
In this model Y1, Y2, X1 and X2 are endogenous
Yes, that is the question.
So is it the readline() ?
--
View this message in context:
http://r.789695.n4.nabble.com/pause-code-in-R-code-tp4523209p4539615.html
Sent from the R help mailing list archive at Nabble.com.
__
R-help@r-project.org mailing
here is the data (fyi this is made-up data)
x = as.data.frame (
cbind(c(5,8,7,5,8,10,11,6,4,5,20,25,27,30,35,32,28,21,20,34,11,15,18,12,15,12,10,15,19,20),
c(1,1,1,1,1,1,1,1,1,1,2,2,2,2,2,2,2,2,2,2,3,3,3,3,3,3,3,3,3,3),c(1,0,1,1,1,1,0,0,0,1,0,1,0,0,0,0,1,0,0,0,1,1,1,0,1,0,1,0,0,0)))
I was trying to put some inequality constraints in the solnp package. But
whenever I put them into the package. It just showed that there is unused
arguments in it. So can anyone give me a simple example with inequality
constraints ??? i want to know the right form of writing the
I have a 11*1562 set of data. I want to regress using ols all combination of
the ten variables on the eleventh. Store the results of the t test of each
regression and the r2 to compare which combination is better. the
combination can use any amount o variables. How can I do that in R?
--
View
Thanks for help... do anybody also have an idea how to export object l (which
is list class) to .*shp file. I use SpatialPolygons, but output in console is
:
Error in SpatialPolygons(list(l)) :
cannot get a slot (area) from an object of type list
Cheers, OV
On 07.04.2012 17:27, kebrab67 wrote:
I have a 11*1562 set of data. I want to regress using ols all combination of
the ten variables on the eleventh. Store the results of the t test of each
regression and the r2 to compare which combination is better. the
combination can use any amount o
This is, more often that not, statistically speaking a bad idea:
you're likely to get a false positive (https://xkcd.com/882/)
But it sounds like something along these lines should work for you to
get the models :
x - data.frame(rnorm(50), rnorm(50), rnorm(50), rnorm(50))
m - vector(list,
On Apr 7, 2012, at 11:35 AM, wcheckle wrote:
here is the data (fyi this is made-up data)
x = as.data.frame (
cbind
(c
(5,8,7,5,8,10,11,6,4,5,20,25,27,30,35,32,28,21,20,34,11,15,18,12,15,12,10,15,19,20
),
c
(1,1,1,1,1,1,1,1,1,1,2,2,2,2,2,2,2,2,2,2,3,3,3,3,3,3,3,3,3,3
Hello,
I have the following code, but it seems that something must be doing
wrong, because it is giving the results I want.
The idea is to create segments while the value of Commutation is less than
1000.
for example, from the small set of data below
I believe that code was mine,
Dear All
I am trying to find a uniroot of a function within another function (see
example) but I am getting an error message (f()values at end points not of
opposite sign). I was wondering if you would be able to advise how redefine
my function so that I can find the solution. In short my first
On 07-04-2012, at 19:48, jeka12386 wrote:
Dear All
I am trying to find a uniroot of a function within another function (see
example) but I am getting an error message (f()values at end points not of
opposite sign). I was wondering if you would be able to advise how redefine
my function so
On Sat, Apr 7, 2012 at 11:16 AM, David Winsemius dwinsem...@comcast.net wrote:
xyplot(mortality ~ type, data=xdat,
panel=function(x,y){
panel.xyplot(x,y, jitter.x=TRUE)
panel.segments(x0=c(.9, 1.9, 2.9),
On 2012-04-07 08:27, kebrab67 wrote:
I have a 11*1562 set of data. I want to regress using ols all combination of
the ten variables on the eleventh. Store the results of the t test of each
regression and the r2 to compare which combination is better. the
combination can use any amount o
On Fri, Apr 6, 2012 at 9:42 PM, help ly help.ly2...@gmail.com wrote:
Dear All,
I would like to make a quadratic with a plateau model in R. Is there a
package in R doing this? The bentcableAR package seems won't work.
The link below describes what I am looking for in R exactly:
Hello,
using the - assignment operator I do not understand why the following does
not work.
l - list()
l
list()
l$arg1 - test
error in l$arg1 - test : Objekt 'l' not found
?- says: The operators - and - cause a search to made through the
environment for an existing definition of the variable
On Sat, Apr 7, 2012 at 1:30 PM, Mark Heckmann mark.heckm...@gmx.de wrote:
Hello,
using the - assignment operator I do not understand why the following does
not work.
l - list()
l
list()
l$arg1 - test
error in l$arg1 - test : Objekt 'l' not found
?- says: The operators - and - cause a
Hello All,
I'm getting my workflow switched over to Sweave, which is very cool.
However, I collaborate with folks (as many of you must as well) who use
Word to Track Changes amongst a group while crafting a paper. In the
simplest case, there will just be two people (one Sweave user and one
Word
Thanks! I'll try to stick to that advice!
Maybe there is a better way... Here is what I want:
I want to save some default settings for a package.
The user can change these using a function similar to par().
I do not want to use options() here as it will be quite a lot of parameters.
I was
I cannot reply in detail at the moment, but we have a proposal to the
Google Summer of Code this year which will address the collaboration
issue. The basic idea is to write everything in plain text with
markdown (including R code), and compile the file with the knitr
package, then convert it to
would it be possible to separate the parts generated by R from the
actually written paths? As creating small doc's by sweave which contain
the generated reports and embedding them via OLE into the static written
text. However, this could require sending around multiple documents
around each
On Sat, Apr 7, 2012 at 2:01 PM, Mark Heckmann mark.heckm...@gmx.de wrote:
Thanks! I'll try to stick to that advice!
Maybe there is a better way... Here is what I want:
I want to save some default settings for a package.
The user can change these using a function similar to par().
I do not
On Sat, Apr 7, 2012 at 3:58 PM, Gabor Grothendieck
ggrothendi...@gmail.com wrote:
On Fri, Apr 6, 2012 at 9:42 PM, help ly help.ly2...@gmail.com wrote:
Dear All,
I would like to make a quadratic with a plateau model in R. Is there a
package in R doing this? The bentcableAR package seems won't
On Sat, Apr 7, 2012 at 5:01 PM, Mark Heckmann mark.heckm...@gmx.de wrote:
Thanks! I'll try to stick to that advice!
Maybe there is a better way... Here is what I want:
I want to save some default settings for a package.
The user can change these using a function similar to par().
I do not
The eleventh variable is log(wage) and the 10 other variables are actually
all at least theorically significant to the regression. like wether the
person is black, immmigrant, it´s experience, age, if he works in a city...
But i need only take a few variables out of this ten. Actually taking only
Peter showed how to get the minimums from a list or data frame using
sapply, here is a way to copy your 1440 vectors into a single list
(doing this and keeping your data in a list instead of separate
vectors will make your life easier in general):
my.list - lapply( 1:1440, function(x) get(
On Apr 7, 2012, at 3:58 PM, kebrab67 wrote:
The eleventh variable is log(wage) and the 10 other variables are
actually
all at least theorically significant to the regression. like wether
the
person is black, immmigrant, it´s experience, age, if he works in a
city...
But i need only take a
The emp.hpd function in the TeachingDemos package will do this
(assumes a single interval result, either unimodal or multimodes but
the valleys between don't drop far enough to split the interval). I
am sure there are similar functions in other packages as well.
On Fri, Apr 6, 2012 at 12:39 PM,
One other minor glitch I just noticed:
x - seq_along(x)
should read
x - seq_along(y)
cheers,
Rolf
On 08/04/12 10:08, Gabor Grothendieck wrote:
On Sat, Apr 7, 2012 at 3:58 PM, Gabor Grothendieck
ggrothendi...@gmail.com wrote:
On Fri, Apr 6, 2012 at 9:42 PM, help
Yes, I had fixed that on my end but didn't notice I had copied the old
version. I think its obvious enough that anyone will fix it
themselves. Regards.
On Sat, Apr 7, 2012 at 7:23 PM, Rolf Turner rolf.tur...@xtra.co.nz wrote:
One other minor glitch I just noticed:
x - seq_along(x)
Dear Jenn,
Could you please provide a reproducible example. In your case, you could for
example provide the data (using dput(dat)).
Regards
Søren
-Original Message-
From: r-help-boun...@r-project.org [mailto:r-help-boun...@r-project.org] On
Behalf Of Jenn Barrett
Sent: 6. april 2012
On 12-04-07 4:51 PM, Henrik Bengtsson wrote:
On Sat, Apr 7, 2012 at 1:30 PM, Mark Heckmannmark.heckm...@gmx.de wrote:
Hello,
using the- assignment operator I do not understand why the following does not
work.
l- list()
l
list()
l$arg1- test
error in l$arg1- test : Objekt 'l' not found
?-
I see what you mean as I go trgouh my courses (chapters related to type i and
II errors) Is there any way in the therory for me to choose which variables
in a set of variables explain more my Y. Or should I base my reflexion
solely on economic theory and discriminate between variables?
--
View
I am not asking anyone to deliver anything, to the best of my understanding
all that I requested in terms of functionality is more than likely done in
quantstrat as you point out. Issues is: quantstrat has extremely poor
documentation. Yes, there are samples, but they are certainly not simple
and
On Apr 7, 2012 7:55 PM, sysot1t syso...@gmail.com wrote:
I am not asking anyone to deliver anything, to the best of my
understanding
You asked for documented code examples of the items you listed in your
original email.
all that I requested in terms of functionality is more than likely done
Thank you David, the bwplot option does what I need:
x11(height=8,width=11)
bwplot ( mortality~ attend|type,
pch=95,cex=5,col=2,
par.settings=list(
box.rectangle = list(col = transparent),
box.umbrella = list(col = transparent),
plot.symbol = list(col = transparent)
),
On Apr 7, 2012, at 10:29 PM, wcheckle wrote:
Thank you David, the bwplot option does what I need:
x11(height=8,width=11)
bwplot ( mortality~ attend|type,
pch=95,cex=5,col=2,
par.settings=list(
box.rectangle = list(col = transparent),
box.umbrella = list(col = transparent),
plot.symbol =
Hi Maria, you need to name your history file somefilename.Rapp.history,
example
savehistory(file = history.Rapp.history)
/nkschmidt
Maria Gouskova wrote
Dear R users,
I am having a minor but annoying issue with R.app. It doesn't retain
the history information from the previous
yes, I would have expected documented samples for something simple.. as I
said, I am not asking for anything complex but something rather simple that
I can use to learn from and build upon.. not a highly complex example that
provides me with little explanation of what is going on... quanstrat is
On 08-04-2012, at 03:31, nkschmidt wrote:
Hi Maria, you need to name your history file somefilename.Rapp.history,
example
savehistory(file = history.Rapp.history)
/nkschmidt
You replied to a message almost two years old, which to my mind is pretty
useless.
Original message is here:
66 matches
Mail list logo
