similar to: spatial logit help

Hello R users... how to align this Zmat (triangle image) in X axis? I would like that the triangle's base become in the X axis and the triangle's height become in the Y axis Is there some trick for make this? Thanks. Cleber ######################## my test and try f <- function(x,y){ z=1-x-y z[ z < (-1e-15) ] <- NA return( -100*x+0*y+100*z ) } x = seq( 1, 0,

I am having a problem with the package SSPIR. The code below illustrates it. I keep getting the message: "Error in y - f : non-conformable arrays." I tried to tweak the code below in many different ways, for example, substituting rbind for cbind, and sometimes I get a different error message, but I could not find a variation of this code that would work. Any help will be greatly

Hi all, I'd like to plot the Real and Imaginary parts of some f(z) as two different surfaces in wireframe (the row/column axes are the real and imag axes). I know I can do it by, roughly speaking, something like plotz <- expand.grid(x={range of Re(z)}, y={range of Im(z), groups=1:2) plotz$func<-c(Re(f(z),Im(f(z)) wireframe(func~x*y,data=plotz,groups=groups) But that seems like a

I am trying to write the following code in R. The code works in S+ and i am trying to do the program in R. x=discrim(admit~gpa+gmat,prior=c("uniform"),data=data.mm) i wrote the following in R: x=lda(admit~gpa+gmat,data=data.mm) i could not figure out how to write prior=c("uniform") in R. I would get an error every time. I think that it has something to do with

At the risk of ridicule for my deficient linear algebra skills, I ask for help using the solve.QP function to do portfolio optimization. I am trying to following a textbook example and need help converting the problem into the format required by solve.QP. Below is my sample code if anyone is willing to go through it. This problem will not solve because it is not set up properly. I hope I

Hi, I am plotting values of log(hazards ratio) as a function of two predictors, using the plotting function filled.contour(). Here is a simple simulated example of this: x <- seq(0,1, length=20) y <- seq(0,1, length=20) z <- outer(x,y, function(x,y) x^2 + y^2 ) zmat <- matrix( rexp(n=400, rate = z+0.001), 20, 20) filled.contour(x, y, log(zmat),

Hi all I have a question about correct usage of persp(). I have a simple neural net-based XOR example, as follows: library(nnet) xor.data <- data.frame(cbind(expand.grid(c(0,1),c(0,1)), c(0,1,1,0))) names(xor.data) <- c("x","y","o") xor.nn <- nnet(o ~ x + y, data=xor.data, linout=FALSE, size=1) # Create an (x.y) surface and predict over all points d <-

Hey All, I downloaded winebottler so that I can run the GMAT Prep software and was successful in downloading it but i am now running into a few problems. For examples, if there is a "next" icon on the window, I will need to click on it and hit the space bar simultaneously to move on to the next page. I was able to do this on the writing assessment section of the practice test but on the

Dear R-Users, I'm using SSPIR package for a spatio-temporal application. Is it possible to modify the structure of the involved matrixes (Fmat, Gmat, Vmat,Wmat)? I want to create a model like this #y(t)=k*theta(t)+epsilon(t) #theta(t)=h*theta(t-1)+eta(t) #epsilon(t) N(0,V) V=sigma2*I #eta(t) N(0,W) W=sigma2_eta where the state variable theta has dimension 1(p=1) and at

I have been trying to specify a different color palette to the image method in plot.Design. My model has crossed two rcs() arguments and one two-level gender argument. The goal which appears to have been mostly achieved is to produce separate bivariate plots for men and women The call to plot does produce a level plot but it appears only with the default color palette despite various

Hi, All: Is there a simple way to detect complex conjugates in the roots returned by 'polyroot'? The obvious comparison of each root with the complex conjugate of the next sometimes produces roundoff error, and I don't know how to bound its magnitude: (tst <- polyroot(c(1, -.6, .4))) tst[-1]-Conj(tst[-2]) [1] 3.108624e-15+2.22045e-16i

I'm currently running into a little trouble with the kfilter method, and would love some clarification if you are able to offer it. When trying to run kfilter, I've been running into errors that seem to result from having mismatched dimensions. Specifically, the dimension of my observations is 2, while the dimension of the state space is 4. In the filterstep function (file sspir_kfs.R),

Hello, I apologize for such a basic question, but I have been trying to do this in multiple packages without much success. I am trying to set up a state space model for Kalman filtering. I am using package dlm. The DLM is specified by: observation: y(t) = F(t)*theta(t) + v(t) state: theta(t) = G(t)*theta(t-1) + w(t) I have no problem setting up a simple example where F is constant. I am

Bugs item #13676, was opened at 2007-09-05 14:01 You can respond by visiting: http://rubyforge.org/tracker/?func=detail&atid=218&aid=13676&group_id=35 Category: Incorrect behavior Group: current Status: Open Resolution: None Priority: 3 Submitted By: Bela Babik (teki321) Assigned to: Nobody (None) Summary: xrc listctrl windows custom data crash Initial Comment: I have attached the

Hello, I would like to instrument backend generated code to recover all load and store instructions' operands. What would be the best way to implement this instrumentation. I was thinking of creating a MachineFunctionPass, get all load and store instructions and add an instruction that writes the operands of Load/Store instructions to an address. How can I control the address of

Test =) soddsjbmshdwqwu -- Test, yep.

From: Gene Cumm <gene.cumm at gmail.com> COM32: Add directory functions getcwd(), opendir(), readdir() and closedir(). This depends on the patch that I just submitted creating the COMBOOT API calls. Signed-off-by: Gene Cumm <gene.cumm at gmail.com> --- The intention is to create the library calls with (hopefully) POSIX compliance such that a small application written to use these

Folks, The contour() function wants x and y to be in increasing order. I have a situation where I have a grid in x and y, and associated z values, which looks like this: x y z [1,] 0.00 20 1.000 [2,] 0.00 30 1.000 [3,] 0.00 40 1.000 [4,] 0.00 50 1.000 [5,] 0.00 60 1.000 [6,] 0.00 70 1.000 [7,] 0.00 80 0.000 [8,] 0.00 90

More scholarships are posted at http://cambodiajobs.blogspot.com ** *MBA scholarships at Queensland University of Technology* ** BGSB offers a small number of scholarships for superior applicants from *all countries* commencing in the Master of Business Administration (MBA), Master of Business Administration (Major) and Graduate Certificate in Business Administration programs. The scholarships

Hello all, I am wondering if there is a way to do a spatial error probit/logit model in R? I can't seem to find it in any of the packages. I can do it in MATLAB with Gibbs sampling, but would like to confirm the results. Ideally I would like to use this model to predict probability of parcel conversion in a future time period. This seems especially difficult in a binary outcome model