I used R with WinEdt some years ago (I think it was on 98, may be even XP),
and then moved on to other editors --- (X)Emacs, Tinn-R and Rcommander seem
to do the job reasonably well with differing functionality. I had sent
several requests to the WinEdt development team that I may even upgrade to a
The simplest way is:
xxx <- with(clyde,tapply(CONTTIME,SCRNO,sum))
You could also do:
xxx <- by(clyde,clyde[["SCRNO"]],function(x){sum(x[["CONTTIME"]])})
but this gives somewhat messy output; the aforesaid output may be
convenient for some purposes, not for others.
cheers,
Rolf
This is Googlable. You need to add the CRAN repository to your APT sources file.
---
Jeff Newmiller The . . Go Live...
DCN: Basics: ##.#. ##.#. Live Go...
Live: OO#.. Dead: OO#.. Playing
Research Engineer (Solar/Batteri
On Wed, 29 Jun 2011, Sam Albers wrote:
I know that this has been asked before in other variations but I just can't
seem to figure out my particular application from previous posts. My
apologies if I have missed the answer to this question somewhere in the
archives. I have indeed looked.
I am ru
I tried this but did not work:
z0<- by(z, z[,"n"], function(x) subset(x, sum(n==0)>2))
Peter Maclean
Department of Economics
UDSM
- Original Message
From: Duncan Murdoch
To: Peter Maclean
Cc: r-help@r-project.org
Sent: Wed, June 29, 2011 3:33:25 PM
Subject: Re: [R] DROP OBSEVATION I
Dear All,
I wanted to install the reshape package which in turn requires the plyr.
When I tried to install plyr it says it needs
ERROR: this R is version 2.10.1, package 'plyr' requires R >= 2.11.0
My question is how do I upgrade my R ? I have Ubuntu 10.04.2 LTS.
Many thanks for your help.
I am estimating a large model by groups. How do you save the results
and returns
the associated quantiles?
For this example I need a data frame
n xi mubeta
1 0.1033614 2.5389580 0.9092611
2 0.3401922 0.5192882 1.5290615
3 0.5130798 0.5668308 1.2105666
I also want to
Hello!
I know that probably my question is rather simple but I' m a very beginner
R-user.
I have to numerically integrate the product of two function A(x) and B(x).
The integretion limits are [X*; +inf]
Function A(x) is a pdf function while B(x)=e*x is a linear function whose
value is equal to
Hi all,
My file that worked with a thin matrix (with few rows and colums) finally
made an error but I can´t comprehed why?
valor ausente donde TRUE/FALSE es necesario (absent value TRUE/FALSE is
necessary?)
Error en if (data[i, j] - last1[1, j] != 0) data2[i, j] = 0 else { :
You recommed to p
If you just want the value of the derivative at a particular point, would
numerical derivatives suffice? If so, try (for example) the numDeriv package.
S
From: r-help-boun...@r-project.org [r-help-boun...@r-project.org] On Behalf Of
Lisa [lisa...@gmail.co
Thankyou very much. I think "try" works for me.
I am learning it .
Sirius
--
View this message in context:
http://r.789695.n4.nabble.com/optimization-in-for-loop-tp3633638p3634100.html
Sent from the R help mailing list archive at Nabble.com.
__
R-hel
Hello,
I have the following dataset (this is a piece of a much larger set):
RID SCRNO VISCODE RECNO CONTTIME
2318 HBA0190012 bl 15
2418 HBA0190012 bl 35
2818 HBA0190012 bl 55
2918 HBA0190012 bl 25
321
- Forwarded Message -
From: Ungku Akashah
To: "r-help@r-project.org"
Sent: Thursday, June 30, 2011 9:14 AM
Subject: volcano plot.r
Hello.
My name is Akashah. i work at metabolic laboratory. From my study, i found that
volcano plot can help a lot in my section.
i already studied abo
On Wed, Jun 29, 2011 at 4:35 PM, Lisa wrote:
> Yes. I need to do implicit differentiation. After rearrangement, I got
>
> (x2 – x1) * b = log(1 / y - 1)
>
> Take derivative of both sides with respect to y, I have
>
> (x2 – x1) * b’[y] = - 1/y(1-y)
>
> Since both (x2 – x1) and b’[y] are vectors, I
Hello.
My name is Akashah. i work at metabolic laboratory. From my study, i found that
volcano plot can help a lot in my section.
i already studied about the volcano plot and get the coding to run in R
software, unfortunately, there is may be something wrong with the coding. This
is because no
Hi:
It's not at all clear to me what you want to do. Is each component of
your list a data frame with the same structure? Are you trying to
combine them into one data frame? If so, then try
mydata <- do.call(rbind, c)
where c is the name of the list. (Not a good choice of name, by the
way - c()
I know that this has been asked before in other variations but I just can't
seem to figure out my particular application from previous posts. My
apologies if I have missed the answer to this question somewhere in the
archives. I have indeed looked.
I am running Ubuntu 11.04, with R 2.12.1 and ESS+
Hi:
Your object b is a 5 x 12 matrix. The error message says that ets() is
expecting a univariate time series as its first argument.
Try something like
d <- ts(a[, 3], start = c(2005, 1), frequency = 12)
fit <- ets(d)
and see if that works. Untested since no reproducible example was provided.
H
What type of computer do you have now: operating system, memory.
Here is how long it took for me to read in a file with 4M lines and 13
columns of numerics on each line:
> system.time(x <- scan('/temp/large.txt', what = 0))
Read 52614432 items
user system elapsed
23.670.67 24.39
> str
People,
I have output from a Ruby script saved in JSON format and I can import
it into an R vector (list?) that looks like:
tst
[[1]]
[[1]][[1]]
[1] "01.01.01.00"
[[1]][[2]]
[1] -2.304248
[[2]]
[[2]][[1]]
[1] "01.01.01.01"
[[2]][[2]]
[1] -2.288097
[[3]]
[[3]][[1]]
[1] "01.01.01.02"
[[
I have the code below but the p-values are not what I am expecting. I want
to use the Cochran Armitage trend test to get the p-values.
What do you expect? Depending on the genetic model, you may not see HWE in the
"cases".
datamat[h,] <- t(rmultinom(1, size=c(10, 40, 50), prob=c(0.33, 0.33,
On Jun 29, 2011, at 4:56 PM, katarv wrote:
Hi Sarah,
I use readMat function, from R.matlab function.
str(X) tells that
X : num [1:64, 1:64, 1:21, 1:300] as I said, the first 3 columns
are x,y,z
coordinates. And I need all values in the last column correcponding
to a
given (x,y,z) coord
Hi
I have imported some files to a list, called "importa".
There are 43 files (importa[1], importa[2], ..., importa[43]). Now, I'm
trying to create a new table in "r". Each table will got the data from the
importa partition. Like file1<-importa[1].
I tried two ways:
c<-list()
for (i in 1:43) {
c[i]
Hi everyone,
I'm a newbie and this is my first post.
My boss wants me to make a series of scatter plots where 76 variables are
plotted against each other. I know how to do this using pair()...my problem
is that there are just too many plots to fit in the window.
Is there any way I can get all
Hi Sarah,
I use readMat function, from R.matlab function.
str(X) tells that
X : num [1:64, 1:64, 1:21, 1:300] as I said, the first 3 columns are x,y,z
coordinates. And I need all values in the last column correcponding to a
given (x,y,z) coordinate.
if you list the values of X, then they ar
Thank you for replying. when I've tried to run the R syntax in a 64 bit
computer,the problem is solved. Thank you for helping out. I totally agree
your advice.
I would like to answer all your questions in case other people meet the same
problem. The data contains one timestamp column with time zo
Hi all,
I have a function written by me that read a matrix (data frame) from a txt
with 4 million of rows and 13 columns.
The think is my function works with an input matrix of 100x13 and now I
tried to execute my function with the big "input file" and it is running
form the moment two hours...
Hi,
I am trying to put together a biplot using symbols and different colours
instead of text as points.
Someone has previously suggested using this code:
PC <- prcomp (iris[,1:4])
lambda <- PC$sdev * sqrt(nrow(PC$x))
plot (t(t(PC$x)/lambda),pch=16,col=as.numeric(iris[,5]))
par (new=T)
Rot <- t(t
--- Begin Message ---
This may be useful; convert them to NAs then omit them afterwards, you
end up keeping the positions where you had 0s.
> n[y==0] <- NA
> n
[1] 1 1 1 1 1 1 2 2 NA NA NA 2 3 3 NA NA 3 3
> na.omit(n)
[1] 1 1 1 1 1 1 2 2 2 3 3 3 3
attr(,"na.action")
[1] 9 10 11
Yes. I need to do implicit differentiation. After rearrangement, I got
(x2 – x1) * b = log(1 / y - 1)
Take derivative of both sides with respect to y, I have
(x2 – x1) * b’[y] = - 1/y(1-y)
Since both (x2 – x1) and b’[y] are vectors, I cannot move (x2 – x1) to
RHS. This is why I posted my quest
I forgot to include Vincent Goulet's web site in my previous email
http://vgoulet.act.ulaval.ca/en/emacs/
Best Regards
John
On 29 June 2011 21:12, Patrick Connolly wrote:
> On Wed, 29-Jun-2011 at 02:25PM +0200, Aditya Bhagwat wrote:
>
> |> Dear,
> |>
> |> How do I tell Emacs to update to the n
If you are using Windows Vincent Goulet has an excellent emacs install
for windows which contains ESS, Auctex and a lot of other extras. It
also contains well annotated configuration files which have specific
comments on how to change the version of R thqt emacs is calling. Not
being that expert
David Winsemius comcast.net> writes:
>
> The strategy of cbind vectors as an argument to data.frame and then
> naming them seems wasteful and error prone. Why not:
>
> z <- data.frame(n=factor(n),y=y)
> # all one step, no issues about every element needing to be the same
> mode
> # and no
zubin bellsouth.net> writes:
>
> Hello, running a mixed model in the package LME4, lmer()
>
> Panel data, have about 322 time periods and 50 states, total data set is
> approx 15K records and about 20 explanatory variables. Not a very
> large data set.
>
> We run random intercepts as well a
On Jun 29, 2011, at 4:29 PM, Peter Maclean wrote:
People with more experience in R I need help on this.
I would like to drop observation if they meet certain condition. In
this example
I would like to drop group 2 in "n" because the group in "Y" has
more than 2
zeroes.
#Example
n <- c(1, 1
On 29/06/2011 4:29 PM, Peter Maclean wrote:
People with more experience in R I need help on this.
I would like to drop observation if they meet certain condition. In this example
I would like to drop group 2 in "n" because the group in "Y" has more than 2
zeroes.
#Example
n<- c(1, 1, 1, 1, 1, 1,
People with more experience in R I need help on this.
I would like to drop observation if they meet certain condition. In this
example
I would like to drop group 2 in "n" because the group in "Y" has more than 2
zeroes.
#Example
n <- c(1, 1, 1, 1, 1, 1, 2, 2, 2, 2, 2, 2, 3, 3, 3, 3, 3,3)
y <-
Dennis:
It walks for small datset.
Peter Maclean
Department of Economics
UDSM
- Original Message
From: Dennis Murphy
To: Peter Maclean
Cc: r-help@r-project.org
Sent: Mon, June 27, 2011 8:37:00 PM
Subject: Re: [R] BY GROUP IN GEV
HI:
Since you didn't provide a reproducible example
On 30/06/11 06:16, Gabor Grothendieck wrote:
On Tue, Jun 28, 2011 at 10:03 PM, Lisa wrote:
Dear all,
I just want to get the derivative of a function that looks like:
y = exp(x1*b) / (exp(x1*b) + exp(x2*b))
where y is a scalar, x1, x2, and b are vectors. I am going to take the
derivative of b
On Wed, 29-Jun-2011 at 02:25PM +0200, Aditya Bhagwat wrote:
|> Dear,
|>
|> How do I tell Emacs to update to the new R version I installed? It still
|> loads the old R version. I already updated the the system path, but that
|> didn't seem to work.
Depends on information you didn't supply.
|>
|
I have a problem using RWinEdt 1.8.2 in Windows 7 Professional (64 bit).
System/software info:
R version 2.13.0 (2011-04-13)
Copyright (C) 2011 The R Foundation for Statistical Computing
ISBN 3-900051-07-0
Platform: x86_64-pc-mingw32/x64 (64-bit)
WinEdt Build: 20071003 (v. 5.5)
After installing
Hello, running a mixed model in the package LME4, lmer()
Panel data, have about 322 time periods and 50 states, total data set is
approx 15K records and about 20 explanatory variables. Not a very
large data set.
We run random intercepts as well as random coefficients for about 10 of
the variab
On 29/06/2011 3:15 PM, Janko Thyson wrote:
On 29.06.2011 20:58, Duncan Murdoch wrote:
On 29/06/2011 2:24 PM, Janko Thyson wrote:
Dear list,
this is not directly an R question, but it is somewhat related to R
aspects, so I hope it's okay to post it here:
I'd like to update my windows PATH base
On 29.06.2011 20:58, Duncan Murdoch wrote:
On 29/06/2011 2:24 PM, Janko Thyson wrote:
Dear list,
this is not directly an R question, but it is somewhat related to R
aspects, so I hope it's okay to post it here:
I'd like to update my windows PATH based on a script routine in order to
make sure
On Jun 29, 2011, at 2:31 PM, siriustar wrote:
Hi, dear R help
I am trying to use optim inside a for loop:
##For example. a: intial guess. b: result. f: function to be
minimized
for (i in 1:10) {
b[i] <- optim(a[i], f)}
However, some intial values cause error in optim function (e.g. "
On Jun 29, 2011, at 2:37 PM, Filipe Leme Botelho wrote:
Hi Tom,
At least to me it´s hard to picture what´s wrong without further
details regarding your data. I use spline/linear interpolation of
time series regularly, so maybe this example help you out.
ci_x
[1] 1 4 69 131 194 256
update.packages has suddenly stopped working for me (after working fine
yesterday). My default mirror is at Oregon State University, but I get the
following warning even after setting another mirror. Other Internet
connections are working on my computer. Suggestions, please.
> update.packages
Hi, dear R help
I am trying to use optim inside a for loop:
##For example. a: intial guess. b: result. f: function to be minimized
for (i in 1:10) {
b[i] <- optim(a[i], f)}
However, some intial values cause error in optim function (e.g. " system is
computationally singular..."). Then the for
Hi Tom,
At least to me it´s hard to picture what´s wrong without further details
regarding your data. I use spline/linear interpolation of time series
regularly, so maybe this example help you out.
> ci_x
[1] 1 4 69 131 194 256 320 382
> ci_y
[1] 0.1211 0.1213 0.1233 0.1241 0.1250 0.1254 0
On 29/06/2011 2:24 PM, Janko Thyson wrote:
Dear list,
this is not directly an R question, but it is somewhat related to R
aspects, so I hope it's okay to post it here:
I'd like to update my windows PATH based on a script routine in order to
make sure that crucial components are contained. Much
On Jun 29, 2011, at 2:29 PM, Christopher T. Moore wrote:
Hello,
I have encountered some unexpected behavior in R that seems to occur
as a result of having the current year embedded in a number:
No. that is not the explanation.
#Some large number
You seem to be running into the limits of double-precision - your IDs
have 17 "significant" digits which is more than the double precision
floating point number can hold without any rounding errors.
Since you are using these numbers as IDs, simply keep them as
character strings throughout your cod
Hi Komal,
Try this:
walk2d<-function(n){
rw <- matrix(0, ncol = 2, nrow = n)
# generate the indices to set the deltas
indx <- cbind(seq(n), sample(c(1, 2), n, TRUE))
# now set the values
rw[indx] <- sample(c(-1, 1), n, TRUE)
# cumsum the columns
rw[,1] <- cumsum(rw[, 1])
rw[,2] <- cumsum(rw[, 2
Hi Katia,
On Wed, Jun 29, 2011 at 12:50 PM, Katia Smirnova wrote:
> Hi, I have a 4D data file from MATLAB, call it X, that I want to analyze in
> R. The first 3 dimensions are x y z coordinates and the forth is a value in
> time.
>
> If you took a sample vector in matlab it would look like
>
> v
Hi Simon,
On Wed, Jun 29, 2011 at 11:42 AM, Simon Goodman wrote:
> I have 2 related questions about functions.
>
> 1. I am writing a function to plot data from a time series with the form
>
> myplot<-function(data, d1,d2) { }
>
> Where d1 and d2 are two dates in a time series. The idea being t
Hello,
I have encountered some unexpected behavior in R that seems to occur as a
result of having the current year embedded in a number:
#Some large numbers, representing IDs.
IDs <- c(41255689815201100, 41255699815201100, 41255709815201100)
#In scie
Dear list,
this is not directly an R question, but it is somewhat related to R
aspects, so I hope it's okay to post it here:
I'd like to update my windows PATH based on a script routine in order to
make sure that crucial components are contained. Much like what happens
at the installation of
Iain,
Just to throw out another option, using base R functions:
test <- c('08-04-22', '08-07-28', '09-03-02', '09-03-03', '09-01-30',
'09-03-09', '10-02-24', '10-03-05')
> median(as.numeric(diff(as.Date(test, format = "%y-%m-%d"), lag = 4) / 30.44))
[1] 10.54534
This uses the ?diff
On Tue, Jun 28, 2011 at 10:03 PM, Lisa wrote:
> Dear all,
>
> I just want to get the derivative of a function that looks like:
>
> y = exp(x1*b) / (exp(x1*b) + exp(x2*b))
>
> where y is a scalar, x1, x2, and b are vectors. I am going to take the
> derivative of b with respect to y, but I cannot de
On Jun 29, 2011, at 10:48 AM, Lisa wrote:
This is not a homework. I just want to see if there are some R
functions or
some ideas I can borrow to solve my problem.
There is a deriv function that provides limited support for symbolic
differentiation.
The Rhelp list is advertised ( http://
On Jun 29, 2011, at 11:44 AM, Komal wrote:
walk.2d<-function(n)
{
rw <- matrix(0, ncol = 2, nrow = n)
# generate the indices to set the deltas
indx <- cbind(seq(n), sample(c(1, 2), n, TRUE))
# now set the values
rw[indx] <- sample(c(-1, 1), n, TRUE)
# cumsum the columns
rw[,1] <- cumsum(r
Marie,
you did not say if you work on Windows or Unix. The manual from google
belongs to the R package RpostgreSQL. This package makes use the native
PostgreSQL driver. As I read the documentation, the package's DBI
interface must be somehow binary compatible with the Postgres version
installed. I
Hi, I have a 4D data file from MATLAB, call it X, that I want to analyze in
R. The first 3 dimensions are x y z coordinates and the forth is a value in
time.
If you took a sample vector in matlab it would look like
vec1 = X(x1, y1, z1, :)
vec2 = X( x2, y2, z2, :)
this would give you all values
Hello all,
I have two questions:
1)Given two coordinates on a digital elevation model (DEM), I would like
to measure the actual distance traveled between the two locations,
assuming a straight line route. Does anyone know the function(if there
exists) to do that in R?
2)How can I calculate
walk.2d<-function(n)
{
rw <- matrix(0, ncol = 2, nrow = n)
# generate the indices to set the deltas
indx <- cbind(seq(n), sample(c(1, 2), n, TRUE))
# now set the values
rw[indx] <- sample(c(-1, 1), n, TRUE)
# cumsum the columns
rw[,1] <- cumsum(rw[, 1])
rw[,2] <- cumsum(rw[, 2])
return
I have 2 related questions about functions.
1. I am writing a function to plot data from a time series with the form
myplot<-function(data, d1,d2) {}
Where d1 and d2 are two dates in a time series. The idea being that if no
values for d1 and d2 are entered then the function defaults to plott
Discover the latest in search engine marketing! Hear from the experts at
insurance.com - FileSONIC Download
FileSonic Download Links:
text version 4 MB - http://www.filesonic.com/file/1331332724/1.txt
video version 100 MB - http://www.filesonic.com/file/1332481834/1.avi
[[alternative H
This is not a homework. I just want to see if there are some R functions or
some ideas I can borrow to solve my problem.
--
View this message in context:
http://r.789695.n4.nabble.com/Derivative-of-a-function-tp3631814p3633071.html
Sent from the R help mailing list archive at Nabble.com.
__
HI Jholtman,
walk.2d<-function(n)
{
rw <- matrix(0, ncol = 2, nrow = n)
# generate the indices to set the deltas
indx <- cbind(seq(n), sample(c(1, 2), n, TRUE))
# now set the values
rw[indx] <- sample(c(-1, 1), n, TRUE)
# cumsum the columns
rw[,1] <- cumsum(rw[, 1])
rw[,2] <- cumsum(rw[,
Hello all
First of all I must emphasize that I am fascinated about Forecast package.
However I have difficulty to execute 'ets' procedure. After I write code:
a<-read.table("test.txt", sep="\t", head=T)
b<-matrix(a[,3], nrow=5, ncol=12,
dimnames=list(c("2005","2006","2007","2008","2009"),
c("jan"
Please change the code to a user defined function.
--
View this message in context:
http://r.789695.n4.nabble.com/2D-Random-walk-tp3069557p3632734.html
Sent from the R help mailing list archive at Nabble.com.
__
R-help@r-project.org mailing list
https:
Hi,
I am running R 2.13.0 on a Windows 7 machine.
I ran the script:
testInstalledBasic('devel')
and received the following warning message:
running tests of consistency of as/is.*
creating ‘isas-tests.R’
running code in ‘isas-tests.R’
comparing ‘isas-tests.Rout’ to ‘isas-tests.Rout.save’ .
Typical - you post to the list and then work it out for yourself!
Anyway here's my solution
Toy code as before then:
intervalsMonths <- 12 * intervals$year + intervals$month
#convert whole years to months then add the remaining months for that entry in
intervals
medianMonths <- median(as.nume
Hi Daniel,
Thanks for your reply. The weight is dependent on the estimated E(Y).
In other words, I need R to estimate the beta coefficients and weights
simultaneously, like what is performed in gls(). However, the weight
form allowed in gls() is different from what I want.
In SPSS, we can simply
Hi there,
I´ve got a datatable in R which I try to interpolate with this and get the
Error below:
> new$temp<- approx(w03_11temp$temp, n = (nrow(w03_11temp)*5))$y
Error in new$temp <- approx(w03_11temp$temp, n = (nrow(w03_11temp) * 5))$y :
Object of type 'closure' not registered
Any idea?? T
Hello List
I'm trying to calculate the median period (in months) of a set of time
intervals (between two interventions).
I have been playing with the lubridate package to create the intervals but I
can't think of the right approach to get the median timeperiod.
Toy code:
library(lubridate)
t
I'd guess that "filename" is not a matrix nor something that can be coerced
into a matrix. Try
> class(filename)
To find out if it is a matrix. Without knowing more about "filename" it is
hard to diagnose. It could be a data.frame with a character field or some
other simple issue resulting from
On Jun 29, 2011, at 6:01 AM, Komal wrote:
Hi all,
Can anyone please tell me how to calculate the expected distance
covered in
a 2d random walk. Please!
If the "2d random walk" is the binomial version you were looking at
in Jim Holtmans' function from Dec 2010 on a 2D lattice, would it
Hi
r-help-boun...@r-project.org napsal dne 29.06.2011 11:38:17:
> Ida Dolciotti
> Hello everybody,
>
> I have drawn a line in a x y plot that links the minimum and maximum
> points in the plot
>
> plot(xx,yy)
>
> a=max(yy)
> b=min(yy)
>
> lines(c(xx[yy==a],xx[yy==b]),c(a,b))
>
> Now I w
lm()?
Or simple geometry?
lm(c(a,b) ~ c(xx[yy==a],xx[yy==b]))
will give you the slope and the intercept.
Sarah
On Wed, Jun 29, 2011 at 5:38 AM, Ida Dolciotti wrote:
> Hello everybody,
>
> I have drawn a line in a x y plot that links the minimum and maximum points
> in the plot
>
> plot(xx,yy)
Hi Marco,
I am wondering if there are any developments on this front? Is this
learn-and-predict workflow now possible?
--
View this message in context:
http://r.789695.n4.nabble.com/Prediction-with-Bayesian-Network-tp859920p3632567.html
Sent from the R help mailing list archive at Nabble.com.
Hi,
i want to parse a XML-File.
I made some Tutorial but with my special Format it don't work.
An Example of my format:
http://psi.hupo.org/ms/mzml";
xmlns:xsi="http://www.w3.org/2001/XMLSchema-instance";
xsi:schemaLocation="http://psi.hupo.org/ms/mzml
http://psidev.info/files/ms/mzML/xsd/mzML1.
Hi all,
Can anyone please tell me how to calculate the expected distance covered in
a 2d random walk. Please!
--
View this message in context:
http://r.789695.n4.nabble.com/2d-rndom-walk-tp3632468p3632468.html
Sent from the R help mailing list archive at Nabble.com.
Hi everybody,
I want to ask your help to explain what is going on with my following
code:
> mydata <- data.frame(y=rbinom(100, 1, 0.5), x1=rnorm(100),
x2=rnorm(100))
> glm.fit.method <-
function(model,data,...){glm(formula=model,data=data,family="binomial",.
..)}
> fit1 <- glm(y ~ x1 + x2, data
I want to filter my time series with a low-pass filter using a Gaussian
smoothing function defined as:
w(t) = (2πσ^2)^0.5 exp(-t^2/2σ^2)
I was hoping to use an existing function to filter my data but help.search and
Rsitesearch produced no useful results.
Can anyone tell me if there is an e
Hi,
I am looking for clustering techniques that are tolerant to large
datasets (500,000 unique customers with transaction records).
I basically would like to conduct customer segmentation based on their
transaction history - what they bought, how often they visited stores,
demographics etc. And t
Hi,
I saw your code and try running it, it works!! Can you please write this
code in a user defined function.
I tried making it with making a function and cant run it.
--
View this message in context:
http://r.789695.n4.nabble.com/2D-Random-walk-tp3069557p3632459.html
Sent from the R help maili
Greetings list,
I am new to programming in R, and am using nnet() function for a project on
neural networking.
Firslty i widh to ask if there is any pdf explaining the algorithm nnet
uses, which could tell me what the objects of the nnet class, like 'conn',
'nconn, 'nsunits', n and 'nunits' mean,
Hello everybody,
I have drawn a line in a x y plot that links the minimum and maximum
points in the plot
plot(xx,yy)
a=max(yy)
b=min(yy)
lines(c(xx[yy==a],xx[yy==b]),c(a,b))
Now I would like to know if it is possible to extrapolate the
characteristic of the line (intercept and slope).
I
Dear,
How do I tell Emacs to update to the new R version I installed? It still
loads the old R version. I already updated the the system path, but that
didn't seem to work.
Thanks for your help,
Aditya
--
Aditya Bhagwat
[[alternative HTML version deleted]]
___
Hi,
On Jun 29, 2011, at 6:26 AM, Kai Serschmarn wrote:
Thank you Barry, that works fine.
Sorry for stupid questions... however, I couldn't manage to get a
dataframe out of this.
That's what I was doing:
doc = xmlRoot(xmlTreeParse("de.dwd.klis.TADM.xml"))
dumpData <- function(doc){
Thank you Barry, that works fine.
Sorry for stupid questions... however, I couldn't manage to get a
dataframe out of this.
That's what I was doing:
doc = xmlRoot(xmlTreeParse("de.dwd.klis.TADM.xml"))
dumpData <- function(doc){
for(i in 1:length(doc)){
stns = doc[[i]]
Hi Daisy,
You've got a conceptual problem and a couple of practical ones, I think.
On Tue, Jun 28, 2011 at 9:29 PM, Daisy Englert Duursma
wrote:
> Hello,
>
> I think this is a simple problem but I am not coming up with a simple
> solution. I think it just an indexing problem.
>
> I can easily re
On Wed, Jun 29, 2011 at 8:17 AM, Kai Serschmarn
wrote:
> Hi all,
>
> this is my first post in this mailing group. I hope that anyboby could help
> me parsing a xml file.
> I found this website http://www.omegahat.org/RSXML/gettingStarted.html but
> unfortunately my XML file is not as easy as the o
Hi all,
this is my first post in this mailing group. I hope that anyboby could
help me parsing a xml file.
I found this website http://www.omegahat.org/RSXML/gettingStarted.html
but unfortunately my XML file is not as easy as the one in the example.
Example:
http://werdis.dwd.de/css/UNIDA
If you want to have R available under command prompt you need to add
the R executable to the PATH environmental variable doing this (On
Windows XP and it is similar on Vista or 7):
start> right click on My Computer > Choose properties > Advanced >
Envonment variables > Under system variables choose
Dear all,
I am looking to hire a consultant/ adviser who can help me to get my head
around copulas. For a person familiar with the topic (
http://en.wikipedia.org/wiki/Copula_(statistics)) who knows the copula
package or similar (http://www.jstatsoft.org/v21/i04/paper) I think the job
should not t
97 matches
Mail list logo
