similar to: Bug in stars.R (PR#739)

This message is in MIME format. The first part should be readable text, while the remaining parts are likely unreadable without MIME-aware tools. Send mail to mime@docserver.cac.washington.edu for more info. --1465796871-645734172-974620647=:18687 Content-Type: TEXT/PLAIN; charset=US-ASCII Hi all, I've discovered and (hopefully) fixed a bug in the stars function I wrote some time ago.

Hi all, I've got a density graph made with the following commands: win.graph(width=13,height=6) par ( fin=c(13,3) ,mai=c(1,1,0.5,0.5) ,mfrow=c(1,2) ,cex.axis=1.5 ,cex.lab=1.5) dens<-density(DATA1.y[2,]-mean(DATA1.y[2,]),kernel="gaussian") xlimit<-range(dens$x) ylimit<-range(dens$y) hist( DATA1.y[2,]-mean(DATA1.y[2,]) ,xlim=1.1*xlimit ,xlab=expression(q[e])

Hello r-help, I'm fitting a model with lm() and using the orthogonal polynomials from poly() as my basis: dat <- read.csv("ConsolidatedData.csv", header=TRUE) attach(dat) nrows <- 1925 Rad <- poly(Radius, 2) ntheta <- 14 Theta <- poly(T.Angle..deg., ntheta) nbeta <- 4 Beta <- poly(B.Beta..deg., nbeta) model.1 <- lm( Measurement ~ Block + Rad + Theta + Beta

R 1.4.0 ?stars has scale: logical flag: if `TRUE', the columns of the data matrix are scaled independently so that the maximum value in each column is 1 and the minimum is 0. If `FALSE', the presumption is that the data have been scaled by some other algorithm to the range [0,1]. but the code has if (scale) { x <- sweep(x, 2,

Hi, I don't know much about the RGL package, and I have read the documentation and tried some parameters, with no luck... I would like to generate a movie from a 3D object (code below), where the vortex A is closer to the observer, and then the object rotates and the B vortex gets closer. I would like to capture this movie to a file. By the way, I am not being able to insert unicode text

I thought about that star graph again, and realized that it would be quite a handy thing for visualizing cyclic data like time or compass direction. Here is a cleaned up (and renamed) version to do a polar plot that starts at the right and goes counterclockwise or a 24 hour clock plot that starts at the top and goes clockwise. There are probably other varieties that would be interesting. Jim

Dear List members I am trying to find a way to generate a random point in a circle centred in a geographical location. So far I have used the following formula (see code below): random_x = original_x + radius*cos(angle) random_y = original_y + radius*sin(angle) where radius is a random number between 0 and the radius of the circle and angle is between 0 and 360 degrees The code bellow works

I would like to set the y axis limit of an xyplot using the object 'ylimit', but receive this error: [1] 990 Error in extend.limits(limitlist[[i]], axs = axs) : improper length of lim I get the same error if I use ylim. library(lattice) trellis.device(col = FALSE, theme = lattice.getOption("col.whitebg")) name <- "Variable name" symbols <-

###Dear R users ###I have been using SensoMineR package from CRAN for most of my work in sensory data analysis and from my usage experience, I encountered some areas for improvement and considered ###modifying the function in SensoMineR package for my personal use. I felt that it could be useful to share this to the community for enabling adoption by other users where they might require a

I think there are two bugs in aov() that shows up when the right hand side of `formula' contains both `-1' and an Error() term, e.g., aov(y ~ a + b - 1 + Error(c), ...). Without `-1' or `Error()' there is no problem. I've included and example, and the source of aov() with suggested fixes below. The first bug (labeled BUG 1 below) creates an extra, empty stratum inside

Hello, I am quite new to R and although I have read some material (apparently not enough) I'm getting very frustrated. I'm trying to do some simple circular statistics in R (CircStats, primarily circ.mean) and I'm confused about the results. What is a negative radian or degree mean in R? I thought it would be a coterminal angle, but that doesn't make sense with my data

Dear R helpers, I try to use the 'delaunayn' function in the 'geometry' package for Delaunay triangulation in 2 dimensions. For the four following points, I get a warning message : > coord=matrix(ncol=2,byrow=TRUE,c(622633,7073452, + 621228,7073517, + 621879,7071762, +

I believe you are right, but can you please explain why anyone would want to fit this model? It differs only in the coding from aov(y ~ a + b + Error(c), data=test.df) and merely lumps together the top two strata. There is a much simpler fix: in the line if(intercept) nmstrata <- c("(Intercept)", nmstrata) remove the condition (and drop the empty stratum later if you

Hello, I am a post-doc of the Federal University of Santa Catarina State (UFSC). Last year, used EnQuireR for hirarchical cluster analisis and end up very well. I formated my computer couple months ago and installed R again as version x64 3.0.2. have new data which ENMCA function of EnQuireR package is not running. R seems to be fine as it runs funcions of other packages. Even MCA

I'm having trouble understanding the output from as.windrose(). For one thing, data on a boundary between sectors seem to be left out of the counts. I assume that explains the missing point in the output below (angle 45). Shouldn't one side of each sector interval be open, to include values such as my 45 in the example? Also, why does the angle 180 in my input apparently not result in

I am merely trying to reproduce Figure 1.2 of Chris Chatfield's 6th edition of his The Analysis of Time Series: An Introduction (page 2). The S-PLUS code is on pages 305-306. I am almost there but I am having a heck of a time trying to modify and change the x-axis per the book. The book shows the x-axis with 10 tick marks, correctly positioned, and labeled "Jan 53", ...,

On 21/10/2020 8:45 a.m., Balint Radics wrote: > Hello, > > Could someone suggest a package/way to make a 3D raster plot of the Earth > (with continent boundaries), and then make a "cut" or "slice" of it such > that one can also visualize some scalar quantity as a function of the > Radius/Depth across that given slice ? > > Formally, I would have a

Hello, I am using the following code to create ppp files from csv data and map shape files, but I am getting some errors which I have been unable to fix by searching them online: library(spatstat) library(maps) library(maptools) NYC2<-readShapePoly("nybb.shp") # this is a map of the NYC boroughs without waterways and no census tract divisions (but it does include lines separating

Hi R users, I would like to construct a sort hybrid vector/scatter plot. My data is in the following format: 3-column x,y,z data-frame in which every row is a separate data-point. The x & y columns are coordinates, and the z column contains orientation data (range 0-180 degrees, with East=0 & North=90). I need to set each x,y, point to have the alignment in z. Hence my

Hi all, I have currently been using the book Modern Applied Statistics with S from Venables and Ripley. At chapter 6 on Linear Statistical models I wanted to produce the plot as shown by "Figure 6.1" using Whiteside's data. xyplot command seems not to work on my version of "R"(version 1.7.0) running on different environments (Win XP, Win 98 and Win 2000). I then used the