search for: bremson

...a function matSummary.spss. I think I ought to be able to call >matSummary(ret) to run the function, but only >matSummary.spss(ret) will work. What am I doing wrong here? This seems like a simple problem yet I've been able to find nothing in the archives about this. Joel -- Joel Bremson Graduate Student Institute for Transportation Studies - UC Davis http://etrans.blogspot.com [[alternative HTML version deleted]]

...r :: n n = 10 allocate(a(n,n+1)) end I then compiled: >g95 -c allloc.f >R CMD SHLIB allloc.o Here is what happens when I try to dyn.load it in R 2.1.0 > dyn.load("allloc.so") Error in dyn.load(x, as.logical(local), as.logical(now)) : unable to load shared library '/home/jbremson/dev/f_code/allloc.so': /home/jbremson/dev/f_code/allloc.so: undefined symbol: _g95_filename Is there any way to deal with this? Is it somehow possible to reach the fortran 95 through C and avoid this whole problem? Regards, Joel Bremson UC Davis [[alternative HTML version deleted]]

...r :: n n = 10 allocate(a(n,n+1)) end I then compiled: >g95 -c allloc.f >R CMD SHLIB allloc.o Here is what happens when I try to dyn.load it in R 2.1.0 > dyn.load("allloc.so") Error in dyn.load(x, as.logical(local), as.logical(now)) : unable to load shared library '/home/jbremson/dev/f_code/allloc.so': /home/jbremson/dev/f_code/allloc.so: undefined symbol: _g95_filename Is there any way to deal with this? Is it somehow possible to reach the fortran 95 through C and avoid this whole problem? Regards, Joel Bremson UC Davis [[alternative HTML version deleted]]

...w up to my question from yesterday. I want to do in R what is called a "3d scatter plot with drop lines" in S-PLUS. Basically, it's a 3dscatterplot with lines connecting the x-y grid to the z points. The lines give a better perspective on the shape of the data surface. How to? Joel Bremson UC Davis Statistics

I've got a 4x4 matrix of points from a 2-way ANOVA I'd like to plot. The x,y correspond to the treatment groups and look like this ((1,1),(1,2),(1,3),(1,4),(2,1),...). The z is the 4x4 matrix. How can I get persp to grid the x,y axis with only the numbers 1-4 on both? Regards, Joel Bremson UC Davis Statistics Dept.

...uppress -L/sw/lib -L/usr/local/lib -o $@ $^ \ -L$(RLIB_LOC)/lib -lR ###EOF#### The -L lib dirs are not correct. On a *nix platform I would do something like this sh -x R CMD SHLIB ... to get at the R internal link information but I can't get that to work on Cygwin. Regards, Joel -- Joel Bremson Graduate Student Institute for Transportation Studies - UC Davis http://etrans.blogspot.com [[alternative HTML version deleted]]

...al data looks like this: 1 2 3 1 1 1 2 2 2 3 3 I can easily invert the output mathematically with: q = sapply(m,function(x){ x + 2*(median(unique(m))-x)}) (m is composed of sequential integers starting at one) ,but it seems as though something wrong is happening with read.spss. Any ideas? Joel Bremson Graduate Student UC Davis [[alternative HTML version deleted]]

...es as follows in bash: export PATH=~/bin/:$PATH export F77=gfortran export LD_LIBRARY_PATH=~/bin/irun/lib export GFORTRAN_STDIN_UNIT=-1 Now I'm trying to write a Makevars file for my test package and not quite sure how to do it. I've tried FF=gfortran GFORTRAN_STDIN_UNIT=-1 FLIBS=/home/jbremson/bin/irun/lib/libgfortran.a but when running an R CMD check on the package I see that it's still using g77 to compile: -----output---- ... ** libs g77 -fPIC -g -O2 -c estimate.f -o estimate.o estimate.f: In subroutine `estimate': estimate.f:20: forall (i = 1:nxrows) beta(i) = i * 2 ^ Inva...

...off by running the full rtest data set I also don't have any problems with R. This package runs without problems on a G4/Tiger OS X mac with R-2.1.1, although we do compile it with gfortran on that platform. This is just a report. If anyone has any questions let me know. Best, Joel -- Joel Bremson Graduate Student Institute for Transportation Studies - UC Davis http://etrans.blogspot.com

...ould get a core dump in order to get some more clues about what is going on, but no luck. I'm relatively new to fortran programming and have never used gdb before, although I'm familiar with other debuggers. Any tips on how to proceed from this point would be appreciated. Regards, Joel Bremson Statistics Graduate Student UC Davis [[alternative HTML version deleted]]

...suppress -L/sw/lib -L/usr/local/lib -o $@ $^ \ -L$(RLIB_LOC)/lib -lR ###EOF#### The -L lib dirs are not correct. On a *nix platform I would do something like this sh -x R CMD SHLIB ... to get at the R internal link information but I can't get that to work on Cygwin. Regards, Joel -- Joel Bremson Graduate Student Institute for Transportation Studies - UC Davis http://etrans.blogspot.com [[alternative HTML version deleted]]

Hi all, If I write v = vector(mode="numeric",length=10) I'm still allowed to assign non-numerics to v. Furthermore, R figures out what kind of vector I've got anyway when I use the mode() function. So what is it that assigning a mode does? Joel

...FORTRAN_FILES:%.f=%.o) g95 -L/Library/Frameworks -o $@ $^ -L/Library/Frameworks/R.framewo\ rk/Resources -----end makefile------------------------ Here is the relevant output of R CMD CHECK: * checking for working latex ...sh: line 1: latex: command not found NO * using log directory '/Users/jbremson/dev/bpkg.Rcheck' * using R version 2.1.0, 2005-05-12 * checking for file 'bpkg/DESCRIPTION' ... OK * this is package 'bpkg' version '1.0-1' * checking if this is a source package ... OK * Installing *source* package 'bpkg' ... ** libs g95 -c -g dbest_dbase_class...

...lled a "3d scatter plot with drop lines" in S-PLUS. >> >>Basically, it's a 3dscatterplot with lines connecting the x-y grid to >>the z points. >>The lines give a better perspective on the shape of the data surface. >> >>How to? >> >>Joel Bremson >>UC Davis Statistics >> >>______________________________________________ >>R-help@stat.math.ethz.ch mailing list >>https://stat.ethz.ch/mailman/listinfo/r-help >>PLEASE do read the posting guide! http://www.R-project.org/posting-guide.html > > >-- >...

...Fortran("add",as.numeric(a),as.numeric(b),c=as.numeric(0)) return(out$c) } > x = myadd(3,7) [1] 10 And you've done it! Be vanquished naysayers, simple addition will trouble us no more! (I thought a very simple tutorial on integrating Fortran code into R might be of use.) Joel Bremson jbremson@wald.ucdavis.edu Statistics Graduate Student UC Davis [[alternative HTML version deleted]]