similar to: Interpolation of azimuth values

Hello, Once again, I posted a message without a subject line. Sorry.... here is the question again. Is there a simple way to modify the circ.mean function in the CircStats package to include a vector of weights to obtain a weighted average angle? Thanks! Martin -- Martin Biuw Sea Mammal Research Unit Gatty Marine Laboratory, University of St Andrews St Andrews, Fife KY16 8PA Scotland Ph:

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

Hi all, I would like to know a way of calculating the distance between 2 points when I already have the longitude and latitude of the points. For example : Point 1 : 52? 9'54.00"N 4?25'8.40"E Point 2 : 52?27'46.80"N 4?33'18.00"E Distance between point 1 and point in km .... Is there any functions already available for this ? Any help will be much

Is there packaged code to convert geographical coordinates (e.g., longitude, latitude, elevation) to Cartesian coordinates in 3-space? I can see how to do this using 1. a spherical-to-Cartesian conversion like pracma::sph2cart(tpr) http://cran.r-project.org/web/packages/pracma/ 2. a geographical-to-spherical conversion. This seems to involve (in roughly increasing order of difficulty or

Hi all, My apologies for the confusion ... but I am re-naming this thread again because the proposal doesn't come from Ambisonia. As the founder of Ambisonia I have only been 'mediating' between some members of the ambisonic community ... so it would be wrong to associate the spec with Ambisonia. To clarify. This is what the spec proposes: - that Mapping = 1 means the contents are

On 2/28/07, Ivo Emanuel Gon?alves <justivo@gmail.com> wrote: > On 2/28/07, Ralph Giles <giles@xiph.org> wrote: > > Well, there are todo pages at wiki.xiph.org, but I meant more in the > > community folklore sense. My point is a roadmap doesn't help much unless > > there are people committed to making things happen. That's been the > > problem with a

On 2/28/07, Ivo Emanuel Gon?alves <justivo@gmail.com> wrote: > On 2/28/07, Ralph Giles <giles@xiph.org> wrote: > > Well, there are todo pages at wiki.xiph.org, but I meant more in the > > community folklore sense. My point is a roadmap doesn't help much unless > > there are people committed to making things happen. That's been the > > problem with a

Hello R-helpers, What I have: I am working with vectors not matrice: #Basic equations A <-c(-20,-9,-2) x <-c(0.17,0.22,0.61) B <- crossprod(A,x) # R matrix multiplication works with vectors A%*%x==B # Is true... Question: If x is unknown and A and B are known, how do I solve for x in R? solve(A,B) won't work because A is not a square matrix solve(A,B) Error in solve.default(A,

I''m new to Ruby, Rails, and this list. I''ve tried searching for an answer but it all seems to be over my head. I have all the books but again they seem to leave me on my own to figure out the stuff that really matters. Where is the best place to find out the details? Exactly what are these much celebrated "naming conventions"; i.e. where do I find an exhaustive

Dear all, I am using the function 'D' in the 'stats' package to perform symbolic derivation. This works very well and it is much faster than e.g. Mathematica (at least for my purposes). First, I would like to thank the development team for this excellent function. However, I run into trouble in some cases, particularly when I am to do some operations on long expressions

> On Apr 11, 2017, at 10:25 PM, Madhur Amilkanthwar via llvm-dev <llvm-dev at lists.llvm.org> wrote: > > I am interested in knowing more. > 1. What benchmarks does LLVM community use for compile-time study? I see CTMark, but is that the only one being analyzed? CTMark is not cast in stone. Its purpose is for the community to have a trackable proxy for the overall llvm test

I am examining the following nlme model. asymporig<-function(x,th1,th2)th1*(1-exp(-exp(th2)*x)) mod1<-nlme(fa20~(ah*habdiv+ad*log(d)+ads*ds+ads2*ds2+at*trout)+asymporig(da.p,th1,th2), fixed=ah+ad+ads+ads2+at+th1+th2~1, random=th1+th2~1, start=c(ah=.9124,ad=.9252,ads=.5,ads2=-.1,at=-1,th1=2.842,th2=-6.917), data=pca1.grouped) However, the two random effects (th1 and th2)

My apologies, the data can now be found at: url <- "http://www.econ.uiuc.edu/~roger/research/ebayes/velo.d" x <- scan(url,skip = 1) If I could get each of the histograms to mimic what is produced by hist(x, 100, freq = FALSE) I?ve experimented with xlim, ylim, without success so far... url: www.econ.uiuc.edu/~roger Roger Koenker email rkoenker at uiuc.edu

Hello all, My question is how do I write a multilevel regression model of individual responses to environmental predictors that accounts for social interactions between individuals. i.e.; 1) Individual response is nested within a group response. 2) The same environmental predictors apply to both hierarchical levels but, 3) Lower level slope/intercept are also affected by high-level response.

Folks, I have a dataframe (snippet shown below). > demo.df[1:10, 1:6] dirn county year exp exp.wave r3 1 43901 Cuyahoga 2006 0 0 56 2 49098 Pickaway 2006 0 0 77 3 44164 Portage 2006 0 0 85 4 44610 Wayne 2006 1 1 76 5 45120 Wayne 2006 0 0 82 6 49593 Scioto 2006 1 1 89 7 46516 Crawford 2006 0 0

Dear all, I want to solve a score function by using Newton-Raphson algorithm. Is there such a fucntion in R? I know there's one called optim, but it seems only doing minimizing or maximizing. Thanks, Jimmy

How to use "contour" function? I type "contour(th1,th2,SumofSquares,levels=c(seq(1000,4000,200)))" But I didn't see contour in the plot, I guess it is because my data is so big that contour is out of range. Does anybody know how to fix this problem? here th1 is a 1x101 array(158741.8-298529.6) th2 is another 1X101 array(-0.0058352269-0.2093140935) SumofSquares is a

I am trying to do some comparisons of density estimators using lattice. The code below attempts to plot the same histogram in each panel and then overplots a kernel estimate with different bandwidths. Finding packet.number() was a bit arduous, but seems to do what I want. My concern now is that close examination of the 4 histograms reveals that they are different even though they use the same

Hi everyone! I meet one problem when embedding R in C code, when I run the the R code in one child thread , it always print error info: Error: C stack usage is too close to the limit I also try to set R_CStackLimit = (uintptr_t)-1 to disable the C stack check as the R-exts doc say, but it still does not work, the error info still exist. Besides it is interesting that if i

Can some kind person point out my error here? I'm trying to set up a grid for a countour plot of a likelihood function. > u <- rnorm(20,9.5,2.5) > # sample of size 20 from N(9.5,2.5^2) > loglik <- function(th1,th2) { + n <- length(u) + -(n/2)*log(2*pi*th2^2)-0.5*sum((u-th1)^2/th2^2) + } > x <- seq(4.5,14.5,len=50) > y <- seq(0.5,6,len=50) > f <-