On Mon, Nov 4, 2013 at 7:27 AM, Enzo Cocca 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 = 10)
>
>
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 br
On Tue, 5 Nov 2013 16:22:26 -0700
"Angela Dwyer" 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
__
R-h
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 : http://www.londonm
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
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 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
> would be more than thankf
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 lo
[[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
http://www.branah.com/a
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
2013
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 functions
First of all, use "readBin" to verify you get the desired data back.
Second, that '00' is, I believe the 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
> migh
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 r
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), dat
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 strin
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 converte
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
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, E
(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 wrote:
> Hi all,
> the following returns
On Nov 6, 2013, at 10: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 stam
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 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 wondering
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 SN
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 column
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 a
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 the
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 prob
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 wrote:
> Is this what you are after:
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 t
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 th
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) :
t
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 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
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
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 bo
On Wed, Nov 6, 2013 at 5:44 PM, Ye Lin 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 companies d
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 unab
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 fin
On Nov 6, 2013, at 11:09 AM, Silvia Espinoza 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 there a paid v
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:
splittedStrings<-str_split(datafr
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 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
>
> > -Or
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 th
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 and
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 offer
The changepoint package might give you a way to do this.
-
David L Carlson
Department of Anthropology
Texas A&M University
College Station, TX 77840-4352
-Original Message-
From: r-help-boun...@r-project.org
[mailto:r-help-boun...@r-project.org] On Beh
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 regression
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 indivi
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.
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 runnin
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 cl
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
On 11/07/13 10:57, David Winsemius wrote:
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.
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 wrote:
> On 11/07/13 10:57, David Winsemius wrote:
>
>
>>
>> I think you need to add a statistician to your [PhD] committee. The
>> diffic
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 data(
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)
> data('studentnets.m
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 h
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 che
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
> 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 cod
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"
wrote:
>I am not trying to fit a horizontal line at every unique value of y. I am
>trying f
> 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.
Bil
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 general
>>> 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
>> wr
79 matches
Mail list logo
