Displaying 20 results from an estimated 1000 matches similar to: "using ifelse with surprising results"
2009 Apr 16
2
Reading in a large number of dbf files
good morning
This question is not a stats question per say but a data management and
lattice plotting problem. I apologize now if I'm asking an inappropriate
question to this gracious group.
I'm need to bring in approximately 100 *.dbf files into R but I'm having
difficultly understanding several examples I've tracked down regarding this
procedure and could benefit from your
2009 Feb 02
2
Upgrading to TINN - R 2.1.1.6
I know this has been addressed before, but I'm still confused by the
solution.
When running TINN - R (ver 2.1.1.6) with R (2.8.0) I have lost the
functionality of TINN-R
It returns an error
source(.trPaths[5], echo=TRUE, max.deparse.length=150)
I modified the RProfile.site file as follows:
options(pager="interal")
.trPaths=c(
'C:/Documents and
2009 Dec 08
2
Upgrading To 2.10 from 2.6.2
Hello
I have a Linux machine (Ubuntu 8.04 hardy, Gcc version 4.2.4
(i486-linux-gnu) currently running R 2.6.2. I'd like to upgrade to 2.10.
First Question): What is the appropriate way to remove the old version of
R?
Part 2.
After downloading r-base_2.10.0.orig.tar.gz and opening the archive. I
ran the ./configure routine.
It failed claiming that it could not find the F77 compiler.
My
2011 Nov 03
3
Plotting skewed normal distribution with a bar plot
Hi,
I need to create a plot (type = "h") and then overlay a skewed-normal
curve on this distribution, but I'm not finding a procedure to accomplish
this. I want to use the plot function here in order to control the bin
distributions.
I have explored the sn library and found the dsn function. dsn uses known
location, scaling and shape parameters associated with a given input
2009 Mar 31
3
labeling panels in lattice plots
I am using windows XP with R 2.8.1
I am generating a lattice plot of annual rain patterns using the following
function:
> xyplot(rain.stats$min+ rain.stats$max + rain.stats$ave ~
rain.stats$month |rain.stats$year,
lty = 1, data = rain.stats, type = c("l","l", "l"), col =
c("red", "blue", "green"), distribute.type =
2010 Oct 19
2
Clustering with ordinal data
Hello
I've been asked to help evaluate a vegetation data set, specifically to
examine it for community similarity. The initial problem I see is that the
data is ordinal. At best this only captures a relative ranking of
abundance and ordinal ranks are assigned after data collection. I've
been trying to find a procedure in R that can handle ordinal based
classification and so far have
2010 Jul 09
2
Ctree Question
Hello,
I've been using ctree and have developed a 55 node - 28 terminal solution.
As can be imagined, the plot is difficult to travel down each of the major
branches.
I've read the help files for ctree I saw where terminal nodes can be color
coded.
plot(airct, type = "simple")
> plot(airct, terminal_panel = node_boxplot(airct, col = "blue", + fill =
hsv(2/3,
2010 Feb 24
2
mlogit is not an S4 object error
Hello,
I've been getting the following error when using the mlogit function from
the mlogit package
This is one of the examples provided in the Package "mlogit" January 27,
2010 description
data("Fishing", package="mlogit")
Fish <- mlogit.data(Fishing, varying = c(4:11), shape="wide",
choice="mode")
summary(mlogit(mode ~ pr + ca - 1,
2010 May 10
2
Installing randomForest on Ubuntu Errors
Hello,
I've tried to install randomForest on a Ubuntu 8.04 Hardy Heron system.
I've repeatedly rec'd the error:
> install.packages("randomForest", dependencies = TRUE)
ERROR: compiliation failed for package 'randomForest'
** Removing '/home/admuser/R/i486-pc-linux-gnu-library/2.6/randomForest'
The downloaded packages are in
2009 May 18
2
Superscripts and subscripts in trellis graphics
I'd like to annotate an xyplot with an R^2 value, but can not find the
syntax to define the superscript format for the text.
I'd appreciate suggestions,
xyplot(SharkSloughEggs.df$Sharkeggs.rel +
SharkSloughEggs.df$SharkHatched.rel + SharkSloughEggs.df$SharkFlooded.rel +
HSI.shark$MEAN ~ Year,
pch=c(1,9,5,4), lty=c(1,6,3,4), lwd = c(2), col = c("black",
2010 Jun 21
2
ctree
Hello,
This is a re-submittal of question I submitted last week, but haven't rec'd
any responses.
I need to extract the probabilities used to construct the barplots
displayed as part of the graph produced by plot("ctree").
For example,
library(party)
iris.ct <- ctree(Species ~ . , data = iris)
plot(iris.ct)
Instead of a simple example with only 4 terminal nodes, my
2009 Apr 28
3
help working with date values
My data contains a variable "observation_date" and it contains values as:
1985-09-02
1985-09-15
1985-07-31
1985-09-02
I need to process data annually rather than daily, therefore I'm trying to
1) either extract the first 4 digits from this field and use them as a new
variable "year" or 2) keep the variable as it is and process the analysis
using the first 4 digits of
2009 Jan 30
2
help with time series
Hello everyone
I'm working with R 2.8.1 on a windows machine
I have a question regarding time series analysis
The first question is how does R expect the input file to be structured?
I'm working with a *.txt file similar to the abbreviated one here:
Date,stage
4/2/1953,7.56
4/3/1953,7.56
4/4/1953,7.54
4/5/1953,7.53
4/6/1953,7.5
4/7/1953,7.47
4/8/1953,7.44
4/9/1953,7.41
4/10/1953,7.37
2009 Apr 09
2
numbers not behaving as expected
If someone can explain this odd behavior I'd appreciate it.
I have a data.frame Cell.ave (attached and created via dput(Cell.ave,
"Cell.ave")
which contains three columns of parameters year, month and AveRain.
I need to subset the data.frame by months such that
DrySeaonMonths are 1,2,3,4, 11, and 12 or Jan - April, November and
December
and
WetSeasonMonths are 5:10 or
2008 Sep 18
2
graphing netCDF files
Hello
I'm working with a large hydrological data set stored in a netCDF format.
The file stores x and y coordinates in the UTM projected coordinate system,
yet when I use image to graphically display the z variable, the image is
distorted in the sense that it does not plot the map in the correct spatial
organization.
I'm wondering if I need to define the projection of the netCDF file
2012 Apr 23
2
automating a script to read a file
Hi,
The following script (which I did not develop) is used to calculate and
plot a skewed normal curve. The script currently requires the user to
input six parameters, rather than reading these directly from a file.
I've been spinning wheels here, trying to figure out how to modify the
script to automate it. I have four data sets, each in excess of 300
records that I need to process.
My
2010 Jun 25
1
variograms and kriging
Hello
Trying to develop variograms and kriged surfaces from a point file. Here is
what I've done so far.
library(gstat) # also loads library(sp)
library(lattice)
soilpts$x <- soilpts$UTM_X
soilpts$y <- soilpts$UTM_Y
soil.dat <- subset(soilpts, select=c(x, y, Area, BulkDensity, LOI, TP, TN,
TC, Total_Mg))
dim(soil.dat)
[1] 1292 7
coordinates(soil.dat) <- ~ x+y
2008 Nov 13
1
Plotting from a list
Hello,
I am working with version 2.7.2 on a PC and have 2.8.0 available, but have
not upgraded completely yet.
I have 6 species for which I've run 6 unique environmental management
scenarios against for comparison purposes. Each scenario is run for 36
years (1965 - 2000). Some of the species have subpopulations and other do
not. For those with subpopulations (A, B, C, D, E, and F) I have
2009 Mar 24
1
Help with lattice
Hello.
I'm working on a windows XP machine with R 2.8.1
I'm working with lattice and I don't understand how to set the size of the
text for the labels. Can anyone suggest a solution?
I'm establishing a melted dataframe first and then setting a key and
calling the plot as follows;
> key.variable <- list(space = "right", text =
list(levels(wbirdm$variable)),
2009 Apr 23
1
Setting lattice par parameters
Hello
I'm plotting a large suite of barcharts and need to modify the size of the
text for both the yaxis and xaxis labels.
I've tried using the following:
> trellis.par.set(list(par.ylab.text = list(cex = 0.65)),
trellis.par.set(list = par.xlab.text = list(cex = 0.65))))
On inspection, however after I invoke this line,
> trellis.par.get("par.ylab.text")
>