On Mon, Nov 4, 2013 at 7:27 AM, Enzo Cocca enzo@gmail.com wrote:
yes barry I really need this.
I tried to use raster or rgdal but with poor results.
I have this function:
VGM_PARAM_A3 - gstat(id=bos_bison,
formula=combusto~1,locations=~coord_x+coord_y, data=archezoology_table,
nmax =
Hi,
After helpstart() and the first jump into a package, the nest jump will
keep trying endlessly Then I go back to the R window, where ' Making
'packages.html' ... done' is being displayed but no prompt. So I bring
back the prompt by ^C^C. Then, o wonder, the prompt comes back and in
the
On Tue, 5 Nov 2013 16:22:26 -0700
Angela Dwyer angela.dw...@rmbo.org wrote:
You didn't forget to load the library did you? The bit of output you
provide doesn't show a library(adehabitat) line. That needs to be run
before the function can be found.
JWDougherty
Is there a way which I can use a hidden function say f() of a package say
foo, (that is, f()is non exported in the NAMESPACE of foo) within the
package foo without using foo:::f()
Mikis
Prof Mikis Stasinopoulos
d.stasinopou...@londonmet.ac.uk
Companies Act 2006 :
Dear R project officials,
I have found that in R 3.0.1 version writeBin function of base package
might not work correctly. For command writeBin(100,raw()) it answers 31 30
30 00 the last double 0 is differs from http://www.branah.com/ascii-converter
there ascii codes are 31 30 30. So is it
Hello,
hope all are well. I just wanted some help with finding some of the
packages that could be used for Fraud and Anomaly detection purposes. I
would be more than thankful if some packages might be referenced which
could be in the telecommunications industry.
Regards,
Hossam Hassanien
take a look at zoo and outliers package. may be help you
On Wed, Nov 6, 2013 at 3:00 AM, Hossam Hassanien hossam.hassan...@gmail.com
wrote:
Hello,
hope all are well. I just wanted some help with finding some of the
packages that could be used for Fraud and Anomaly detection purposes. I
Dear Experts:
I am very appreciate your comments and help!
Actually I am a new comer from MATLAB. If the function
can see global variables, then it may output wrong results without
any error messages. For example, there is a gloabl variable named
v, and I write one funciton with one
[[alternative HTML version deleted]]
__
R-help@r-project.org mailing list
https://stat.ethz.ch/mailman/listinfo/r-help
PLEASE do read the posting guide http://www.R-project.org/posting-guide.html
and provide commented, minimal, self-contained,
On 06/11/2013 3:35 AM, Harutyun Khachatryan wrote:
Dear R project officials,
I have found that in R 3.0.1 version writeBin function of base package might not work correctly. For
command writeBin(100,raw()) it answers 31 30 30 00 the last double 0 is differs from
Could you please approve this question and let it be sent in R_help?
addendum on the question: it also fails when using a data.frame in guess,
so I have no clue
guess(grnn, data.frame(x1=c(2), x2=c(4)))
Error in (X - Xa) %*% t(X - Xa) :
requires numeric/complex matrix/vector arguments
On 11/4/2013 10:45 AM, Sergio Fonda wrote:
Hi,
I'm not able to get information about the following question:
is the variables standardization a default option in manova() (stats package)?
Or if you want to compare variables with different units or scales and
rather different variances, you have
Dear Zhong-Yuan Zhang,
R is lexically scoped. Pretending that you're using a different programming
language is probably a bad idea.
The findGlobals() function in the codetools package, which is part of the
standard R distribution, can help you locate references to global variables
(and
First of all, use readBin to verify you get the desired data back.
Second, that '00' is, I believe the EOF character you'll find at the end
of any file.
Harutyun Khachatryan wrote
Dear R project officials,
I have found that in R 3.0.1 version writeBin function of base package
might not
See
http://cran.r-project.org/mirror-howto.html
Best,
Uwe Ligges
On 06.11.2013 13:34, Abhinav Kashyap wrote:
[[alternative HTML version deleted]]
__
R-help@r-project.org mailing list
https://stat.ethz.ch/mailman/listinfo/r-help
PLEASE do
You can reproduce the problem by having a data.frame (or anything else) in your
environment:
new - data.frame(ï..VAR1 = rep(c(TRUE,NA,FALSE), c(10,2,8)),
random=rep(1:3,len=20), clustno=rep(c(1:5),len=20),
validatedRS6=rep(0:1,len=20))
model1- glmer(validatedRS6 ~ random + (1|clustno),
Sorry, Carl, but you missed the boat on both responses.
Using readBin to read what they wrote won't help the OP if they don't
understand what they are writing. Nor is byte 0 the EOF marker on any
operating system I have ever used. (It does happen to be the string terminator
for in-memory
Hi all,
the following returns the hour and the minutes
paste(DataSet$TimeStamps[selectedInterval$start,4],
DataSet$TimeStamps[selectedInterval$start,5],sep=:)
[1] 12:3
the problem is that from these two I want to create a time stamp so 12:03. The
problem is that the number 3 is not converted
On Nov 6, 2013, at 8:25 AM, Alaios wrote:
Hi all,
the following returns the hour and the minutes
paste(DataSet$TimeStamps[selectedInterval$start,4],
DataSet$TimeStamps[selectedInterval$start,5],sep=:)
[1] 12:3
the problem is that from these two I want to create a time stamp so 12:03.
Dear administrators,
I tried to stop receving mails from your webpage to my current mail (
gaiarr...@usal.es). I try to do it through your webpage but i was not
able. I would appreciate to be removed from the mailing list.
Thanks in advance.
--
Mario Garrido Escudero
Dpto. de Biología Animal,
(Assuming I understand) tons of ways of doing this.
So I'll just point out the
?nchar
function, which you can use to count characters in your tail end and
paste a 0 if there's only one, e.g. via ifelse() .
-- Bert
On Wed, Nov 6, 2013 at 8:25 AM, Alaios ala...@yahoo.com wrote:
Hi all,
the
On Nov 6, 2013, at 10:25 AM, Alaios ala...@yahoo.com wrote:
Hi all,
the following returns the hour and the minutes
paste(DataSet$TimeStamps[selectedInterval$start,4],
DataSet$TimeStamps[selectedInterval$start,5],sep=:)
[1] 12:3
the problem is that from these two I want to create a time
Hi all,
I have estimated a random panel model using plm function.
I have a question about the vector of resduals obtained with the
object $residuals.
example:
data(Produc, package = plm)
zz - plm(log(gsp) ~ log(pcap) + log(pc) + log(emp) + unemp,
model=random, data = Produc, index =
Hello all,
I just noticed the following behavior of plot:
x - c(1,2,9)
plot(x ~ x) # this is just like doing:
plot(x)
# when maybe we would like it to give this:
plot(x ~ c(x))
# the same as:
plot(x ~ I(x))
I was wondering if there is some reason for this behavior.
Thanks,
Tal
You can reproduce the problem by having a data.frame (or anything else) in
your
environment:
I left out called 'new' in the above statement. The example is correct.
Bill Dunlap
Spotfire, TIBCO Software
wdunlap tibco.com
-Original Message-
From: r-help-boun...@r-project.org
On Nov 6, 2013, at 10:40 AM, Tal Galili tal.gal...@gmail.com wrote:
Hello all,
I just noticed the following behavior of plot:
x - c(1,2,9)
plot(x ~ x) # this is just like doing:
plot(x)
# when maybe we would like it to give this:
plot(x ~ c(x))
# the same as:
plot(x ~ I(x))
I was
It probably happens because plot(formula) makes one call to terms(formula) to
analyze the formula. terms() says there is one variable in the formula,
the response, so plot(x~x) is the same a plot(seq_along(x), x).
If you give it plot(~x) , terms() also says there is one variable, but
no response,
Hi,
I have a dataset containing approximately 2 million SNPs. The data is in PLINK
format. The data conversion was successful except for a warning message Nan
introduced due to coercion. My problem is that I get an error message when I
perform a PCA:
Principal Component Analysis (PCA) on SNP
Dear experts,
I’ve been on this for weeks now, and couldn’t find a solution..Sorry for the
long description. I figured I post many details, so you get the problem
entirely, although it’s not hard to grasp.
**Situation:**
Data frame consisting of 4 million entries (total size: 250 MB). Two
I have two columns for speed ('Smoothed velocity') and Spacing. What I want to
do is to first create the intervals of speed (minimum value=0, max value=
85.53), group the Spacing values falling in a particular Speed interval, find
the average of the Spacing for an interval and finally plot the
All these suggestions of using 'sprintf' might be right but you might
be doing it wrong...
If you are working with times, then use the date/time classes and the
handy functions for working on them. Which means the lubridate
package, most likely.
Are these times part of a calendar time, or are
I'm having trouble running the ode function from the deSolve package.
I am running RStudio under Ubuntu 13.1
I am running ode() on compiled code that returns delta values using the .C
convention. While I can include an example of the code, I suspect that it
will not be helpful, since the problem
Is this what you are after:
n - 1000
x - data.frame(speed = runif(n, 0, 85.53)
+ , spacing = rnorm(n)
+ )
# create interval for speed
int - seq(0, max(x) + 4.5, by = 4.5)
# split the data and find average of spacing
tapply(x$spacing, cut(x$speed, int), mean)
should of had in the script:
int - seq(0, max(x$speed) + 4.5, by = 4.5)
Jim Holtman
Data Munger Guru
What is the problem that you are trying to solve?
Tell me what you want to do, not how you want to do it.
On Wed, Nov 6, 2013 at 12:19 PM, jim holtman jholt...@gmail.com wrote:
Is this
I am not trying to fit a horizontal line at every unique value of y. I am
trying fit the y values with as few horizontal lines by trading off the
number of horizontal lines with the error. The actual problem I am trying
to solve is to smooth data in a time series. Here is a realistic example of
y
Good morning. I am interested in downloading R. I would appreciate if you
can help me with the following questions, please.
1. Is R free, or I have to pay for support/maintenance, or it depends
on the version? Is there a paid version?
2. How safe is it to work with data using R? Is
On 06/11/2013 17:20, Adam Clark wrote:
I'm having trouble running the ode function from the deSolve package.
I am running RStudio under Ubuntu 13.1
I am running ode() on compiled code that returns delta values using the .C
convention. While I can include an example of the code, I suspect that
Hi
I have managed to install MPICH2 and Rmpi on my Windows 7 machine. I can
also run the following code
library(Rmpi)
mpi.spawn.Rslaves()
4 slaves are spawned successfully. 0 failed.
master (rank 0, comm 1) of size 5 is running on: MyMaster
slave1 (rank 1, comm 1) of size 5
Interestingly, fitting an LM with x on both sides gives a warning, and
then drops it from the RHS, leaving you with just an intercept:
lm(x~x,data=d)
Call:
lm(formula = x ~ x, data = d)
Coefficients:
(Intercept)
4
Warning messages:
1: In model.matrix.default(mt, mf, contrasts) :
You can get details at http://www.r-project.org/
But to answer your question: Yes it is free
On Wed, Nov 6, 2013 at 9:09 AM, Silvia Espinoza siles...@gmail.com wrote:
Good morning. I am interested in downloading R. I would appreciate if you
can help me with the following questions, please.
Addendum:
unloading and reloading deSolve.so does indeed fix the problem:
library.dynam.unload(deSolve, libpath=paste(.libPaths()[1], //deSolve,
sep=))
library.dynam(deSolve, package=deSolve, lib.loc=.libPaths()[1])
However, this is a little clunky, and seems like overkill. Does anybody
have an
Dear list users,
I transformed two vectors of seasonal data in ts objects of frequency 4:
y1 - ts(x1, frequency=4, start=c(1952,1))
y2 - ts(x2, frequency=4, start=c(1952,1))
In this way Qtr1 corresponds to Winters, Qtr2 corresponds to Springs and so on.
I would like to plot on the same graph
On Wed, Nov 6, 2013 at 5:44 PM, Ye Lin ye...@lbl.gov wrote:
You can get details at http://www.r-project.org/
But to answer your question: Yes it is free
But there is also a paid version. Send me $1000 and I will send you R
on a USB stick, complete with all the source code.
Seriously, other
On Nov 6, 2013, at 8:33 AM, Mario Garrido wrote:
Dear administrators,
I tried to stop receving mails from your webpage to my current mail (
gaiarr...@usal.es). I try to do it through your webpage but i was not
able. I would appreciate to be removed from the mailing list.
If you were unable
You need to read up on the use of the apply functions: this is just
another 'tapply' call with the 'speed' in it.
n - 1000
x - data.frame(speed = runif(n, 0, 85.53)
+ , spacing = rnorm(n)
+ )
# create interval for speed
int - seq(0, max(x) + 4.5, by = 4.5)
# split the data and find average
On Nov 6, 2013, at 11:09 AM, Silvia Espinoza siles...@gmail.com wrote:
Good morning. I am interested in downloading R. I would appreciate if you
can help me with the following questions, please.
1. Is R free, or I have to pay for support/maintenance, or it depends
on the version? Is
Good idea!
I'm trying your approach right now, but I am wondering if using str_split
(package: 'stringr') or strsplit is the right way to go in terms of speed? I
ran str_split over the text column of the data frame and it's processing for
2 hours now..?
I did:
I'll answer myself:
using strsplit with fixed=true took like 2minutes!
--
View this message in context:
http://r.789695.n4.nabble.com/speed-issue-gsub-on-large-data-frame-tp4679747p4679905.html
Sent from the R help mailing list archive at Nabble.com.
William Dunlap wdunlap at tibco.com writes:
You can reproduce the problem by having a data.frame (or anything
else) in your environment: I left out called 'new' in the
above statement. The example is correct.
Bill Dunlap
Spotfire, TIBCO Software
wdunlap tibco.com
-Original
If you could, please identify which responder's idea you used, as well as the
strsplit -- related code you ended up with.
That may help someone who browses the mail archives in the future.
Carl
SPi wrote
I'll answer myself:
using strsplit with fixed=true took like 2minutes!
--
View this
Hi,
I am a graduate student applying published R scripts to compare the
classification accuracy of 2 predictive models, one built using discriminant
function analysis and one using random forests (webpage link for these scripts
is provided below). The purpose of these models is to predict the
Hi there
I've got a (I think) simple problem, but I just can't figure it out.
I've got 3 dataset (datasets attached below).
In 'Value' there is the measured mean value for every unit (SFZN, SGKN,
SGSN, SHLTN, SIK) and for every treatment (1 to 10).
In 'SE' is the Standard Error for every unit
Hello,
I'm having fun exploring the pretty graphing options in R, although
I'm struggling to figure out how to do some simple things; would be
thankful if someone could point me toward relevant sections of the manual
or provide some starter code to get me going.
I'd like to extend what is
The changepoint package might give you a way to do this.
-
David L Carlson
Department of Anthropology
Texas AM University
College Station, TX 77840-4352
-Original Message-
From: r-help-boun...@r-project.org
[mailto:r-help-boun...@r-project.org] On
Greetings, My question is more algorithmic than prectical. What I am
trying to determine is, are the GAM algorithms used in the mgcv package
affected by nonnormally-distributed residuals?
As I understand the theory of linear models the Gauss-Markov theorem
guarantees that least-squares
Hi experts
I would like to run Hamming window on a time series in r. I am using e1071
package and STFT function in it like this:
s-stft(datalist, win=min(80,floor(length(datalist)/10)),
inc=min(24,floor(length(datalist)/30)), coef=10, wtype=hanning.window)
My expectation is so, that I get a
Hi,
The data (attached) I am looking at consists of measurements of growth rate
at different ages, for individuals in two treatments (control and infected).
What I want to know is whether and when (what age) the growth rate of
infected individuals is higher than the growth rate for control
Dear R help community,
I have a version of the data below generated by counting the number of
individuals in each state every 2-3 time units for a series of cohorts (sets).
States are A, B, C, E and T.
It is known that the following transitions are possible A-B, B-C, C-E, B-T,
C-T and E-T.
The
Hi
I have managed to install MPICH2 and Rmpi on my Windows 7 machine. I can
also run the following code
library(Rmpi)
mpi.spawn.Rslaves()
4 slaves are spawned successfully. 0 failed.
master (rank 0, comm 1) of size 5 is running on: MyMaster
slave1 (rank 1, comm 1) of size 5 is running
On Nov 6, 2013, at 10:07 AM, Henderson, Robin Michelle wrote:
Hi,
I am a graduate student applying published R scripts to compare the
classification accuracy of 2 predictive models, one built using discriminant
function analysis and one using random forests (webpage link for these
On 11/6/2013 6:50 PM, Adam Clark wrote: Addendum:
unloading and reloading deSolve.so does indeed fix the problem:
library.dynam.unload(deSolve, libpath=paste(.libPaths()[1],
//deSolve, sep=))
library.dynam(deSolve, package=deSolve, lib.loc=.libPaths()[1])
However, this is a little clunky, and
On Nov 6, 2013, at 9:19 AM, Sashikanth Chandrasekaran wrote:
I am not trying to fit a horizontal line at every unique value of y. I am
trying fit the y values with as few horizontal lines by trading off the
number of horizontal lines with the error. The actual problem I am trying
to solve is
On 11/07/13 10:57, David Winsemius wrote:
SNIP
I think you need to add a statistician to your committee. The
difficulties you are facing (of which you appear to be unaware) are
not just related to being new to R.
SNIP
Fortune?
cheers,
Rolf
On Nov 6, 2013, at 12:46 PM, Collin Lynch wrote:
Greetings, My question is more algorithmic than prectical. What I am
trying to determine is, are the GAM algorithms used in the mgcv package
affected by nonnormally-distributed residuals?
As I understand the theory of linear models the
Why would you need to? The whole point of :: and ::: is to specify the
origin of a function.
--
View this message in context:
http://r.789695.n4.nabble.com/hidden-functions-tp4679849p4679856.html
Sent from the R help mailing list archive at Nabble.com.
Second (perhaps with the slight addition indicated)
-- Bert
... And **amen!** to the sentiment expressed.
On Wed, Nov 6, 2013 at 2:08 PM, Rolf Turner r.tur...@auckland.ac.nz wrote:
On 11/07/13 10:57, David Winsemius wrote:
SNIP
I think you need to add a statistician to your [PhD]
Hi:
I have installed the NetData package, and want to use the MagAct96-98
- Extracurricular affiliation data. But while loading the data, its
giving an error. Any help?
install.packages(NetData)
library(NetData)
data('studentnets.magact96.97.98', package = NetData)
Warning message:
In
On Nov 6, 2013, at 1:24 PM, email wrote:
Hi:
I have installed the NetData package, and want to use the MagAct96-98
- Extracurricular affiliation data. But while loading the data, its
giving an error. Any help?
install.packages(NetData)
library(NetData)
Seen on StackOverflow. Any seconds?
I would heed the warnings and diagnostics. They are there for a reason. The
Ostrich algorithm does not help you.
-- Dirk Eddelbuettel commenting on StackOverflow when a questioner said he
had not run R CMD check because he suspected other problems would
On 11/07/2013 07:46 AM, c_e_cressler wrote:
Hi,
The data (attached) I am looking at consists of measurements of growth rate
at different ages, for individuals in two treatments (control and infected).
What I want to know is whether and when (what age) the growth rate of
infected individuals is
On 11/07/13 13:19, David Winsemius wrote:
Seen on StackOverflow. Any seconds?
I would heed the warnings and diagnostics. They are there for a reason. The Ostrich
algorithm does not help you.
-- Dirk Eddelbuettel commenting on StackOverflow when a questioner said he
had not run R CMD
On Thu, 7 Nov 2013, Rolf Turner wrote:
On 11/07/13 13:19, David Winsemius wrote:
Seen on StackOverflow. Any seconds?
I would heed the warnings and diagnostics. They are there for a reason.
The Ostrich algorithm does not help you.
-- Dirk Eddelbuettel commenting on StackOverflow when a
Dear John Fox:
I highly appreciate
your help!!! Problems solved.
Best Wishes Always.
2013/11/6 John Fox j...@mcmaster.ca
Dear Zhong-Yuan Zhang,
R is lexically scoped. Pretending that you're using a different programming
language is probably a bad idea.
The findGlobals() function
Possibly see the
strucchange
package
--
Don MacQueen
Lawrence Livermore National Laboratory
7000 East Ave., L-627
Livermore, CA 94550
925-423-1062
On 11/6/13 9:19 AM, Sashikanth Chandrasekaran
sashikanth.chandraseka...@gmail.com wrote:
I am not trying to fit a horizontal line at every
The default functional link for mgcv::gam is log, so I doubt that
your theoretical understanding applies to GAM's in general. When Simon
Wood wrote his book on GAMs his first chapter was on linear models, his
second chapter was on generalized lienar models at which point he had
written over
On 11/07/2013 04:56 AM, Stefano Sofia wrote:
Dear list users,
I transformed two vectors of seasonal data in ts objects of frequency 4:
y1- ts(x1, frequency=4, start=c(1952,1))
y2- ts(x2, frequency=4, start=c(1952,1))
In this way Qtr1 corresponds to Winters, Qtr2 corresponds to Springs and so
y2 - ts(x2, frequency=4, start=c(1952,1))
y2w-ts(y2[seq(1,61,by=4)],frequency=1,start=1952)
I think it is simpler to compute y2w using the window() function, which
figures out the right call to seq() for you:
window(y2, frequency=1, start=c(1952,1)) # use c(1952,2) for springs, etc.
Bill
On Nov 6, 2013, at 5:44 PM, Collin Lynch wrote:
The default functional link for mgcv::gam is log, so I doubt that
your theoretical understanding applies to GAM's in general. When Simon
Wood wrote his book on GAMs his first chapter was on linear models, his
second chapter was on generalized
The default functional link for mgcv::gam is log, so I doubt that
your theoretical understanding applies to GAM's in general. When Simon
Wood wrote his book on GAMs his first chapter was on linear models, his
second chapter was on generalized lienar models at which point he had
written over
79 matches
Mail list logo
