>data<-read.csv("file.csv",header=T)
比如说第2列是性别 1是男 0是女
那么你可以强制转换,也可以用函数
>sex<-factor(data[,2]) 或者 >sex<-as.factor(data[,2])
如此sex就变成了因子型的
你可以把0,1改成female ,male
>levels(sex)<-c("female","male")
最后带回去
>data[,2]<-sex
实元k_end,k_start,uz_int中的uz_int并没有在公共块中设置成real*8,而在子程序
中相应的哑元S确设置成了real*8,这可能就是不一致的原因。(不知道你把uz_int设置成什么?)
再有
kxy_integral(k_end,k_start,uz_int)
应该是
call
kxy_integral(k_end,k_start,uz_int)