Displaying 20 results from an estimated 1000 matches similar to: "wrong fill colors in polygon-map"
2013 Jan 04
1
SpatialPolygon with the max value gets no color assigned in spplot function when using "at" parameter
Hi,
I would like to do coloring of map regions based on the region values
"weight". The approach I am taking is first to break regions into equal
intervals,
classIntervals(spdf$weight,4)$brks #4 intervals in this case
and coloring all regions within the interval with the same color
col = brewer.pal(4,"RdYlGn"))
The max "weight" is as well the boundary of the
2019 Mar 03
2
bug: sample( x, size, replace = TRUE, prob= skewed.probs) produces uniform sample
When `length( skewed.probs ) > 200' uniform samples are generated in R-devel.
R-3.5.1 behaves as expected.
`epsilon` can be a lot bigger than illustrated and still the uniform distribution is produced.
Chuck
> set.seed(123)
>
> epsilon <- 1e-10
>
> ## uniform to 200 then small
> p200 <- prop.table( rep( c(1, epsilon), c(200, 999-200)))
> ## uniform to 201
2010 Dec 06
10
R crashes when making graphs
hi,
i am running R with JGR and Deducer under Ubuntu Mint Lucid. since i updated to R 2.12 i can no longer make graphs - R just crashes.
the full error report was posted to launchpad (https://bugs.launchpad.net/ubuntu/+source/openjdk-6/+bug/675905) without any success so far.
anyone with a similar problem or helpful suggestions?
thanks!
kat
2010 Sep 16
2
use same breaks and colors, but the displayed scale are different-image.plot()
Hi all,
I want to put several figures in a one figure for easy comparison, so i
need to use the same methods to plot these figures. The following is an
example. I also list my method, but it does not work.
#Example data
x<- 1:10; y<- 1:10; z<- outer( x,y,"+");z2<- outer( x,y,"-")
#Quick view them
image.plot(x,y,z) #relatively larger value
image.plot(x,y,z2)
2010 Nov 17
2
Problem with plotting diagnostics - Error in object$coefficients : $ operator is invalid for atomic vectors
hi all,
after fitting a multiple linear regression
model <- lm(y ~ a + b+ c+d)
i wanted to plot diagnostics
plot(model)
but get the error message
Error in object$coefficients : $ operator is invalid for atomic vectors.
which does not make a lot of sense, since there is no "$" - i am working with an attached dataset.
can anyone help, please??
thanks a lot,
kat
2009 May 19
3
Remove objects names like character String
Hi,
how can I use rm() on objects named like:
paste("site",i,"_data",sep="") while looping
through i?
I tried rm(paste("site",i,"_data",sep="")) but I get the error that
rm() must contain names or
text strings which is confusing me as I thought paste() would create
something like that...?
Thanks,
Katharina
--
Time flies
2007 Sep 27
5
New R website: R-Cookbook.com
R Community,
I've put together a website that I thought this mailing list might be
interested in: http://www.r-cookbook.com
It's a (free) community-driven content management system for R
"recipes", or working examples. Some of the features of the site are
code highlighting, recipe ratings, recipe comments, personal "recipe
boxes" to save your favorite
2012 Oct 17
2
loop of quartile groups
Greetings R users,
My goal is to generate quartile groups of each variable in my data set. I
would like each experiment to have its designated group added as a
subsequent column. I can accomplish this individually with the following
code:
brks <- with(data_variables,
cut2(var2, g=4))
#I don't want the actual numbers, I need a numbered group
data$test1=factor(brks,
2006 Dec 08
5
Remove " from a string
Hi all!
I have lots of functions called in the following pattern
'NameOfFunctionNumber' where the name always stays the same and the number
varies from 1 to 98.
Another function which I run in advance returns the number of the function
which has to be called next.
Now I want to combine 'NameOfFunction' with the 'Number' returned so that i
can call the desired function.
I
2009 Aug 12
1
CCF for hourly time series?
Hello,
I have a dataframe containing various time series (not time series objects though!)with hourly time steps. I?d like to perform ccf for I need to know the correlation factors for different lags.
Here is an example:
x<-as.POSIXct(c("2008-12-25 16:00:00", "2008-12-25 17:00:00", "2008-12-25 18:00:00", "2008-12-25 19:00:00", "2008-12-25
2009 Jun 17
3
lattice: axis ticks, axis alignment and remove axis from plot
Hi there,
I'm a bit confused concerning the axis tck setting in the lattice
package as the ticks on left sided axis aren't
drawn at all with the following setting:
dados <- data.frame(varsep = factor(rep(1:2,10)),
i = runif(20))
library(lattice)
my.theme <- list(
axis.components = list(left = list(tck = 1, pad1 = 1, pad2 =
2), top = list(tck = 0, pad1
2009 Jun 24
3
distinguish regression lines in grouped, black and white lattice xyplot
Hi,
I've got the following problem which I cannot think of a solution right now:
if got a lattice xyplot in black and white and a grouping variable
with many (more than 8
values) and I plot it as regression lines (type="r"), just like this
one (not reproducable but that's
I guess not the point here):
xyplot(log(AGWB) ~ log(BM_roots), data=sub_agwb_data, groups=species,
2010 Dec 03
3
recode letters to numbers
Hello,
I have a dataframe assigning various scores on around 20 variables to a list of
countries. The scores are rated on a scale of (D, C, B, A) and there are also
some not rated ones (NR) and others are left blank (NA). I now wanted to
transfer the scores into numeric values (such as NR=0, D=25, C=50, B=27, A=100),
while also leaving the NA in the dataframe. Can I transform the whole
2010 Mar 18
5
language!!!
Hi,
there´s something that really bothers me about R and after hours and hours of internet research, I´m still stuck with the same problem: I installed R and it is in Spanish, as the system on my work. I would really like it in English and there seems to be no way to change that!?! I´m doing an abroad internship, so my Spanish is far from perfect and my tutorial is in English.
Is there any
2011 Nov 29
1
Hmisc break points error
I am making frequency histograms using the histbackback function on my 2
datasets. However when I try to use the brks function:
foo<-histbackback(log(fie11),
log(fie86),ylim=c(0,9),probability=FALSE,axes=TRUE,ylab=("log10 Parcel
Size"),brks=16)
The graphic results in a 'NA' label for the y axis (no intervals are
returned)
Also when I use 'summary(foo)' the
2009 May 19
2
create string of comma-separated content of vector
Hi,
how do I create a string of the comma-separated content of a vector?
I've got the vector i with several numeric values as content:
>str(i)
num 99
and want to create a SQL statement to look like the following where
the part '(2, 4, 6, 7)' should be
the content of the vector i:
select * from [biomass_data$] where site_no in (2, 4, 6, 7)
Here my approach (which doesn't
2009 Aug 08
1
linear model: Test difference between coefficients and given values (t.test?)
Hi there,
I've got a question which is really trivial for sure but still I have
to ask as I'm not
making any progress solving it by myself (please be patient with an
undergraduate
student):
I've got a linear model (lm and lmer fitted with method="ML").
Now I want to compare the coefficients (slope, intercept, not the
random effects)
of both models with a given value (e.g.
2009 Jun 26
1
predicted values after fitting gamma2 function
Question: after fitting a gamma function to some data, how do I get
predicted values? I'm a SAS programmer, I new R, and am having
problems getting my brain to function with the concept of "object as
class ...". The following is specifics of what I am doing:
I'm trying to determine the pdf from data I have created in a
simulation.
I have generated frequency counts
2012 Jun 06
2
package zoo, function na.spline with option maxgap -> Error: attempt to apply non-function?
Hello,
I'm trying to use na.spline (package zoo) to fill some missing data in a time series.
this works fine, however, if I apply the 'maxgap' argument, I always get the error:
<------
Error in na.spline.vec(x., coredata(object.), xout = xout., ...) : attempt to apply non-function
------>
I couldn't find a similar error for this case in the mailing lists and zoo vignette,
2007 Nov 12
2
separate window for help/newbie
hi,
is there a way to make r open the help file in a separate window?
i'm quite new to linux(kubuntu) and used to work with r under windows...
thanks,
katharina