similar to: Matrix formatting

I'm trying to write differential equations into R code. Can anyone help me with the following equation: (dy/dt) = v (dv/dt) = y - y^3 + c(sin t) where c is a constant. I've been using the "deriv()" function and it seems to work out the best, but I'm having trouble converting the above formula to R code. Thank you, Tim Sharac Sign up for Internet Service under $10

I'm trying to convert R code into C++, 2 conflicts have come about; a) is.missing is not a function in C++, are there any substitutes? b) need help generating random numbers bound by a normal distribution, with a mean and std dev. as inputs. Tim Sharac _________________________________________________________________

Dear All, I am given a set of files names as: velocity1.txt velocity2.txt and so on. I am sure there must be a way to read them automatically in R. It is really taking me longer to read them than to analyze them. Anybody has a suggestion to help me out with this? Many thanks Lorenzo

Hiall, I''m wondering what''s common practice when dealing with money values stored in databases. I like the idea of storing all the values as integers, i.e. all the values in eurocents not euros. Then I need to multiply all values the users enter by 100 before storing them in the database. Now the question is, what''s the most efficient and dryest way to do this?

Full_Name: Jerome Asselin Version: 1.6.2 OS: RedHat Linux 7.2 Submission from: (NULL) (142.103.173.179) Should be an easy fix... Consider the examble below: plot(0,0) legend(0,0,c("Hello!","Hi!"),pch=1:2,lty=1:2,trace=T) It gives the following trace: > plot(0,0) > legend(0,0,c("Hello!","Hi!"),pch=1:2,lty=1:2,trace=T) xchar= 0.05178 ;

Dear all, I am trying to use the Matrix package to do some calculations on rather large and sparse matrices. An example of such a matrix is given below. mig<-0.2 side<-10 np<-side^2 mig.mat<-matrix(0,np,np) diag(mig.mat[1:(np-side),(side+1):np])<-mig/4 diag(mig.mat[(side+1):np,1:(np-side)])<-mig/4 diag(mig.mat[-np,-1])<-mig/4

Jerome: No, there are no plans to JIT to SPU. That's considerably more complicated -- you'd have to figure out when to JIT to the SPU and live with all of the constraints that the SPU imposes (data reformatting, r/w DMA, ensure your code lives in 256K unless you can manage to interface with the virtual I-cache work.) Basically, it's not trivial and it doesn't quite fit into the

Hi all, As promised, I have started a howto document on setting up Samba with NT ACL support, based on experiences in my workplace, as documentation for this seems to be on the short side at the moment. It is available at the following address: http://www.bluelightning.org/linux/samba_acl_howto I have yet to convert it to proper DocBook format, and it needs more material (probably

Hi All, Currently my plot shows the y-axis in scientific notation (1e07 and so on). I want to be able to display this in dollars such that it shows $10,000,000 (including the commas). How do I do this. Also with the xlabel and ylabel. I've specified: 'title('Cash vs Time',xlab='Period',ylab=''); Im hoping that this will also not display anything on the y-axis.

Hi, I''ve some btrfs fs which are encrypted with loopaes. I decrypt them in a script during bootup. Now I want to determine if the decrypted fs is a btrfs fs. If not I missspelled the password and I can handle that error in my script. Thanks for any help! Felix -- To unsubscribe from this list: send the line "unsubscribe linux-btrfs" in the body of a message to

Full_Name: Jerome Asselin Version: 1.6.1 OS: RedHat Linux 7.2 Submission from: (NULL) (142.103.173.179) Hello, I get a segmentation fault when I run the following code. I wouldn't expect meaningful results because my response variable contains only missing values. However, I would expect something like a regular error (not a segmentation fault). library(survival) data <-

Hi list, After kudos, time comes again with problems. This time, still on the same setup as before : - Linux PDC with ldapsam, ran by RH9, OpenLDAP 2.0.27 (stock RH9 RPM+Solaris RootDSE patch), Samba 3.0.1rc1 recompiled from SRPM ; - Linux BDC is the same ; The PDC and BDC are working Ok, so I won't include the smb.conf from these. - Solaris 9 domain member (jersey) gets Posix accounts

Hello im running sage line50 version 11 on a freebsd machine with samba version 2x My smb.conf for share looks like this [sage] path = /40gb/sage writeable = yes browseable = no create mask = 0777 directory mask = 0777 veto oplock files = /*.mdb/*.MDB/*.ldb/*.LDB locking = yes share modes = yes The problem is when two users for example try and update customer records one

Hello, I have a beginner's Q: if i want to know all the functions provided by a package, what is command for that? in another word, Is there a command to list all the commands available in a packege? thanks a lot, yan

Dear R users, I am interested in uni- bi- multimodality tests, for analysing reaction times data. I was lead to Hartigan's dip test (Ann. Statistics, 13, 1985, pp. 70-84, Applied Statistics, 34, 1985, 320-325). Not being a programmer I am unable to translate the Fortran code given in ref. 2 into a R function. I'd be glad to learn that someone already did it, or has devised a better

I ran this script in a source file on my Mac: library(lattice) year <- 1900:2000 dollars <- (year-1899)^2 plot(year,dollars) quartz("year") histogram(~dollars) The first plot appears in Quartz 2. The second quartz window, named year, opens but the histogram doesn't appear. However, when I copy and paste this script directly into the R console, both quartz windows (Quartz 2

Thank very much for this answer, so my last question will be: is it possible to use the LLVM JIT on a PS3 with Yellow Dog 6.2 distribution, instead of the LLVM interpreter, by using the PPE as it seems to be similar to 64-bit PowerPC processors? 2009/12/18 Scott Michel <scooter.phd at gmail.com> > Jerome: > > No, there are no plans to JIT to SPU. That's considerably more

When trying to browser users or groups on the server I see these messages in the log file. [2004/09/02 10:40:15, 0] lib/smbldap.c:smbldap_search_suffix(1101) smbldap_search_suffix: Problem during the LDAP search: (Size limit exceeded) [2004/09/02 10:40:15, 0] passdb/pdb_ldap.c:ldapsam_setsampwent(1173) ldapsam_setsampwent: LDAP search failed: Size limit exceeded We are unable to browse

Hi, we have tested Samba 3.0.0 and 3.0.1 with LDAP-Support (--with-ldap) on Solaris 8 and it worked fine. The machine authenticates against an OpenLDAP server. Patch 108993-23 is applied and we use native Sun LDAP client modules. On Solaris 9 we ran into problems with secondary groups. Users cannot access files if the rights are based on a secondary group and if this information is stored on

Hi all, Ive got a database with 10 columns (different variables) for 100 subjects, each column with different # of NA's. I'd like to know if it is possible to use a function to exclude the NA's using only a specific column, lets say: Data2 <- omit.exclude(Data1$column1) ??, then Data3 <- omit.exclude(Data1$column2) and so on I tried the code above but with no results Thanks