Displaying 20 results from an estimated 1100 matches similar to: "local average"
2004 Mar 09
1
vector extraction
Hello,
I could need some help on this one:
>From the data.frame "Test.dataset2" below (TSCS data for 151
"countries.to.map" for "year" 1973-95; each "country.to.map" is described by
a unique code), I would like to extract a vector "color" that for each
"country.to.map" takes on the value of "dv" (a categorical variable with
2005 Nov 23
3
Infinite recursion in S3 methods crashes R on windows (related to PR#8203?)
Hi,
Infinite recursion in S3 methods seem to crash R on Windows 2000 (R
terminating with the ("Rgui.exe has generated errors...") message,
rather than throwing an error. This happens with both Rgui and Rterm.
The following toy example triggers this:
myf <- function(x, ...)
UseMethod("myf")
myf.default <- function(x, ...)
myf(x)
myf(1)
...R crashes...
Which I
2005 Apr 05
1
summing columns using partial labels
I have a dataset of the form
Year tosk.fai tosk.isd tosk.gr ....... tosk.total hysa.fai
hysa.isd ...
and so on. I want to sum all the columns using the first four letters in
the columns label(e.g. 'tosk', 'hysa' etc.). How can you do that? Also,
the sums should be without the '.total'column (e.g. 'tosk.total') as
this serves as a check that
2005 Apr 20
2
Suggestions for manipulating formula objects
I'm trying to manipulate/change a formula prior to passing it to another
function. A simplified example:
User passes formula to my function: y~x
My function does: lm(transform(y)~x)
Here, transform() is added to the model's response.
What is the best way to accomplish this?
2005 Apr 06
1
read.table with header and text data
Hi R!
I am reading in a text file which has one column of alpha data and 5 columns
of numeric data.
There is a header row.
I would like the alpha data column to just be character rather than factor.
Is there a way to do this, please? I'm thinking that it might be I() but
can't figure out exactly how.
Thanks,
Laura
mailto: lauraholt_983 at hotmail.com
R 2.0.1 Windows
2005 Apr 13
1
how to separate a string
hello,
i wonder how is string represent in R. if i have a string s= "hello", how can i refer to first character in the string s?
also if i have s1 = "hello.1", s2 = "ok.1", how can i separate the s1 into "hello" "1" and s2 into "ok" and "1"? I have tried to use the substring function, but i don't where i can get the index
2005 Apr 13
1
i param in "for" loop does not takes zeros?
Hi all
Is there any reason why the parameter i in a "for" loop ignores a value of
zero? For example
sim=c()
p=.2
for(i in 0:5)
{sim[i]=dbinom(i,5,p)
}
sim
[1] 0.40960 0.20480 0.05120 0.00640 0.00032
In this example the quantile i= 0 was ignored since
dbinom(0,5,p)
[1] 0.32768
The same behaviour occurs if I use a while loop to perform the same
calculation:
sim=c()
p=.2
i=0
2005 May 04
1
Questions about the intersection area under two kernel densities
Hi there,
I am working on a project which needs the value of the interaction area under two distributions( eatimated by kernel density estimators).
For example:
x<-rnorm(100,0,1)
y<-rnorm(100,0.2,1)
density(x) # This produces the summary of dependent variable and independent variable.
How can I get the individual values of variables and reform a curve to calculate the area under
2005 Apr 06
2
HP-UX and IRIX recent builds? Any other rare platforms?
I am revising the section in R-admin on platforms. When I asked about
locales earlier this year I got no reply about HP-UX and IRIX, so I will
presume that no one has built R on those OSes recently. If you have,
please send me the OS number and the flags you used.
Platforms where I have seen recent reports:
Linux
MacOS X
Solaris 8, 9, 10 (any Solaris 7 users please let me know)
AIX 5.1
2005 Apr 07
2
vectorized approach to cumulative sampling
Hi All,
I need to sample a vector ("old"), with replacement, up to the point
where my vector of samples ("new") sums to a predefined value
("target"), shortening the last sample if necessary so that the total
sum ("newsum") of the samples matches the predefined value.
While I can easily do this with a "while" loop (see below for example
2005 Apr 04
3
plotting mathematical notation and values substitution
Dear R-users,
I'm trying to add a title on a plot with both mathematical notation and
values substitution. I read the documentation and search the mailing list
but I was not able to solve my problem. Actually, there is a message by
Uwe Ligges on June 2003 which addresses a question very close to mine, but
the code provided doesn't work. The code is the following:
# I add this to let you
2005 Apr 19
1
controlling the x axis of boxplots
v 2.0.1 (sooooh old!) on Win2k
I think I know the answer to this but I can hope ...
I have data for continuous variables (measures of residents) by a
categorical variable in range (1,22), the units in which they live.
I want to plot these data with a pair of boxplots one above another
with same x-axis (1,22) using par(mfrow=c(2,1)) and then plotting
first for the women then for the men.
2005 Apr 13
3
terminate R program when trying to access out-of-bounds array element?
I want R to stop running a script (after printing an
error message) when an array subscript larger than the
length of the array is used, for example
x = c(1)
print(x[2])
rather than printing NA, since trying to access such
an element may indicate an error in my program. Is
there a way to get this behavior in R? Explicit
testing with the is.na() function everywhere does not
seem like a good
2005 Apr 10
2
Plotting the occassional second label
Dear useRs,
I'm trying to plot spray quantities against dates, and label the points
on the plot. Basically quite simple, but sometimes two chemicals have
been used and are listed in separate rows in the table as below; then
the labels are written on top of each other.
> spray
SprayDate PD Trt Qwater Qai
1 2005-03-09 Spidermites Pegasus 1300 1.04
2
2005 Apr 05
3
How to do aggregate operations with non-scalar functions
Hi,
I have a data set, the structure of which is something like this:
> a <- rep(c("a", "b"), c(6,6))
> x <- rep(c("x", "y", "z"), c(4,4,4))
> df <- data.frame(a=a, x=x, r=rnorm(12))
The true data set has >1 million rows. The factors "a" and "x"
have about 70 levels each; combined together they subset
2005 Apr 04
1
need any advises for code optimization.
Dear colleagues,
I have the following code. This code is to 'filter' the data set.
It works on the data frame 'whole' with four numeric columns: a,b,d, and c.
Every row in the data frame is considered as a point in 3-D space.
Variables a,b, and d are the point's coordinates, and c is its value.
This code looks at every point, builds a cube 'centered' at this
2005 May 11
2
time zones, daylight saving etc.
Hi, I have a whole bunch of data, which looks like:
15/03/2003 10:20 1
15/03/2003 10:21 0
15/03/2003 12:02 0
16/03/2003 06:10 0
16/03/2003 06:20 0.5
16/03/2003 06:30 0
16/03/2003 06:40 0
16/03/2003 06:50 0
18/03/2003 20:10 0.5
etc. (times given on a 24 hour clock)
and goes on for years. I have some code:
2005 Oct 13
3
Optim with two constraints
Hi R-list,
I am new to optimization in R and would appreciate help on the following
question. I would like to minimize the following function using two
constraints:
######
fn <- function(par,H,F){
fval <- 0.5 * t(par) %*% H %*% par + F%*% par
fval
}
# matrix H is (n by k)
# matrix F is (n by 1)
# par is a (n by 1) set of weights
# I need two constraints:
# 1.
2005 Apr 11
0
S4 group methods don't dispatch for "Summary"?
Hi all,
I am having a little problem with S4 group generics (apologies if I
get some terminology wrong below). I'm finding I can set methods for some
group generic functions, but not others:
> setClass("foo", representation(data="numeric"))
> x <- new("foo", data=1:10)
## Setting "Ops" works just as expected:
> setMethod("Ops",
2004 Mar 08
1
drawing filled countries according to data using map('world')? - follow up
Hello,
this is a follow up on my previous inquiry regarding the use of the map
library (Becker and Wilks 1993).
Using the 'world' database I would like to draw filled countries in a world
map so that the filling colors of each country corresponds to the value of a
policy variable "fix.float" at a specific "year" (the goal is to visualize a
policy diffusion pattern